Определение частоты колебаний с нелинейностью в модели

[Jesse 977]

Всем привет. 
При расчёте собственных частот (СЧ) мы обычно пренебрегаем нелинейностями, которые в той или иной степени присутствуют в любой конструкции.
А если заказчик требует оценить СЧ конструкции, где нелинейность - часть её работы?
К примеру, рассмотрим колебания бака на опоре, который стоит на земле.

Или даже ещё проще пример: консоль, которая при колебаниях влево-вправо будет напираться на жёсткие стенки (скоммуниздил картинку из смежной темы).

Как бы вы посчитали частоты такой конструкции? Ну т.е. понятно, что обычный линейный модальник тут не канает. Linear Pertubation тоже, т.к. у нас при колебаниях зона контакта меняется.
Получается остаётся только нелинейный динамический расчет? Задаём перемещение, "отпускаем" (рассматриваем свободные колебания), смотрим как меняется амплитуда. И в конце FFT? Кто-нибудь так делал?
Насколько физически корректно говорить о частоте колебаний в данном случае?

[Fedor 1 625]

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

[AlexKaz 1 848]

2 часа назад, Jesse сказал:

Насколько физически корректно говорить о частоте колебаний в данном случае?

Всё корректно, если строить не плоский спектр по всей выборке отклика, а водопад, нарезая отклик на куски, чтобы смотреть а будет ли изменяться спектр в моменте при непосредственно контакте.

Шаг по времени по теореме Котельникова-Найквиста-Шеннона надо брать <= 0.5/f, где f - максимальная интересуемая частота в Гц.

Демпфирование можно задать стандартно 2-3-4%.

 

Если transient и fft не покажут принципиальной разницы между чисто модальным - на transient можно забить.

Transient  + FFT есть большой смысл применять если есть неупругая деформация, вот тогда частоты могут сильно измениться

[Jesse 977]

56 минут назад, Fedor сказал:

В нелинейном случае частота обычно зависит от амплитуды

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

Мой вопрос немного другой: могу ли я использовать результаты МКЭ по амплитуде из нелинейного динамического анализа, а затем просто по набору данных A(t) сделать дискретное преобразование Фурье? ну или спектрограмму какую-то построить...

Вот мой коллега гидродинамик по такому графику осцилляций давления в насосе, посчитанному в Star-CCM, взял FFT и построил спектр частоты изменения давления...

 

Изменено 1 час назад пользователем Jesse

[AlexKaz 1 848]

И если есть параметрические колебания, как выше написал Fedor. Теоретически, "метроном" (хотя, это не он) на рисунке может впасть в параметрические колебания, найти этот режим через transient проще.

4 минуты назад, Jesse сказал:

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

Если частота растёт, то период колебаеий уменьшается, колебания становятся чаще. А на графике всё наоборот.

[Jesse 977]

23 минуты назад, AlexKaz сказал:

Если transient и fft не покажут принципиальной разницы между чисто модальным - на transient можно забить.

Transient  + FFT есть большой смысл применять если есть неупругая деформация, вот тогда частоты могут сильно измениться

 А как вы моделировали в Трензиенте нагрузку? Вынужденные или свободные колебания? Или вообще ударно делали (ступеньку Хевисайда/дельта функцию)?

[pugach 9]

Цитата

В данном случае со свободными колебаниями консоли со временем площадка контакта будет уменьшаться, и частота должна расти. 

Если тело при колебаниях задевает другое тело (берем пример с картинки), значит оно "хочет" колебаться с большей амплитудой и меньшей частотой. И когда контакт пропадет оно должно колебаться со своей собственной 1 частотой, или я не прав?

1 час назад, AlexKaz сказал:

Демпфирование можно задать стандартно 2-3-4%.

Не знаю, может, отдельную тему про это задать. Я намучился с заданием "стандартного" демпфирования в workbench.

Тестил на колебаниях заделанного стержня. Пробовал считать альфа и бетта, но в результате получал значительно большее затухание - логарифмический декремент колебаний ~0,35 вместо 0,12.

Может это быть из-за того, что при моих условиях в стержне были небольшие (до 1 МПа) напряжения?

Ведь НП-031 предлагает следующую таблицу, получается стандартные 2% относятся к стандартным условиям по нагрузкам на конструкцию?

[Jesse 977]

8 минут назад, pugach сказал:

Если тело при колебаниях задевает другое тело (берем пример с картинки), значит оно "хочет" колебаться с большей амплитудой и меньшей частотой. И когда контакт пропадет оно должно колебаться со своей собственной 1 частотой, или я не прав?

по идее наоборот: когда тело контачит с жёстким телом, то его частота завышается.. Если свободно болтается, то частота ниже должна быть..)) Ну т.е. как обычно рассуждаем в МКЭ: больше суммарная жёсткость модели - больше частота 
 

3 часа назад, Jesse сказал:

Насколько физически корректно говорить о частоте колебаний в данном случае?

вот тут всё же кажется, что если мы посчитаем transient с нелинейностью а потом возьмём FFT с графика амплитуды, то это уже халтура будет)) Ну т.е. когда мы считаем модальник обычный, то это линейный идеализированный случай, когда мы пренебрегли всеми возможными нелинейностями. Иначе говоря, малая нелинейность которая присутствует всегда и везде сделает колебания лишь слегка ангармоничными (очень слабые отклонения от чистой синусоиды при измерении колебаний в точке)..

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

[pugach 9]

10 минут назад, Jesse сказал:

по идее наоборот: когда тело контачит с жёстким телом, то его частота завышается.. Если свободно болтается, то частота ниже должна быть..)) Ну т.е. как обычно рассуждаем в МКЭ: больше суммарная жёсткость модели - больше частота

Я это и имел ввиду, вроде :)