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

Расчет центробежного вентилятора. Методика


MotorManiac

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

Добрый день.

Возникли вопросы по расчету центробежного вентилятора.

1. Первый расчет выполнял в стационарной постановке. Использовал функцию frame motion, решатель coupled с функцией  pseudo transient.

В принципе расчет получился, ошибка по сравнение с экспериментом составила 12%

1.jpg.be32ffaa4351e4ee9c68270e6048b8c3.jpg2.jpg.5f18e390b48c1adc226d7dbc393de527.jpg5.jpg.4c416a105162c471ec6d17d2641fa434.jpg

 

2. Хочу произвести расчет в нестационарной постановке с теми же самыми условиями. Использую функцию mesh motion, решатель simple.

В ходе решения получаются такие графики

3.jpg.e7f6c740182ade02121d9f63d660a57c.jpg4.jpg.cda3b022cf90a6014947d00e2f4c9f7e.jpg

 

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

 

 проект     1.rar

 

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


13 minutes ago, MotorManiac said:

В ходе решения получаются такие графики

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

Это к вопросу по поводу пилы на transient графиках.

 

 

 

Но круто, что наконец за нормальный CFD взялся)

28 minutes ago, MotorManiac said:

 проект     1.rar

что-то ни .cas, ни меш файлов нету.

Сам попробуй открыть свой же проект в новой папке...

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

что-то ни .cas, ни меш файлов нету.

Сам попробуй открыть свой же проект в новой папке...

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

 

 

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

вроде запускается

Ну раз так, то ок. Жди того, у кого вроде запустится.

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

выходной патрубок бы удлинить не мешало. прицепи исходную геометрию, кажется, она была в формате x_t.

 

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

прицепи исходную геометрию, кажется, она была в формате x_t.

 

N 315 7-22.x_t

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

посмотрите эту модель, может быть, чем-то поможет. разбил ее на две части, при помощи 7zip можно собрать в одно целое.
дождался массовый расход на выходе в стационарном режиме, получилось, в нестационарном долго решать, нет достаточных мощностей :)

в SIMPLE решает не хуже и быстрее, чем в COUPLED.

 

Geom_steady.zip.001

one more

 

Geom_steady.zip.002

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

посмотрите эту модель

Спасибо за пример.

Может быть подскажете еще по расчету. Как мне посмотреть значения ( мин, макс, среднее.) функции на определенном интервале  ? Наподобие как во флоусим.

6.thumb.jpg.62b5165951172ddec831d68a940d6468.jpg

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

извините, досадная ошибка в координатах MRF зоны, я забыл перевести метры в миллиметры. теперь расхождение едва ли 1% по массовому расходу. файлы прицепляю еще раз.
 

Geom_steady_2.zip.001


еще один

 

Geom_steady_2.zip.002

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

 

 

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

 

Geom_sliding_mesh_2.zip.001


и еще один

 

 

Geom_sliding_mesh_2.zip.002

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

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

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

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

хорошо :)

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

удачи.

 

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

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

Да я обратил на это внимание.  Одинаковые по площади поверхности (интерфейсы) получаются разбивкой поверхностей доменов на одинаковые контактные поверхности ?  А в стационарной постановке на уровне генерации сетки флюент сам автоматически создает интерфейсы, видимо их можно как то отключить ?

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

Хочется увидеть колебания давления

а про модель турбулентности вопросов не возникает?

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

а про модель турбулентности вопросов не возникает?

Я использовал к-е 

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

разве не k-w SST?

при стационарном расчете. в переходном попробовал еще на к-е сделать расчет.

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • Orchestra2603
      Да, в рот мне ноги... Ну, считает же Ансис, считает! Блин.. Как до вас донести простой мэсеж. Мне често, пофигу, делает он там себе факторизацию или еще что-то. Но то, что он нормально считает Ланцошем для случая с вырожденной матрицей жесткости - это факт! А раз он считает, то либо (А) вы хрень несете про факторизацию, либо (Б) он такую факторизацию делат сам без проблем. Третьего тут быть не может, и выходит, что в любом случае ваши тезисы оказываются несостоятельными.   Какая вам нужна факторизация? давайте! заказывайте! На какие вам надо матрицы разложить? с какими свойствами? давайте ваш заказ!   Вы на мой пердак не зарьтесь! Смотрите, чтоб ваш функционировал исправно! А то того и гляди, фекализация случится.  
    • ДОБРЯК
      Вам похоже ничего не нужно, а для алгоритму Ланцоша нужно сделать численную факторизацию. Вам только нужно побольше букв и слов написать. :=) Всё подгорел пердак? :=)
    • Orchestra2603
      При чем здесь это? Речь не идет про эффективные или неэффективнеы алгоритмы. Вы утвержаете, что это просто невозможно.   При чем здесь первые или не первые, все или не все... Какое это имеет значение? У Ансис есть алгоритм Ланцоша, он находит столько, сколько запросите. В Маткаде - он да, находит все. Наверное, можно и не все, но я хз, как это делать. Я не понимаю, какое это имеет отношение к разговору.   Неа Ткните пальцем и объясните, в чем ошибка в конкретном рассуждении   Как я должен решать методом Гаусса задачу на собственные значения? Мне же не решения СЛАУ нужны! Их то, ясен пень, бесконечное множество для системы с вырожденной матрицей. Мне собственные значения нужны! Мне нужно базис найти, чтобы матрица (или матрицы в обобщенном случае) оператора была диагональной.   Сначала скажаите мне, зачем ее дать! В чем идея заключается? Вы все слюной брызжете, но ничего по делу еще не сказали. Какой-то нормальный контраргумент от вас можно толковый услышать?   Как еще изволите вас обслужить? Вы - полегче, уважаемый!
    • ДОБРЯК
      Вам нужно десять раз написать одно и тоже, чтобы вы поняли. Почитайте документацию к программе там же всё написано. Мне не сложно копировать...:=) Те алгоритмы которые вы используете находят все собственные числа и вектора. Для разреженных матриц это приводит к их заполняемости. Смекаете о чем я говорю. Или нужно еще разжевать? Во всех.   Решайте методом Гаусса. Матрица 2х2 в чем проблема сделать численную факторизацию любым методом? :=) Больше дела, меньше слов и букв...
    • mannul
      Скорость резания от 8 до 12 м/мин, подача 0,04 мм/зуб, глубина резания - 1,5 мм.
    • Orchestra2603
      " Стою на асфальте я, в лыжи обутый: То ли лыжи не едут, то ли я ..." (с)   1. Вы утверждаете, что для расчетов на СЗ и СВ в случае вырожденной матрицы жесткости нужно провести ее факторизацию. Без нее - никак! Непонятно, какую и зачем, но надо. Ок, допустим.  2. @Jesse,  я, @Fedor, да и вы сами показываете, что таки СЗ определяются для вырожденной матрицы. Значит, либо ваше утверждение про факторизацию неверно, либо она таки делается тихонько себе за кадром. 3. Еще одно наблюдение... Матлаб говорит, что он делает QZ разложение для матриц, у которых нет положительной определенности, т.е. для которых не работает разложение Холецкого. При этом опять же с собственные числа с одной вырожденной матрицей считаются нормально. Т.е. даже если вдруг и необходима такая факторизация (хотя я сомневаюсь, что вот прям обязательно она нужна), то делается она без особых трудностей.   В каком месте я неверно рассуждаю?   Так и не обязательно их записывать для расчета собственных колебаний. Ну... Можно конечно, но тогде некоторые частоты и формы, которые такой симметрией могут и не обладать, вы потеряете в расчете, поскольку введенные условия симметрии их исключат просто из модели.    
    • gudstartup
    • gudstartup
      у меня рободриллы и чпу и тем не менее picture нет наверно она к роботам относится
    • boomeeeer
      Если компания отечественная, то тоже нет смысла. Роботы, ЧПУ, лазеры и рободрилы
    • gudstartup
      это вряд ли если ось Y угловая а какая у вас редукция у вращающейся оправки? по 2мм нержавейки снимать это нелегкая задача для вашего маленького станочка
×
×
  • Создать...