Презентация на тему MSC.Nastran 102 2001 - 03

КОЛЕБАНИЙ… 3 - 3
ОБЗОР ТЕОРИИ……………………………………………………………………………… 3 - 4
ОСНОВНЫЕ

СВОЙСТВА СОБСТВЕННЫХ КОЛЕБАНИЙ…..………………………. 3 - 8
ДОПОЛНИТЕЛЬНЫЕ СВОЙСТВА ФОРМ СОБСТВЕННЫХ КОЛЕБАНИЙ.……… 3 - 12
МЕТОДЫ ВЫЧИСЛЕНИЙ………….………………………………………………………. 3 - 14
ТЕОРИЯ МЕТОДА Штурма……………………………………………………………………3 - 16
МЕТОД Ланцоша………….…………………………………………………………………… 3 - 17
ИНТЕРФЕЙС ПОЛЬЗОВАТЕЛЯ ДЛЯ МЕТОДА Ланцоша.…………………………….. 3 - 18
ИНТЕРФЕЙС ПОЛЬЗОВАТЕЛЯ ДЛЯ ДРУГИХ МЕТОДОВ……………………………… 3 - 19
УПРАВЛЕНИЕ РЕШЕНИЕМ ПРИ АНАЛИЗЕ СОБСТВЕННЫХ КОЛЕБАНИЙ……… 3 - 21
ВИДЫ ВЫЧИСЛЯЕМЫХ ВЕЛИЧИН..……………………………………………………… 3 - 22
ПРИМЕР №1 – МОДАЛЬНЫЙ АНАЛИЗ ПЛОСКОЙ ПЛАСТИНЫ.…………………… 3 - 23
“ГЕОМЕТРИЯ” ПЛАСТИНЫ……………………….………………………………………… 3 - 26
ВХОДНОЙ ФАЙЛ ДЛЯ ПРИМЕРА №1……..……………………………………………… 3 - 29
РЕЗУЛЬТАТЫ РЕШЕНИЯ ПРИМЕРА №1..……………………………………………….. 3 - 30

Например, если ротационную машину предполагается установить на некоторое основание,

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

(1)
Представим решение в форме
(2)
(Физически это означает, что все координаты изменяются синхронно – форма деформаций системы остается постоянной, меняется только ее амплитуда.)
Из уравнения (2): (3)

Подставляя уравнения (2) и (3) в уравнение (1), имеем


что (после деления на eiωt) упрощается к виду

Это – формулировка задачи о собственных значениях.

, то из уравнения (4):

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

или

где

“присоединенными” массами, тогда будет N

которые являются решениями собственной задачи.. Эти ω’s(ω1,ω2,…,ωn) – частоты собственных колебаний, характеристические частоты, фундаментальные частоты или резонансные частоты.
Собственный вектор , ассоциирующийся с частотой , называется нормальной модой или модальной формой (формой собственных колебаний). Нормальная мода соответствует некоторой форме деформаций конструкции.
При колебаниях форма конструкции в любой момент времени является линейной комбинацией её нормальных мод.

и действительные (что характерно для конечно-элементных моделей обычных конструкций),

справедливы следующие условия ортогональности:

и

а также


Собственные частоты (ω1,ω2,…) измеряются в рад/c. Для этого можно использовать также Гц (колебание/с), причем

колебаний жесткого тела.




Если конструкция не полностью закреплена, т.е. существует

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

Например,

являются одними и теми же модами колебаний.

,

с помощью выбранного метода. В MSC.Nastran предусмотрены три метода

нормализации (исключая метод Ланцоша):
Единичная обобщенная масса (по умолчанию)


Единичное значение наибольшего для каждой моды компонента набора A-set
Единичное значение для указанного компонента (не рекомендуется)
При использовании метода Ланцоша возможна нормализация по методу единичной обобщенной массы и методу единичного наибольшего компонента.

силы и напряжения в них зависят от величины деформации

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

Соотношение деформация - напряжение

Соотношение перемещение - сила

Соотношение перемещения – энергия деформации

деформации

Модальные напряжения

Модальные силы

Модальную энергию деформации


Дополнительные модальные характеристики могут быть

получены в табличной форме путем задания соответствующих параметров в Case Control Section или в матричной форме путем использования модуля DRMS1 на языке DMAP.

собственных значений:
Последовательные (tracking) методы (см. Приложение B)
Собственные значения (или

собственные частоты) определяются поочередно методом итераций. Возможен выбор между двумя вариантами метода “inverse power method”: INV и SINV. Эти методы более приемлемы при необходимости определения небольшого количества собственных частот. В целом, метод SINV более надежен, чем метод INV.
Методы трансформации (см. Приложение B)
Собственная задача

преобразуется в форму

где

с использованием метода Гивенса (Givens) или метода Хаусхолдера (Householder).

После этого все собственные значения определяются одновременно с помощью алгоритма QR. Предусматривается по два варианта методов Гивенса и Хаусхолдера: GIV, MGIV, HOU и MHOU. Эти методы более эффективны, если необходимо определить большое количество собственных значений.
Метод Ланцоша (рекомендуемый метод)
Этот метод является комбинированным (tracking-transformation) методом.

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

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

метода Штурма
Применение метода Гивенса и алгоритма QL
Возможность использования для модального анализа и анализа устойчивости
Нормализация только по массе и наибольшему компоненту
В отличие от метода GDR (см. Приложение A), нет необходимости в скалярных переменных (точках) типа QSET, ASET и т.п., (если только не выполняется анализ методом модального синтеза)

103
Case Control Section
METHOD (инициализация оператора EIGRL)
Bulk Data Section
EIGRL (метод

Ланцоша)

элементов
ELSTRESS (или STRESS)
ESE
EKE
ELFORCE (или FORCE)
STRAIN
Специальный оператор
MODES – результаты для

каждой моды выводятся как результаты
для отдельного “случая” расчета

найти первые десять собственных частот и форм колебаний плоской

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

(DFMSYN)
DECOMP ORDERING METHOD CHOSEN: DEFAULT, ORDERING METHOD

USED: MMD
*** USER INFORMATION MESSAGE 5010 (LNCILD)
STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION.
TRIAL EIGENVALUE = 9.433958D+07, CYCLES = 1.545849D+03 NUMBER OF EIGENVALUES BELOW THIS VALUE = 3
*** USER INFORMATION MESSAGE 5010 (LNCILD)
STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION.
TRIAL EIGENVALUE = 1.912348D+09, CYCLES = 6.959910D+03 NUMBER OF EIGENVALUES BELOW THIS VALUE = 10
TABLE OF SHIFTS: (LNNRIGL)
SHIFT # SHIFT VALUE FREQUENCY, CYCLES # EIGENVALUES BELOW # NEW EIGENVALUES FOUND
1. 9.4339576E+07 1.5458490E+03 3 10
2. 1.9123482E+09 6.9599097E+03 10 0
1 NORMAL MODES EXAMPLE APRIL 8, 1998 MSC.Nastran 4/6/98 PAGE 6

0 SUBCASE 1
0
E I G E N V A L U E A N A L Y S I S S U M M A R Y (READ MODULE)
BLOCK SIZE USED ...................... 7
NUMBER OF DECOMPOSITIONS ............. 2
NUMBER OF ROOTS FOUND ................ 10
NUMBER OF SOLVES REQUIRED ............ 7
1 NORMAL MODES EXAMPLE APRIL 8, 1998 MSC.Nastran 4/6/98 PAGE 7

0 SUBCASE 1

R E A L E I G E N V A L U E S
MODE EXTRACTION EIGENVALUE RADIANS CYCLES GENERALIZED GENERALIZED
NO. ORDER MASS STIFFNESS
1 1 7.056994E+05 8.400591E+02 1.336996E+02 1.000000E+00 7.056994E+05
2 2 1.878432E+07 4.334088E+03 6.897916E+02 1.000000E+00 1.878432E+07
3 3 2.811467E+07 5.302327E+03 8.438915E+02 1.000000E+00 2.811467E+07
4 4 1.931709E+08 1.389859E+04 2.212030E+03 1.000000E+00 1.931709E+08
5 5 2.234434E+08 1.494802E+04 2.379052E+03 1.000000E+00 2.234434E+08
6 6 2.328801E+08 1.526041E+04 2.428770E+03 1.000000E+00 2.328801E+08
7 7 6.845924E+08 2.616472E+04 4.164244E+03 1.000000E+00 6.845924E+08
8 8 9.602215E+08 3.098744E+04 4.931804E+03 1.000000E+00 9.602215E+08
9 9 1.365611E+09 3.695418E+04 5.881439E+03 1.000000E+00 1.365611E+09
10 10 1.855386E+09 4.307420E+04 6.855471E+03 1.000000E+00 1.855386E+09

APRIL 8, 1998 MSC.Nastran 4/6/98 PAGE 8
0 SUBCASE 1
1 NORMAL MODES EXAMPLE APRIL 8, 1998 MSC.Nastran 4/6/98 PAGE 9
USING LANCZOS
0 SUBCASE 1
EIGENVALUE = 7.056994E+05
CYCLES = 1.336996E+02 R E A L E I G E N V E C T O R N O . 1

POINT ID. TYPE T1 T2 T3 R1 R2 R3
1 G .0 .0 .0 .0 .0 .0
2 G -3.515604E-14 -2.083550E-14 -9.748799E-01 -1.090191E+00 4.011783E+00 .0
3 G -5.976152E-14 -6.029845E-14 -4.170314E+00 -1.593278E+00 8.679599E+00 .0

4 G -6.872012E-14 -1.060940E-13 -9.445755E+00 -1.531076E+00 1.230244E+01 .0
5 G -6.341396E-14 -1.494814E-13 -1.636251E+01 -1.362968E+00 1.522316E+01 .0
.
.
.
54 G -2.994638E-14 -1.751692E-13 -6.384793E+01 2.199729E-01 2.066229E+01 .0
55 G -2.713030E-14 -1.739597E-13 -7.419337E+01 1.691507E-01 2.070036E+01 .0
.
.
.
.
.
.
EIGENVALUE = 1.878432E+07
CYCLES = 6.897916E+02 R E A L E I G E N V E C T O R N O . 2

POINT ID. TYPE T1 T2 T3 R1 R2 R3

1 G .0 .0 .0 .0 .0 .0
2 G 2.689105E-14 2.049585E-14 -6.305859E+00 5.289161E+00 2.135918E+01 .0
3 G 4.441315E-14 5.120965E-14 -1.933128E+01 1.932312E+01 2.758354E+01 0
4 G 5.540902E-14 8.667080E-14 -3.398484E+01 3.400354E+01 2.819633E+01 .0
5 G 6.112096E-14 1.273261E-13 -4.831832E+01 4.840071E+01 2.654789E+01 .0
6 G 6.240562E-14 1.721406E-13 -6.144451E+01 6.153912E+01 2.363943E+01 .0
7 G 5.830800E-14 2.121037E-13 -7.283444E+01 7.294494E+01 1.992698E+01 .0
.
.
.
53 G -8.359115E-14 2.867785E-13 8.914771E+01 8.923321E+01 -1.121007E+01 .0
54 G -8.422337E-14 3.180499E-13 9.390352E+01 9.367941E+01 -7.136850E+00 .0
55 G -8.453551E-14 3.525693E-13 9.704356E+01 9.580685E+01 -5.385391E+00 .0