FFT график при исследовании пульсации давления

[rasta89 113]

Всем привет.

Стоит задача исследовать пульсацию давления, определить ее частоту.

Задачу решал при помощи модели DES, без уравнения энергии для ускорения. Замерял давления: 1. AreaAve(TotalPressure) полная поверхность выходного отверстия из детали 2. там же на поверхности точка, замер в точке через surface-vertex-ave -total pressure.

Задача решалась transient анализе. Для теста прогнал на трех вариантах времени шага: 300 шагов с временем шага 0,001с, 300/0,01 и 300/0,1.

При помощи FFT получал график зависимости частота/средний импульс давления. Вот тут и проблема. Не определяется частота, она всегда идет от нуля на понижение.

Приложил по 4 скрина на каждый расчет. Два это картина давления на поверхности и точке и два это полученные из них частоты:

 

[spoiler=300/0,001]

[/spoiler]

[spoiler=300/0,01]

[/spoiler]

 

[spoiler=300/0,1]

[/spoiler]

 

Так выглядит сама модель

 

Подскажите пожалуйста, что пошло не так касательно ускорения на графике FFT?

Что-то спойлер не сработал

[rasta89 113]

В общем проблема была в том, что постпроцессор флюента тупил получается.

Сформировал график в CFD POST там показало эти самые 0,33 Гц и импульс 0,04 МПа

Изменено 27 сентября 2022 пользователем rasta89

[Jesse 958]

14 минут назад, rasta89 сказал:

постпроцессор флюента тупил получается.

так преобр-е Фурье вы прям в Ансисе делали получается? я то думал это какое-то стороннее ПО

 

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

Сформировал график в CFD POST там показало эти самые 0,33 Гц и импульс 0,04 МПа

странно что других мелких пиков нет, а всего 1. Графики по времени не похожи на "квази"синус. Мб использовать STFT (оконку)? 

[rasta89 113]

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

так преобр-е Фурье вы прям в Ансисе делали получается? я то думал это какое-то стороннее ПО

да, встроенным инструментом

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

странно что других мелких пиков нет, а всего 1. Графики по времени не похожи на "квази"синус. Мб использовать STFT (оконку)? 

можно подробнее пожалуйста?

[Jesse 958]

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

можно подробнее пожалуйста?

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

[rasta89 113]

5 часов назад, Jesse сказал:

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

Да вот в постпроцессоре флюента он там один и потом очень мелко нисходящий

А вот в CFD post там уже один пик да.

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

[karachun 2 038]

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

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

Вы же уже давно занимаетесь численными методами. Вы же прекрасно знаете что в расчетных пакетах ничего "по идее" не выбирается. Все нужно задавать самому.

 

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

[rasta89 113]

4 часа назад, karachun сказал:

Вы же уже давно занимаетесь численными методами. Вы же прекрасно знаете что в расчетных пакетах ничего "по идее" не выбирается. Все нужно задавать самому.

 

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

"По идее" - я имел в виду я задал этот отрезок сам при формировании графика, отсек часть графика при котором поток жидкости только становился ( старт задачи).

Нужно поработать над своим изложением мыслей)

[rasta89 113]

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

"По идее" - я имел в виду я задал этот отрезок сам при формировании графика, отсек часть графика при котором поток жидкости только становился ( старт задачи).

Нужно поработать над своим изложением мыслей)

Первый график с плоскости выхода, второй с точки

[rasta89 113]

Решил проверить график в программе origin.

Загнал туда три варианта графиков

1. Сбор информации с точки 1.5-4 сек, клгда течение установилось

2. Сбор информации с поверхности выхода жидкости 1.5-4 сек, когда течение установилось

3. Две волны с поверхности 2.1-3с

И сравнил с посьпроцессором флюента и CFD post. Теперь незнаю чему верить.

Внизу приложил скрины последовательно

1. Выбранный промежуток

2. Результат в постпроцессоре флюента

3. Результат в cfd

4. Результат в origin

Также приложил файлы графиков.

Как понять что из чего правда?

Спасибо.

 

Вариант 1

 

 

Вариант 2

 

 

 

 

 

Вариант 3

 

 

 

Файлы

short_area_fluent_2.1-3slong_area_1.5-4.5slong_vertex_1.5-4.5s

[rasta89 113]

В общем так.

Разобрался.

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

2. Постпроцессор флюета не обманывает, просто отображает график по-своему

3. Постпроцессор CFD показывает картину более понятной и наглядной. Так что если флюент показал непонятную Вам картину идите в CFD POST. Только там есть нюанс, амплитуду показывает параметр Power spectral

Скрины раскладки Фурье в CFD POST и ORIGIN

 

[Arkur 15]

Небольшой дилетантский вопрос. Почему используете окно прямоугольное, а не Хэннинга?

[soklakov 1 475]

28.09.2022 в 08:47, rasta89 сказал:

Вариант 1

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