Перейти к публикации

Расчет вибрации при одновременном действии усилий на нескольких частотах


Рекомендованные сообщения

Добрый день, форумчане

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

Ссылка на сообщение
Поделиться на других сайтах


UnPinned posts

Если задача линейная, то сначала считаете на одной частоте и находите виброускорение на ней, а потом на другой.

Ссылка на сообщение
Поделиться на других сайтах
2 часа назад, Борман сказал:

Если задача линейная, то сначала считаете на одной частоте и находите виброускорение на ней, а потом на другой.

И потом можно сложить с учетом фазы?

Ссылка на сообщение
Поделиться на других сайтах
1 минуту назад, Nikk24 сказал:

И потом можно сложить с учетом фазы?

Можно думаю. 

Но смысла в этом будет мало если разные частоты.

Ссылка на сообщение
Поделиться на других сайтах

Не знаю какую программу вы используете.

Посмотрите прилепленные файлы - пластина нагруженная в 2-х точках сосредоточенными силами

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

Это -Frequency Response Analysis - я использtest_1.jpgtest_2.jpgtest_3.jpgtest_4.jpgtest_5.jpgtest_6.jpgtest_7.jpgtest_8.jpgNormal Moded.txtFRAD.bdfовал Direct но можно и Modal

Обратите внимание , что силы действуют одновременно

Ссылка на сообщение
Поделиться на других сайтах
3 часа назад, Nikk24 сказал:

И потом можно сложить с учетом фазы?

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


Что касается фаз, но в контексте прочности на них обычно забивают насколько знаю.. То есть берут обычно только АЧХ (действительную часть) комплексного отклика, и работают с амплитудным спектром или спектром мощности. Разности фаз нужны чтоб по отклику восстановить исходный сигнал. Короче это больше для тех кто с изображениями работает и смежными фичами.. 
Буду рад кто поправит насчёт фаз, самому интересно..)

Ссылка на сообщение
Поделиться на других сайтах
10 hours ago, Nikk24 said:

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

В линейных задачах действует принцип суперпозиции в том числе и в частотной области: считайте на набор сил F_i на частоте w_i, потом второй раз на набр сил F_j на частоте w_j. Результат можно смело складывать.

4 hours ago, Nikk24 said:

И потом можно сложить с учетом фазы?

Как правило, фазовые соотношения между откликами от разных сил (например, от работающего двигателя и от каких-то внещних динамических усилий) даже на одной частота неизвестны - это штука случайная. Поэтому, в принципе невозможно точно сказать, как складывать амплитуды гармоник (скажем, A и B ) от двух раченых случаев даже на одной и той же частоте. Может быть и A - B (в противофазе), и A+B, если совершенно синфазный случай. Можно взять сдвиг по фазе как некое среднее на интервале от 0 до П, и суммирующая амплитуда будет sqrt(A^2 + B^2). Если складывть гармоники с разными частотами, то разница фаз будет еще и во времени меняться. Если в результате нужна максимальное суммарное значение, то я бы брал в расчет просто сумму всех амплитуд (это гарантированный в функциональном анализе максимум), но это была бы очень консервативная оценка.

11 hours ago, Nikk24 said:

...в таком случае не применим и нужно использовать трансиент?

Если вам достоверно известны фазовые сдвиги по силам на разных частотах в разных точках, и нужно только снять c time history максимальный суммарнй отклик, то может иметь смысл. Я бы сильно задумался тогда, а что если эти сдвиги по фазе по силам окажутся чуть-чуть не такими, насколько изменится максимальный отклик.

Ссылка на сообщение
Поделиться на других сайтах

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

Транзиентный анализ хорош если у вас НЕПЕРИОДИЧЕСКОЕ нагружение - как правило удар

 

Ссылка на сообщение
Поделиться на других сайтах
23 минуты назад, dbarlam сказал:

Транзиентный анализ хорош если у вас НЕПЕРИОДИЧЕСКОЕ нагружение - как правило удар

а делали когда - нибудь в МКЭ трензиент для НЕПЕРИОДИЧЕСКОГО нагружения в частотной области? Ну типа спектрограммы и всё такое...

Ссылка на сообщение
Поделиться на других сайтах

Вы наверное имеете ввиду SRS - Shock Response Spectrum Преобразовние/представление непериодического сигнала 

в спектральной области

Если вы об этом - тогда ДА -делал

Имелся импульс как функция времени и я его переводил в частотную область и использовал

полученный график для расчета с использованием SRS - получались результаты 

И это вместо того чтобы гонять транзиентную задачу. Много проще и быстрее. 

Есть свои подводные камни , но в принципе подход известный и удобный а главное быстрый

Ссылка на сообщение
Поделиться на других сайтах
37 minutes ago, Nikk24 said:

Проблема, вроде, решена.

Так и как решена в итоге? Расскажите уж.

Ссылка на сообщение
Поделиться на других сайтах
35 минут назад, dbarlam сказал:

- получались результаты 

И это вместо того чтобы гонять транзиентную задачу. Много проще и быстрее. 

а весёлые картинки имеются?)

 

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

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

можно заключить, что фазовые соотношения интересны только когда у нас несколько сил? То есть когда у нас одна сила, то на фазы можно забить (что обычно и делается)..

Ссылка на сообщение
Поделиться на других сайтах

Конкретно сейчас под рукой ничего нет.

Я завтра подготовлю что-нибудь попроще - и поставлю на форуме. 

Ссылка на сообщение
Поделиться на других сайтах
2 часа назад, dbarlam сказал:

Вы наверное имеете ввиду SRS - Shock Response Spectrum Преобразовние/представление непериодического сигнала 

чё т тупанул сразу не вкурил..)
Не, спектралка - это методика, иначе говоря костыль для расчёта сейсмики.
Я про другое.. конкретно,  перевести входной импульс в частотную область, получив тем самым амплитудный входной спектр, меняющийся во времени. И для каждого момента времена (и, соотв-но, для каждого "мгновенного" спектра) вычисляем частотный отклик по обычному Harmonic Response. 

Spectrogram-19thC.pngSpectrogram.png

 

Это чем то похоже на анализ Linear Pertubation, когда для каждого шага в нелинейном решателе определяем собственные частоты или формы потери устойчивости. Но по ходу то, что я выше расписал крайне невыгодно и бессмысленно использовать..)

Изменено пользователем Jesse
Ссылка на сообщение
Поделиться на других сайтах
2 часа назад, Orchestra2603 сказал:

Так и как решена в итоге? Расскажите уж.

Неточно выразился - не задача решена, а определен способ ее решения. А сама она еще будет считаться. 

Ссылка на сообщение
Поделиться на других сайтах

коллега Jessy - Я вас понял . Мы это делали на другом инструментарии. Если мне не изменяет память - VIEWLAB или что-то похожее по

названию. Лет 10 назад. 

Ссылка на сообщение
Поделиться на других сайтах

Это называется LABVIEW - мы использовали его для тестирования вибраций конструкции

Там очень широкий спектр возможностей 

И есть связь с FEA

Ссылка на сообщение
Поделиться на других сайтах
1 час назад, dbarlam сказал:

Это называется LABVIEW

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

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

Там очень широкий спектр возможностей 

И есть связь с FEA

Многоканальная замерялка данных с портов + обрезанный аналог симулинка. Если и есть связь с FEA, то явно лишняя. Загнать временнЫе данные из CAE в Scilab/Matlab и т.п. и получить водопад куда проще. Но это мало кому надо - тот же ANSYS БПФ и водопад умеет сам искаропки.

Изменено пользователем AlexKaz
Ссылка на сообщение
Поделиться на других сайтах
21 hours ago, Jesse said:

можно заключить, что фазовые соотношения интересны только когда у нас несколько сил? То есть когда у нас одна сила, то на фазы можно забить (что обычно и делается)..

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

Ссылка на сообщение
Поделиться на других сайтах

Присоединяйтесь к обсуждению

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

Гость
Ответить в тему...

×   Вставлено в виде отформатированного текста.   Вставить в виде обычного текста

  Разрешено не более 75 эмодзи.

×   Ваша ссылка была автоматически встроена.   Отобразить как ссылку

×   Ваш предыдущий контент был восстановлен.   Очистить редактор

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

  • Сейчас на странице   0 пользователей

    Нет пользователей, просматривающих эту страницу.




  • Сообщения

    • SAPRonOff
      https://dzen.ru/b/ZjnRssPLbUcXVOUZ
    • davidovka
      Creo 9, VS 2019 (v142). Если использовать стандарт языка С++14, все собирается работает. При переключении на C++17 - ошибки (см. скрин - тестовое приложение из папки Creo) Хотел в проекте использовать <filesystem>
    • maxx2000
    • ДОБРЯК
      Нулевых частот может быть любое количество. Вот эта матрица не будет вырождена.    
    • Orchestra2603
      Речь идет о собственных векторах, связанных с нулевым СЗ. Вы же понимаете.
    • ДОБРЯК
      Для сложных систем их может быть и 100 и 200 и 1000...
    • maxx2000
      @Artas28так надо рецепт в массы нести, а то Ванёк ничему не научится
    • Orchestra2603
      Выкидываем не в смыле, что мы ее вырезаем.. Простите мне эту вольность)) Просто там две одинаковые строки. Достаточно рассмотреть только одну. Нет смысла хранить вторую строку в матрице, поскольку она ничего не привносит. Вы можете совершенно формально взять полученный i-й собственный вектор, умножить на матрицу K-lambda_i*M и получиь нулевой вектор.   На самом деле.. Мне прям пришлось прям подумать над аналогией с включением закрепления в статике. И это было полезно и интересно. Так что вам спасибо за "подгон", что называется )) Но все-таки разница есть ...   Вот у вас в статике:   И здесь действительно получается, что если мы зафиксируем, например, положив a=0, то это будет эквивалентно "выкидыванию" 1-го стобца и 2-й строки, и решению уже такой модифицированной невырожденной системы *в данном случае вообще тривиальной). При этом f2 просто станет узловой реакцией и определится как -f1 автоматически: К слову, можно положить любой a, не только нулевой. Результат будет такой же для сил, но прогибы, ясное дело, поменяются.   Но это отличается от нашей задачи на собственные значения. У нас нет внешних сил {f}.  А знaчит и упругие силы [K]{x} в системе должны быть равны нулю. И для нулевой собственной частоты мы по сути ищем такой набор перемещений, при которых в системе отсутствуют упругие силы.  И таких наборов перемещений бесконечно много. Вообще говоря, для трехмерных более сложных систем оказывается, что существует до 6-ти базисных векторов (понятно, что связанных с движениями как твердого целого), и любая линейная комбинация из этих векторов будет решением такой задачи. Но нам не нужно само решение! Нам нужен как раз сам базис пространства решений! Причем, сами вектора в этом базисе нас интересуют с точностью до множителя, ь.е. таких базисов можно напридумывать бесконечно много. Но нам нужен толкьо какой-то один базис, и всего-то. В этой простейшей задачке - это просто один вектор  {1, 1}. Можно взять {2,2} или {e,e} - не имеет значения, нас любой из них строит. И нам этого достаточно.    Обратите внимание, что этот вектор нельзя получить из задачи статики, где {f} ≠ {0}, положив какое-то значение для a, т.е. внедрив какое-то "закрепление" в вашей интерпретации.  Пожтому я убежден, что это другая задача, т дело тут не в закреплениях.   Не понимаю, про какую вы матрицу
    • Viktor2004
      Туда вроде наоборот заносится значение какое-то, но не считывается Смотрю чужие программы, там вроде только 8549 используют. При чем у них в программах циклы поболее моих, которыми я круги чертил O4198 ; #500=#8549 ; IF[ #8549 EQ 0 ] THEN ; #10409=0 ; ENDIF ; IF[ #8549 EQ 134 && #10409 EQ 0 && #10422 NE 0 ] THEN ; #10409=1 ; #10422=FIX[#10422 / 10] ; #10427=1 ; ENDIF ; IF[ #8549 EQ 48 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 0 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 49 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 1 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 50 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 2 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 51 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 3 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 52 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 4 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 53 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 5 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 54 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 6 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 55 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 7 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 56 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 8 ; #10427=1 ; ENDIF ; IF[ #8549 EQ 57 && #10409 EQ 0 ] THEN ; #10409=1 ; #10422=#10422 * 10 + 9 ; #10427=1 ; ENDIF ; IF[ #10422 GT 9999 ] THEN ; #10422=FIX[#10422 / 10] ; ENDIF ; IF[ #500 EQ 249 && #8500 EQ 4170 ] THEN ; #8500=4601 ; M99 ; ENDIF ; IF[ #500 EQ 240 ] THEN ; IF[ #8500 EQ 4174 || #8500 EQ 4173 || #8500 EQ 4177 ] THEN ; #8500=4170 ; M99 ; ENDIF ; IF[ #8500 EQ 4178 ] THEN ; #8500=4170 ; M99 ; ENDIF ; ENDIF ; IF[ #500 EQ 248 && #8500 EQ 4170 && D34.0 NE 0 ] THEN ; #8500=4174 ; M99 ; ENDIF ; IF[ #500 EQ 247 && #8500 EQ 4170 && D30.4 EQ 1 && D29.6 EQ 0 ] THEN ; #8500=4175 ; M99 ; ENDIF ; IF[ #500 EQ 246 && #8500 EQ 4170 && D35.0 EQ 1 && D35.4 EQ 1 ] THEN ; #8500=4178 ; M99 ; ENDIF ; IF[ #8549 EQ 152 && #10409 EQ 0 && #10422 EQ 5656 ] THEN ; #8500=4177 ; M99 ; ENDIF ; M99 ;  
    • nicomed
      Как победить проблему с блокированием сегментов эскиза не знаю. Поэтому пока полумеры... Ввел активацию чертежного листа после разрыва линии. Для разблокирования сегментов после закрытия программы двойным кликом активировать чертежный вид, где надо двигать сегменты эскиза. Добавил настройки. Можно настроить величину генерируемого зазора в  файле options.xml. По толщине и стилю линии, если кто подскажет как прочитать стиль и толщину выбранной кромки ( силуэтной кромки), то тогда можно будет "напосмотреть".   Пару раз натыкался на такой глюк - остается только один сегмент эскиза, грешу на рассинхронизацию работы солида и приложения. EdgeTrim_v3.3.zip
×
×
  • Создать...