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

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

  • добавлена проволока на лопатке (ребро жесткости);

  • в расчете используются два материала, сталь (диск и проволока) и лопатка (алюминий);

  • в расчете используются угловая скорость и температурные нагрузки (200 °C на диске и 300 °C на лопатке и проволоке);

  • ПГУ задано на поверхностях, которые не находятся в одной плоскости, а также с неконформной сеткой.

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

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

Или введи в командну строку команды:

reset
#{h=0.3}
create surface rectangle width {h} height 0.2 zplane 
create vertex 0.07428 0.01871 0 on surface 1 
create vertex 0.000169 0.007849 0 on surface 1 
create vertex -0.065 0.01939 6.939e-18 on surface 1 
create vertex 0.07514 -0.00287 0 on surface 1 
create vertex 0.0006699 -0.02354 0 on surface 1 
create vertex -0.06414 -0.001799 0 on surface 1 
create vertex 0.09311 0.02074 0 on surface 1 
create vertex -0.08031 0.02169 0 on surface 1 
create vertex 0.09314 0.01093 0 on surface 1 
create vertex -0.08028 0.01345 0 on surface 1 
create curve spline location vertex 8 9 10  
create curve spline location vertex 5 6 7  
create curve spline location vertex 8 13 11 5  
create curve spline location vertex 10 14 12 7  
create surface curve 5 7 6 8  
Surface 2  copy move z 0.25 repeat 2 
rotate Surface 3  angle 15  about Z include_merged
rotate Surface 4  angle 35  about Z include_merged 
Surface 3  scale X 1.1 Y 0.9 Z 1.0 
Surface 4  scale X 1.3 Y 0.9 Z 1.0 
create curve vertex 5 20  
create curve vertex 28 24  
create curve vertex 7 21  
create curve vertex 32 25  
create curve vertex 10 22  
create curve vertex 36 26  
create curve vertex 8 19  
create curve vertex 40 23  
create surface curve 6 17 10 19  
create surface curve 7 23 9 17  
create surface curve 23 5 21 12  
create surface curve 19 11 21 8  
create surface curve 9 24 13 18  
create surface curve 22 15 20 35  
create surface curve 25 18 14 20  
create surface curve 22 16 24 31  
create volume surface 5 6 9 11 10 8 7 12 4 2  heal 
move Volume 2  z 0.2 include_merged 
create Cylinder height 0.2 radius 0.21 
rotate Volume 13  angle 90  about Y include_merged 
move Volume 2  x -0.01 include_merged 
create Cylinder height 0.2 radius 0.11 
rotate Volume 14  angle 90  about Y include_merged 
subtract volume 14 from volume 13  
create curve vertex 49 27  
create vertex -0.1 0.03752 0.2066 on curve 49 
create vertex -0.1002 -0.05317 0.2028 on curve 49 
create vertex -0.1 0.03032 0.1837 on surface 20 
create vertex -0.1 -0.04738 0.1819 on surface 20 
create vertex -0.1 0.01485 0.1763 on surface 20 
create vertex -0.1 -0.03257 0.1752 on surface 20 
create vertex -0.1 0.01478 0.1689 on surface 20 
create vertex -0.1 -0.0314 0.1681 on surface 20 
create vertex -0.1 0.02495 0.1623 on surface 20 
create vertex -0.1 -0.04014 0.1613 on surface 20 
create vertex -0.1 0.02489 0.1561 on surface 20 
create vertex -0.1 -0.03865 0.1551 on surface 20 
create vertex -0.1 0.002773 0.1439 on surface 20 
create vertex -0.1 -0.01179 0.1437 on surface 20 
create curve vertex 73 75  
create curve vertex 75 77  
create curve vertex 77 79  
create curve vertex 79 81  
create curve vertex 81 83  
create curve vertex 83 85  
create curve vertex 85 86  
create curve vertex 86 84  
create curve vertex 84 82  
create curve vertex 82 80  
create curve vertex 80 78  
create curve vertex 78 76  
create curve vertex 76 74  
sweep curve 56 57 58 59 60 61 62 63 64 65 66 67 68  along curve 55  
webcut body 13  with sheet extended from surface 22 23 24 25 26 27 28 29 30 31 32 33 34  
rotate Surface all angle 20  about X include_merged 
webcut body 28 with plane yplane rotate 10 about z
rotate Surface all angle -36  about X include_merged 
webcut body 28  with plane yplane rotate 10 about z 
delete volume 29 28  
delete surface 1 to 12, 22 to 34
webcut volume 2  with sheet extended from surface 49  
delete volume 31
rotate Volume all angle 15  about X include_merged 
delete Curve 55
rotate Volume all angle 90  about Y include_merged 

Достраивание модели

Создайте круг радиуса 0.015. Для этого на панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Создать. Из выплывающего списка выберите Круг и задайте необходимые параметры. Нажмите Применить.

Переместите созданную поверхность. На панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Преобразовать, из выпадающего списка выберите Переместить. Задайте необходимые параметры. Нажмите Применить.

Создайте объем при помощи вытягивание поверхностей. Для этого на панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Создать. Из списка геометрических примитивов выберите Протяжка и задайте требуемые параметры. Нажмите Применить.

Проделайте то же самое для создания второго объема.

Объедините два объема. Для этого на панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Логические операции. Из списка геометрических примитивов выберите Объединить и задайте требуемые параметры. Нажмите Применить.

Удалите лишнюю поверхность. На панели команд выберите Режим - Геометрия, Объект - Поверхность, Действие - Удалить. Введи ID поверхности и нажмите Применить

Вычтите один объем из другого. Для этого на панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Логические операции. Из списка геометрических примитивов выберите Вычесть и задайте требуемые параметры. Нажмите Применить.

Сделайте разрез. На панели команд выберите Режим - Геометрия, Объект - Объем, Действие - Разрез. Из выпадающего списка выберите Радиус цилиндра. Задайте необходимые параметры и нажмите Применить.

Проделайте то же самое для радиуса 0.63.

Построение сетки

Задайте схему для построения четырехугольной сетки. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Преимущественно четырехугольная и укажите необходимые параметры. Нажмите Применить схему.

Постройте сетку на поверхности. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и укажите необходимые параметры. Нажмите Построить сетку.

Задайте схему для построения четырехугольной сетки. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Преимущественно четырехугольная и укажите необходимые параметры. Нажмите Применить схему.

Постройте сетку на поверхности. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и укажите необходимые параметры. Нажмите Построить сетку.

Задайте схему для построения четырехугольной сетки. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Преимущественно четырехугольная и укажите необходимые параметры. Нажмите Применить схему.

Постройте сетку на поверхности. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и укажите необходимые параметры. Нажмите Построить сетку.

Задайте схему для построения четырехугольной сетки. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Преимущественно четырехугольная и укажите необходимые параметры. Нажмите Применить схему.

Постройте сетку на поверхности. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и укажите необходимые параметры. Нажмите Построить сетку.

Задайте разбиение на кривых. На панели команд выберите модуль построения сетки на кривых Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выпадающего списка выберите Равномерно. Укажите необходмые значения и нажмите Применить.

Задайте разбиение на кривых. На панели команд выберите модуль построения сетки на кривых Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выпадающего списка выберите Равномерно. Укажите необходмые значения и нажмите Применить.

Задайте разбиение на кривых. На панели команд выберите модуль построения сетки на кривых Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выпадающего списка выберите Равномерно. Укажите необходмые значения и нажмите Применить.

Задайте разбиение на кривых. На панели команд выберите модуль построения сетки на кривых Режим - Сетка, Объект - На кривых, Действие - Построение сетки. Из выпадающего списка выберите Равномерно. Укажите необходмые значения и нажмите Применить.

Задайте схему для построения четырехугольной сетки. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Преимущественно четырехугольная и укажите необходимые параметры. Нажмите Применить схему.

Постройте сетку на поверхности. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и укажите необходимые параметры. Нажмите Построить сетку.

Задайте схему для построения четырехугольной сетки. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Построение сетки. Из выплывающего списка выберите Преимущественно четырехугольная и укажите необходимые параметры. Нажмите Применить схему.

Постройте сетку на поверхности. На панели команд выберите Режим - Сетка, Объект - Поверхностная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и укажите необходимые параметры. Нажмите Построить сетку.

Задайте разбиение сетки. На панели команд выберите Режим - Сетка, Объект - Объемная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и укажите необходимые параметры. Нажмите Задать размер.

Постройте тетраэдральную сетку. На панели команд выберите Режим - Сетка, Объект - Объемная, Действие - Построение сетки. Из выплывающего списка выберите Тетраэдральная и укажите необходимые параметры. Нажмите Построить сетку.

Выполните протяжку сетки. На панели команд выберите Режим - Сетка, Объект - Объемная, Действие - Построение сетки. Из выпадающего списка выберите Протяжка и задайте необходимые параметры. Нажмите Применить схему, Построить сетку.

Проделайте то же самое для объема 35, исходной поверхности 90 и целевой поверхности 115.

Проделайте то же самое для объема 2, исходной поверхности 120 и целевой поверхности 4.

Проделайте то же самое для объема 33, исходной поверхности 102 и целевой поверхности 105.

Задайте схему построения сетки. На панели команд выберите Режим - Сетка, Объект - Объемная, Действие - Построение сетки. Из выплывающего списка выберите Тетраэдральная и укажите необходимые параметры. Нажмите Применить схему.

Постройте сетку. На панели команд выберите Режим - Сетка, Объект - Объемная, Действие - Интервалы. Из выплывающего списка выберите Автоматический размер и укажите необходимые параметры. Нажмите Построить сетку.

Задание материала и свойств блока

Создайте два материала. На панели команд выберите Режим - Материал, Объект - Управление материалами. Перетащите материал «Углеродистая сталь» в поле Материал. Далее создайте второй материал, для этого в средней колонке укажите имя втогорого материала. В колонке свойств откройте список Упругость и перетащите Материал Гука в колонку Свойства материала, затем в левой колонке перейдите в раздел и перетащите Общие → Плотность, Температурные свойства → Температурные изотропные. Задайте необходимые значения. Нажмите Применить.

Создайте первый блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте необходимые параметры. Нажмите Применить.

Задайте свойства первому блоку. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте необходимые параметры. Нажмите Применить.

Создайте второй блок. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Добавить сущность в блок. Задайте необходимые параметры. Нажмите Применить.

Задайте свойства второму блоку. На панели команд выберите Режим - Блоки, Объект - Блок, Действие - Свойства/параметры блока. Задайте необходимые параметры. Нажмите Применить.

Задание граничных условий

Задайте связанный контакт между поверхностями. На панели команд выберите Режим - Граничные условия, Объект - Контакт, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Задайте связанный контакт между поверхностями. На панели команд выберите Режим - Граничные условия, Объект - Контакт, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Задайте связанный контакт между поверхностями. На панели команд выберите Режим - Граничные условия, Объект - Контакт, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Задайте связанный контакт между поверхностями. На панели команд выберите Режим - Граничные условия, Объект - Контакт, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Задайте связанный контакт между поверхностями. На панели команд выберите Режим - Граничные условия, Объект - Контакт, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Закрепите один конец лопатки по всем перемещениям. На панели команд выберите Режим - Граничные условия, Объект - Перемещение, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Приложите угловую скорость. На панели команд выберите Режим - Граничные условия, Объект - Угловая скорость, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Задайте значение температуры. На панели команд выберите Режим - Граничные условия, Объект - Температура, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Задайте значение температуры. На панели команд выберите Режим - Граничные условия, Объект - Температура, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Задайте периодическое условие. На панели команд выберите Режим - Граничные условия, Объект - Периодическое условие, Действие - Создать. Задайте необходимые параметры. Нажмите Применить.

Запуск расчета

Для запуска на расчет в CAE Fidesys выберите на панели команд Режим - Настройки расчета, Настройки расчета - Модальный анализ, Модальный анализ - Общие. Задайте требуемые параметры и нажмите Применить.

Задайте параметры решателя. На панели команд выберите Режим - Настройки расчета, Настройки расчета - Модальный анализ, Модальный анализ - Решатель. Задайте требуемые параметры и нажмите Применить.

Задайте параметры решателя СЧ. На панели команд выберите Режим - Настройки расчета, Настройки расчета - Модальный анализ, Модальный анализ - Решатель СЧ. Задайте требуемые параметры и нажмите Применить, Начать расчет.

Анализ результатов

Собственный частоты выводятся в Командную строку:

Откройте файл с результатами. Это можно сделать тремя способами:

  • Нажмите Ctrl+E;

  • В главном меню выберите Расчёт - Результаты. Нажмите Открыть последний результат;

  • На панели команд выберите Результаты (Режим - Результаты, Результаты - Открыть Результаты).

Для просмотра результатов с преднагружением откройте файл с пометкой в названии preload.

На панели инструментов выберите Фильтры - Алфавитный указатель - Трансформация. Во вкладке Свойства укажите угол поворта модели. Нажмите Применить.

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

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

Отобразите распределение напряжений. Для этого на верхней панели выберите из первого выпадающего списка выберите Напряжения, из второго – Мизес.

Использование консольного интерфейса

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

reset
#{h=0.3}
create surface rectangle width {h} height 0.2 zplane 
create vertex 0.07428 0.01871 0 on surface 1 
create vertex 0.000169 0.007849 0 on surface 1 
create vertex -0.065 0.01939 6.939e-18 on surface 1 
create vertex 0.07514 -0.00287 0 on surface 1 
create vertex 0.0006699 -0.02354 0 on surface 1 
create vertex -0.06414 -0.001799 0 on surface 1 
create vertex 0.09311 0.02074 0 on surface 1 
create vertex -0.08031 0.02169 0 on surface 1 
create vertex 0.09314 0.01093 0 on surface 1 
create vertex -0.08028 0.01345 0 on surface 1 
create curve spline location vertex 8 9 10  
create curve spline location vertex 5 6 7  
create curve spline location vertex 8 13 11 5  
create curve spline location vertex 10 14 12 7  
create surface curve 5 7 6 8  
Surface 2  copy move z 0.25 repeat 2 
rotate Surface 3  angle 15  about Z include_merged
rotate Surface 4  angle 35  about Z include_merged 
Surface 3  scale X 1.1 Y 0.9 Z 1.0 
Surface 4  scale X 1.3 Y 0.9 Z 1.0 
create curve vertex 5 20  
create curve vertex 28 24  
create curve vertex 7 21  
create curve vertex 32 25  
create curve vertex 10 22  
create curve vertex 36 26  
create curve vertex 8 19  
create curve vertex 40 23  
create surface curve 6 17 10 19  
create surface curve 7 23 9 17  
create surface curve 23 5 21 12  
create surface curve 19 11 21 8  
create surface curve 9 24 13 18  
create surface curve 22 15 20 35  
create surface curve 25 18 14 20  
create surface curve 22 16 24 31  
create volume surface 5 6 9 11 10 8 7 12 4 2  heal 
move Volume 2  z 0.2 include_merged 
create Cylinder height 0.2 radius 0.21 
rotate Volume 13  angle 90  about Y include_merged 
move Volume 2  x -0.01 include_merged 
create Cylinder height 0.2 radius 0.11 
rotate Volume 14  angle 90  about Y include_merged 
subtract volume 14 from volume 13  
create curve vertex 49 27  
create vertex -0.1 0.03752 0.2066 on curve 49 
create vertex -0.1002 -0.05317 0.2028 on curve 49 
create vertex -0.1 0.03032 0.1837 on surface 20 
create vertex -0.1 -0.04738 0.1819 on surface 20 
create vertex -0.1 0.01485 0.1763 on surface 20 
create vertex -0.1 -0.03257 0.1752 on surface 20 
create vertex -0.1 0.01478 0.1689 on surface 20 
create vertex -0.1 -0.0314 0.1681 on surface 20 
create vertex -0.1 0.02495 0.1623 on surface 20 
create vertex -0.1 -0.04014 0.1613 on surface 20 
create vertex -0.1 0.02489 0.1561 on surface 20 
create vertex -0.1 -0.03865 0.1551 on surface 20 
create vertex -0.1 0.002773 0.1439 on surface 20 
create vertex -0.1 -0.01179 0.1437 on surface 20 
create curve vertex 73 75  
create curve vertex 75 77  
create curve vertex 77 79  
create curve vertex 79 81  
create curve vertex 81 83  
create curve vertex 83 85  
create curve vertex 85 86  
create curve vertex 86 84  
create curve vertex 84 82  
create curve vertex 82 80  
create curve vertex 80 78  
create curve vertex 78 76  
create curve vertex 76 74  
sweep curve 56 57 58 59 60 61 62 63 64 65 66 67 68  along curve 55  
webcut body 13  with sheet extended from surface 22 23 24 25 26 27 28 29 30 31 32 33 34  
rotate Surface all angle 20  about X include_merged 
webcut body 28 with plane yplane rotate 10 about z
rotate Surface all angle -36  about X include_merged 
webcut body 28  with plane yplane rotate 10 about z 
delete volume 29 28  
delete surface 1 to 12, 22 to 34
webcut volume 2  with sheet extended from surface 49  
delete volume 31
rotate Volume all angle 15  about X include_merged 
delete Curve 55
rotate Volume all angle 90  about Y include_merged 
create surface circle radius .015 yplane 
move Surface 100 x .6 include_merged 
sweep surface 100 zaxis angle 15  keep 
sweep surface 100 zaxis angle -21  keep 
unite volume 33 34  
delete surface 100
subtract volume 33 from volume 2 keep_tool
webcut volume 2  with cylinder radius .57 axis z center 0 0 0  
webcut volume 2  with cylinder radius .63 axis z center 0 0 0  
surface 48 scheme quad_dominant 
surface 48 size auto factor 4
mesh surface 48
surface 90 scheme quad_dominant 
surface 90 size auto factor 4
mesh surface 90
surface 120 scheme quad_dominant 
surface 120 size auto factor 4
mesh surface 120
surface 102 scheme quad_dominant 
surface 102 size auto factor 7
mesh surface 102
curve 111 241  interval 25 ###195
curve 111 241  scheme equal
mesh curve 111 241 
curve 273  interval 15
curve 273  scheme equal
mesh curve 273 
curve 282  interval 10
curve 282  scheme equal
mesh curve 282
curve 252  interval 50
curve 252  scheme equal
mesh curve 252 
surface 67  scheme quad_dominant 
surface 67  size auto factor 2
mesh surface 67 
surface 85  scheme quad_dominant 
surface 85  size auto factor 3
mesh surface 85 
volume 30 size auto factor 2
volume 30 scheme tetmesh
mesh volume 30
volume 13 redistribute nodes off 
volume 13 scheme Sweep  source surface 48   target surface 50  sweep transform least squares 
volume 13 autosmooth target on  fixed imprints off  smart smooth off 
volume 13 redistribute nodes off 
volume 13 scheme Sweep  source surface 48   target surface 50  sweep transform least squares 
volume 13 autosmooth target on  fixed imprints off  smart smooth off 
mesh volume 13
volume 35 redistribute nodes off 
volume 35 scheme Sweep  source surface 90   target surface 115  sweep transform least squares 
volume 35 autosmooth target on  fixed imprints off  smart smooth off 
volume 35 redistribute nodes off 
volume 35 scheme Sweep  source surface 90   target surface 115  sweep transform least squares 
volume 35 autosmooth target on  fixed imprints off  smart smooth off 
mesh volume 35
volume 2 redistribute nodes off 
volume 2 scheme Sweep  source surface 120   target surface 4  sweep transform least squares 
volume 2 autosmooth target on  fixed imprints off  smart smooth off 
volume 2 redistribute nodes off 
volume 2 scheme Sweep  source surface 120   target surface 4  sweep transform least squares 
volume 2 autosmooth target on  fixed imprints off  smart smooth off 
mesh volume 2
volume 33 redistribute nodes off 
volume 33 scheme Sweep  source surface 102   target surface 105  sweep transform least squares 
volume 33 autosmooth target on  fixed imprints off  smart smooth off 
volume 33 redistribute nodes off 
volume 33 scheme Sweep  source surface 102   target surface 105  sweep transform least squares 
volume 33 autosmooth target on  fixed imprints off  smart smooth off 
mesh volume 33
volume 36 scheme tetmesh
undo group begin
volume 36 size auto factor 5
mesh volume 36
undo group end
create material 1 from 'Углеродистая сталь'
create material 2
modify material 2 name 'Алюминий'
modify material 2 set property 'MODULUS' value 7e+10
modify material 2 set property 'POISSON' value 0.33
modify material 2 set property 'DENSITY' value 2640
modify material 2 set property 'ISO_THERMAL_EXPANSION' value 2.47e-05
modify material 2 set property 'ISO_CONDUCTIVITY' value 122
set duplicate block elements off
block 1 add volume 30 33
block 1 name 'Диск_ребро'
block 1 material 1 cs 1 element solid order 1
set duplicate block elements off
block 2 add volume 13 2 35 36
block 2 name 'Лопатка'
block 2 material 2 cs 1 element solid order 1
create contact master surface 63 62 59 58 57 61 56 55 53 52 51 54 60 slave surface 35 36 37 38 39 40 41 42 43 44 45 46 47 type tied offset 0.0 tolerance 0.0005 method mpc
create contact master surface 49 slave surface 90 type tied offset 0.0 tolerance 0.0005 method auto
create contact master surface 115 slave surface 110 type tied offset 0.0 tolerance 0.0005 method auto
create contact master surface 120 slave surface 125 type tied offset 0.0 tolerance 0.0005 method auto
create contact master surface 101 slave surface 107 type tied offset 0.0 tolerance 0.0005 method auto
create displacement  on surface 86 dof all fix  
create angular velocity global
modify angular velocity 1 dof 3 value 3.14
create temperature  on volume 30 value 200
create temperature  on volume 13 35 36 2 33 value 300
create periodic master surface 67 102 slave surface 85 105 general cs 1 sectors 10
analysis type eigenfrequencies elasticity dim3 preload on
eigenvalue find 10 smallest
solver method direct use_uzawa auto try_other off
eigenvalue eigensolver eps auto linearsolver direct