Динамика - падение груза на консоль, явное решение против теории

[imt 59]

38 минут назад, Victor_M сказал:

налипание груза ИМЕЕТ значение

скорость ВАЖНА

41 минуту назад, Victor_M сказал:

демпфирование "гасит" все же не все частоты

Конечно важно. Я про другое.

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

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

[AlexKaz 1 822]

42 минуты назад, NNK сказал:

можно предположить, что максимальная частота на графиках составляет 52631 Гц.

Ну да (если смотреть на график), в соответствие с теоремой Котельникова-Найквиста-Шеннона, в Вашем тексте - Sampling Theorem. Но это совсем не означает, что частота пилы на графике 52.6 кГц. Не пробовали прогнать через FFT? Если скините данные u(t) - могу построить АЧХ. Ну или в абакусе что-нибудь есть чтоб построить АЧХ...

 

42 ... прям как в известном фильме...

Изменено 18 января 2018 пользователем AlexKaz

[NNK 9]

Я попробовал немного упростить задачу - ничего не падает, нагрузкой является мгновенная сила 1000 Н. Коэффициент динамичности в таком случае равен 2. Целевое перемещение 6 мм, напряжение 96 МПа, сила реакции 2000 Н.

Вот что получается:

Перемещение совпадает.

Напряжение достигает 110 МПа.

Реакция достигает 4500 Н.

Применение фильтра по первой собственной частоте (42 Гц) приводит к заниженным значениям.

Что-то я делаю не так...

Изменено 18 января 2018 пользователем NNK

[NNK 9]

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

Ну да (если смотреть на график), в соответствие с теоремой Котельникова-Найквиста-Шеннона, в Вашем тексте - Sampling Theorem. Но это совсем не означает, что частота пилы на графике 52.6 кГц. Не пробовали прогнать через FFT? Если скините данные u(t) - могу построить АЧХ. Ну или в абакусе что-нибудь есть чтоб построить АЧХ...

 

42 ... прям как в известном фильме...

42... да уж, может получиться весьма символично )

К сожалению, Абакус сам по себе не умеет делать FFT и строить АЧХ (по крайней мере, насколько мне известно).

Экселька не пролезает в разрешенный размер файла, попробую ее покромсать на части.

[Борман 2 375]

Предлагаю просто интегрировать с шагом, который соответствует частоте 9Гц... сделайте 20 временных точек на период и хватит.

[NNK 9]

Преобразовал в текстовый файл, надеюсь это не вызовет неудобств. Это запись реакции в заделке от времени. Спасибо за помощь!

Книга1.7z

[soklakov 1 475]

Только что, Борман сказал:

Предлагаю просто интегрировать с шагом, который соответствует частоте 9Гц... сделайте 20 временных точек на период и хватит.

 

23 часа назад, NNK сказал:

Abaqus/Explicit (решатель не принципиален, главное - явный с записью результатов на каждом шаге).

явные решалки сами выбирают шаг. и 20 точек на период не прокатит. тут больше скорость распространения упругой волны решает и сетка.

курант, фридрихс, леви

ленин, партия, комсомол

[Борман 2 375]

Только что, soklakov сказал:

явные решалки сами выбирают шаг. и 20 точек на период не прокатит. тут больше скорость распространения упругой волны решает и сетка.

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

[soklakov 1 475]

6 минут назад, Борман сказал:

но если цель задавить высокую частоту и приблизить численное решение к аналитическому...

то взять неявную решалку?

кстати, первая частота чуть меньше 41 Гц получилась... (40,689 Гц)  - странно, что не 42 или норм?

 

 

54 минуты назад, NNK сказал:

Целевое перемещение 6 мм, напряжение 96 МПа, сила реакции 2000 Н.

По перемещениям и напряжениям так-то близко... а вот реакция... не, ну порядок-то совпадает 

 

Изменено 18 января 2018 пользователем soklakov

[AlexKaz 1 822]

100 000 отсчётов не требуется, достаточно и 10 000. Полная АЧХ на первой картинке внизу. На второй картинке АЧХ в области нижних частот.

Изменено 18 января 2018 пользователем AlexKaz

[imt 59]

15 минут назад, soklakov сказал:

... взять неявную решалку?

А Analysis Settings, можно посмотреть?

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

АЧХ внизу

И какой вывод?...

[soklakov 1 475]

1 минуту назад, imt сказал:

А Analysis Settings, можно посмотреть?

А что интересует?

Геометрическая нелинейность - да.

Демпфирование - да (грустных 2% просто для красоты).

Шаг по времени на графиках виден... можно любой так-то, но лучше 

40 минут назад, Борман сказал:

сделайте 20 временных точек на период и хватит

или более.

[AlexKaz 1 822]

Система многомассовая, частот много. Сопромат рулит, но не в Вашем случае, и софт ловит всё подряд - в данном случае нет ошибки.

Подозреваю, у Вас чересчур много степеней свободы для столь простой задачи - уменьшите размерность и жить станет полегче. А шаг выбрать побольше =)

[imt 59]

6 минут назад, soklakov сказал:

Демпфирование - да (грустных 2% просто для красоты).

А если сделать демпфирование 0 (чтобы было веселей:)

Изменено 18 января 2018 пользователем imt

[soklakov 1 475]

Только что, imt сказал:

А если сделать 0 (чтобы было веселей:)

Да разве ж это весело? (там численного копейку осталось, если что)

Ну не бывает демпфирование меньше) 1-2% есть всегда.

[imt 59]

За что люблю ANSYS: хочешь - красиво, хочешь - весело:)

2 минуты назад, soklakov сказал:

Ну не бывает демпфирование меньше) 1-2% есть всегда.

Я понимаю.

Осталось @AlexKaz попросить сделать AЧХ...

Пожалуйста.

[soklakov 1 475]

Только что, imt сказал:

За что люблю ANSYS: хочешь - красиво, хочешь - весело:)

да Ansys-то тут при чем))) демпфирование везде задать можно. и везде можно не задавать)

1 минуту назад, imt сказал:

Осталось @AlexKaz попросить сделать AЧХ...

так он сделал вроде уже...

@imt меня больше смущает, что по частоте разница 1.5 Гц... я уже и в солидах вроде прикинул, а 42 не получается

[AlexKaz 1 822]

Вас смущает, что на графиках нет всплеска на 9 Гц? Вот он, если взять больше отсчётов

 

И на АЧХ нет 42.

[imt 59]

7 минут назад, soklakov сказал:

демпфирование везде задать можно

А я так и не понял у @NNK есть демпфирование или нет?

Уважаемый @AlexKaz ! Если у @soklakov взять реализацию, где ANSYS дребезжит, что АЧХ покажет?

Прошу прощения за настойчивость...

11 минуту назад, soklakov сказал:

42 не получается

Как я понял, про 42 мы знаем со слов @NNK ...

[soklakov 1 475]

2 минуты назад, imt сказал:

А я так и не понял у @NNK есть демпфирование или нет?

У него промоделирована целая секунда. И  амплитуда колебаний явно не собирается падать... Какой можно сделать вывод?

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

Уважаемый @AlexKaz ! Если у @soklakov взять реализацию, где ANSYS дребезжит, что АЧХ покажет?

Прошу прощения за настойчивость...

ну, для начала мне придется выложить эти отсчеты, чтобы @AlexKaz смог сделать с ними FFT) 

 

только их нет смысла сравнивать