Jump to content

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


Recommended Posts

rasta89

Всем привет.

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

Задачу решал при помощи модели 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]

p1.PNGtp1.PNGp2.PNGtp2.PNG

[/spoiler]

[spoiler=300/0,01]

p1.PNGtp1.PNGtp2.PNGp2.PNG

[/spoiler]

 

[spoiler=300/0,1]

p1.PNGtp1.PNGtp2.PNGp2.PNG

[/spoiler]

 

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

me.PNG

 

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

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

Link to post
Share on other sites


UnPinned posts
rasta89

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

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

photo_2022-09-27_15-36-27.jpg

Edited by rasta89
Link to post
Share on other sites
14 минут назад, rasta89 сказал:

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

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

 

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

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

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

Link to post
Share on other sites
rasta89
5 минут назад, Jesse сказал:

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

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

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

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

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

Link to post
Share on other sites
6 минут назад, rasta89 сказал:

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

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

Link to post
Share on other sites
rasta89
5 часов назад, Jesse сказал:

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

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

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

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

Link to post
Share on other sites
karachun
1 час назад, rasta89 сказал:

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

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

 

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

Link to post
Share on other sites
rasta89
4 часа назад, karachun сказал:

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

 

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

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

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

Link to post
Share on other sites
rasta89
1 час назад, rasta89 сказал:

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

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

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

tp-area.png

tp-vertex.png

  • Нравится 1
Link to post
Share on other sites
rasta89

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

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

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

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

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

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

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

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

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

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

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

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

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

Спасибо.

 

Вариант 1

long_vertex_fluent_1.-4.5s_scr.PNG

long_vertex_fluent_1.-4.5s.PNG

 

long_vertex_1.5-4.5s_fluent_post.PNG

long_vertex_origin_1.5-4.5s.PNG

 

Вариант 2

 

 

 

long_area_1.5-4.5s_fluent_scr.PNG

long_area_1.5-4.5s_fluent.PNG

 

long_area_1.5-6s_fluent_post.PNG

long_area_1.5-4.5s_origin.PNG

 

Вариант 3

short_area_fluent_2.1-3s_scr.PNG

short_area_fluent_2.1-3s.PNG

 

short_area_fluent_2.1-3s_post.PNG

short_area_origin_2.1-3s.PNG

 

 

Файлы

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

Link to post
Share on other sites
rasta89

В общем так.

Разобрался.

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

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

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

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

SVG_hz.PNG

SVG_hz_origin.PNG

 

Link to post
Share on other sites

Join the conversation

You can post now and register later. If you have an account, sign in now to post with your account.
Note: Your post will require moderator approval before it will be visible.

Guest
Reply to this topic...

×   Pasted as rich text.   Paste as plain text instead

  Only 75 emoji are allowed.

×   Your link has been automatically embedded.   Display as a link instead

×   Your previous content has been restored.   Clear editor

×   You cannot paste images directly. Upload or insert images from URL.

  • Recently Browsing   0 members

    No registered users viewing this page.




  • Сообщения

    • gudstartup
      Проанализируйте программу смены и определите M код который отвечает за подъем зонтика но при применении этого кода в MDI нужно либо указать ему сервис режим для магазина либо M код  который активирует цикл ATC. Обычно он указан в начале подрограммы смены инструмента. Если хочется вручную то можете и так попробовать но это зависит от пневмо-схемы если используются два клапана то нужно сначала один закрыть а другой открыть иначе клапан который подает воздух в одну полость цилиндра не закроется и открытие другого ничего не даст. 
    • gudstartup
      ШВП конечно соответствует своему классу точности и ее изготовитель гарантирует отклонения которые укладываются в допуски но для высокоточных станков с измерительной системой  B C и т.д. этого недостаточно поэтому дополнительно применяется коррекция шага. Функция Dual positon feedback имеет диагностические параметры  Если линейка установлена без значительных перекосов в соответствии с инструкцией по инсталляции тогда ее показания использовать для определения рассогласования данных от датчика мотора и от линейки. Проехали шаг винта посмотрели сколько посчитал датчик мотора а сколько линейка и можно вычислить коррекцию шага но до этого все компенсационные импульсы которые генерирует ЧПУ должны быть отключены а люфты в ШВП быть порядка 1-3 мкм.
    • Viktor2004
      Хотя вряд ли. Как потом координаты совмещать. Я думаю можно как-то по оптической линейке продиагностировать
    • Viktor2004
      А отдельно ШВП пока оно не установлено, кто-нибудь диагностирует?
    • A_1
      Используйте цикл G71  
    • soklakov
    • Igoshechka
      Необходимо провести расчёт для одномерной задачи со скачком. Есть стержень с областями  V1 и V2, между ними граница. На границе узлы имеют двойную нумерацию. Есть два способа, первый - прописать в ANSYS условие контакта на границе, но я не знаю как это сделать. (Буду рад, если кто-нибудь подскажет.) Знаю, что условие контакта имеет вид: ui (на левой границе) = ui (на правой границе) - Ui, где Ui - скачок в перемещениях, который задан по условию задачи.   Т.к. не нашёл способ прописать условие контакта в ANSYS, то был выбран способ извлечения глобальной матрицы жёсткости из ANSYS, её модификации и загрузке обратно.   Я получил из ANSYS Mechanical APDL глобальную матрицу жёсткости с помощью следующих двух команд: *DMAT, MatKD, D, IMPORT, FULL, file.full, STIFF *PRINT, MatKD, Kdense.matrix   После чего отредактировал её и хочу загрузить обратно. Какие команды для этого использовать? Также есть ещё один вопрос: как получить из ANSYS вектор нагрузки F? (Который используется вместе с матрицей жёсткости A*U=F.) Заранее спасибо.
    • Fedor
      Утром в воскресенье хороша чашечка кофе... :) 
    • maniokka
      BSV1, огромное спасибо за наводку! Сделал как адаптивный элемент, можно при вставке в сборку регулировать длину, очень удобно. Но проблема с углами осталась, не совсем понимаю что надо сделать, чтобы они менялись. Экспериментировал с различными видами вычитания, но без толку.  К сожалению, поиск по теме особых результатов не принес. Не подскажите, откуда вы вытащили эту гифку? В принципе, это и есть то что мне надо. Только у меня этот брусок стыкуется с небольшим нахлестом, т.е. на торце еще Г-образный вырез. Но думаю это не важно, главное понять сам принцип, а там разберусь. Возможно там где была эта гифка и есть необходимая мне информация?
    • mrvcf1
      Спасибо, попробую, думал она только для сброса.
×
×
  • Create New...