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

расчет сопла Лаваля


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

Занятная задачка.

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

Сейчас новый решатель все это проглатывает не кашляя.

В чем проблема. Проблема в том, что волна давления доходит до сужения сопла, где площадь сечения критики меньше по площади чем оставшаяся часть сечения. ВОЛна отражается от стенки конфузора и идет обратно. Доходит в виде ударной волны до входа, где задано давление и с этого места начинается "Страх и ненависть в Лас-Вегасе".

Все существующие методики, которые были тут перечислены, годятся для дозвукового истечения с ГУ.

Тут же я не готов целовать крест какая скорость у нас в реальности получается на Входе.

Если есть оценка этой скорости, то нужно вместо ГУ Давление задаь Давление с Нормальной скоростью. И не забыть указать скорость.

Продожаем считать.

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


всё-таки я бы не отказывался от ГУ "давление"...

просто поднять его надо плавно - P=P(t)

поднять настолько плавно, чтоб не было ударной волны

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

Поднимал уже.

Попробую конечно же другую зависимость, но судя по поведению - не поможет.

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

Все, я зашел в тупик.

А что вообще за устройство? Можно по подробнее описать, как осуществляется подача газа в сопло?

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

А что вообще за устройство? Можно по подробнее описать, как осуществляется подача газа в сопло?

...посмотри файл-отчет, который я посылал в конце года (Сопло.pdf)...

Не думаю, что это сопло от моего принципиально отличается.

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

вообще, заинтриговали) люблю проблемки) сейчас ковырну задачку)

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

Проект прилагаю.

в этом проекте мне бросилось в глаза:

1) температура вакуума = 300К .... многовато. поставил 3К

2) вернул опорные значения в 300К@100000Па (соотв-но в пространстве заданы -99999,8Па@-297К)

3) свободный выход : заменил сверхзвук на статическое давление (равное НУ). статическое давление автоматом переключится на сверхзвук в нужных ячейках.

4) вход/выход - Статическое давление на входе заменил на Полное давление. Статику можно задать только на входе (т.е. если вдруг пойдет за обратной волной масса, то ничего хорошего не будет). И еще момент - статика там меняется (на мой взгляд) хитрее, чем полное давление, поэтому, по идее, полное - корректнее.

5) перекроил геометрию, т.к. выглядела странно...

попробую посчитать)

UPD по п.4.:

"Статическое давление на входе" и "Полное давление" - это дозвуковые ГУ.

Поэтому либо надо настолько плавно поднимать как давление, так и температуру, чтоб не было ударной волны совсем,

либо шаманить, меняя тип ГУ.

Вариантом тут является Нормальная скорость с давлением. Это ГУ допускает и вход и выход, на сверхзвуке близко к прибиванию скорости и давления, на дозвуке близко к прибиванию давления. В целом, логика зашита странная (не очевидная) - через массовую скорость, но в д.с. это может сыграть на руку.

в общем, пока оставил на нем считать. поставил маленькую скорость и рост давления.

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

4) вход/выход - Статическое давление на входе заменил на Полное давление.

Это можно делать, только если вход дозвуковой.

Все остальное я уже пробовал.

в общем, пока оставил на нем считать. поставил маленькую скорость и рост давления.

Я уже пробовал.

:)

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

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

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

Я уже пробовал.

А плотность таблично задавать?

Начальная плотность очень высокая получается.

М.б. стоит её понизить, задав таблично?

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

Я пока сделал финт ушами.

Поднял опорное давление до 10 Па, а шаг сверху ограничили 1e-9.

Картина течения другая. Теперь поток плавно доходит до конфузора и от него не отражается, а вдоль него образуется особая пристеночная зона течения.

Сейчас снял ограничение по шагу интегрирования и он вырос до 5e-8.

как получу струю, то опорное давлению опущу до исходного и посмотрим что получилось.

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

А разве можно менять опорное давление между шагами??

Мне говорили, что так нельзя делать.

а то, что шаг надо ограничивать - это понятный момент. я бы еще рекомендовал CFL 0.45 ставить.

и еще - там сужение резкое. в итоге у стенки вектор направлен вниз. Но ячейки вытянуты вдоль Х.

Поэтому считаю, что ячейки надо квадратные делать, чтоб шпв в ячейках был примерно одинаковый

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

Сетку я уже давно переделал. Она у меня внутри сопла везде квадратная - выше изображение приводил.

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

Если перестало считать, то значит звоночек.

P.S. Да нет, один хрен отражается и бьет во вход.

post-1864-1335445262.jpg

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

P.S. Да нет, один хрен отражается и бьет во вход.

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

post-26004-1335446343_thumb.png

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

Это нормально. Как раз раньше мы тут сетку адаптировали.

Сейчас решалка все это проглатывает.

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

Это нормально. Как раз раньше мы тут сетку адаптировали.

Сейчас решалка все это проглатывает.

фишка в том, что на эти несколько ячеек не подействовала, например, установка переменной.

в общем, они стремные...

пробую одновременно поднимать давление на входе и сразу во всем докритическом объеме (с помощью установки переменной)

т.е. и там, и там одна и та же формула стоит.

тогда, по идее, не будет волны.

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

Мы не знаем какой режим должен быть на входе.

Мальчики ждут девочек - показаний автора.

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

Мы не знаем какой режим должен быть на входе.

Мальчики ждут девочек - показаний автора.

формально, да.

реально нужен стационар выхода )

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

по идее, это баллон с клапаном.

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

Это нормально. Как раз раньше мы тут сетку адаптировали.

Сейчас решалка все это проглатывает.

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

скорость большая, ячейка маленькая => наименьший шаг по времени там.

И что - можно множитель шага поднимать?

И всё-таки сомнительно, то там не срабатывает установка переменной.

Занулить бы там модификатором скорость и считать дальше спокойно.

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

Если не срабатывает модификатор, то это повод для разбора.

В этом месте решалка должна прокашляться и поскакать дальше.

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

Если не срабатывает модификатор, то это повод для разбора.

В этом месте решалка должна прокашляться и поскакать дальше.

модификатор на таком месте вроде бы сработал 1 раз из трех.

И такие точки - они мигрируют и/или появляются в новых местах....(см.картинку) (скорость там упирается в ограничитель).

Правда, это дело происходило на 1м порядке точности. я этот момент не отследил в проекте, который скачал.

Сейчас поставил 2й порядок и 128 бит - лаги эти сдуло. Правда, они очень скоро и в более жестком развальном виде проявились в других местах ))

жесть )

post-26004-1335474118_thumb.png

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • boomeeeer
      @anykeyto Я не оператор. Настраивал работу этой функции на токарно-фрезерном, где нужно было позиционировать шпиндель, для нарезания пазов. Что такое G07 не знаю. Вот такой пример для фрезерного станка нашёл   G108 B0; переключение на ось B для управления положением; G68.2 X10 Y10 Z20 I0 J90 K0 Положение заготовки (10,10,20) Угол Эйлера для определения характеристической системы координат; G53.2 Управление направлением инструмента   G109 B0 - отключит режим
    • anykeyto
      Ага, спасибо... А в связке с G07.1 это как работает? Есть детали ,кольца, на которых надо делать по 4 паза. Сейчас программист написал прогу под шаровую фрезу. Там рабочие движеня происходят через одновременное смещение по оси Х и Z... фреза убивается буквально на третьем кольце, да и обработка получаестя фиговой.. На сколько я понимаю надо стремиться к тому, чтобы на 5ти осном станке фреза всегда была перпендикулярна касательной ли самой поверхности. Мне посоветовал как раз цилндрческую интерполяцию, когда ось С преобразуется в плоскость. В мануале к HNC описано эт криво, но в примерах программ как раз сначала включается CTOS и потом G07.1. но дальше как запрограммировать понимания пока нет...  P.S. хочу добавить, что на 5-снике раньше работал только как оператор по готовым программам  на Хенденхайне. А сейчас на новом предприятии на новой для меня стойке нужно все делать с нуля и самостоятельно. С одной стороны интересно, с другой стороны спросить не у кого, в инете инфы по стойке нет от слова совсем, уж тем более видео....  чем это отличается от просто G0 C180? Или это уже будет шпиндель позиционироваться, а не стол?
    • AlexKaz
      Различают жёсткие и гибкие ротора. В данном случае вал жёсткий. Для него соответствуюшая формула. 4000 об/мин = 66.67 Гц, это ни о чём, эффекты для такой частоты почти не скажутся на жёстком роторе. Но проблема в том, что в лопаточных машинах (а насос на скрине наверное лопаточная?) собственные кратны числу лопаток. Я не разбирался в вопросе чтобы что-то здесь советовать. Лучше не спешить и внимательно почитать про модальный в лопаточных. Заодно проверить, а какое нагружение идёт в деталях насоса чтобы корректно выполнить престресс.
    • ДОБРЯК
      Хотите еще и матрицу масс изменить. :=) И сколько килограмм воды будете присоединять?  Опыт правильного моделирования, опыт расчета таких конструкций накапливается десятилетиями. Маловероятно, что кто-то на форуме будет вам объяснять как правильно моделировать, какие КЭ использовать, какие граничные условия задавать и т. д. Вы же всегда уклоняетесь от разговора на эти темы, даже для простых моделей.  Даже если перейдете на Ансис и не разберётесь в том, что он считает это не поможет вам написать грамотный отчет и объяснить грамотному заказчику, что вы там насчитали. Считает расчетчик, а программа это инструмент. Это только в рекламе пишут, что купите программу трехмерного моделирования, сделайте трехмерную модель, потом нажмите одну кнопку и на обычном игровом компьютере все посчитаете. Только в жизни все немножко иначе...:=)
    • Jesse
      Вот что у Биргера: Центробежные силы и индуцированный гироскопический момент зависят от квадрата частоты вращения вала   Хотя уже сама критическая частота зависит только от момента инерции колеса (альфа, бэта, гамма - это прогибы, углы поворота и прочая шелуха), т.е. частота вращения в формулу не входит, а всё определяется только моментом инерции Jm Ну у него рассматривается простая одномассовая система (фактически грузик на пружинке). Думаю поискать публикации.. мб где то раскопали аналитическую зависимость первой и высшей СЧ от частоты вращения для сложных систем. Хотя если выяснится что зависит,  всё равно обоснование придумывать придётся.  Эхх...
    • Fedor
      https://www.litres.ru/book/v-n-faddeeva/vychislitelnye-metody-lineynoy-algebry-65999878/  Ну и классику жанра :) 
    • Jesse
      4000 об/мин при наружном диаметре колеса 380 мм. Думаете, на больших оборотах вилка всё равно может разойтись и могут появиться гироскопические эффекты? Даже несмотря на вышеприведенные массово-инерционные оценки и расположение опор?   да я влияние престресс (вращение) особо не увидел... Думал обычным модальником прогнать, Ланцошом (он пошустрее).   Ещё мои раздумья направлены в сторону присоединённой массы воды, что по идее должно занижать частоты. Хотя для такого массивного вала как у меня даже на это дело можно забить, наверно...
    • Handrusik
      Maxx2000 - благодарю Павлуха - это не фриланс Утомлённый солнцем - приятно, что меня помнят
    • Борман
      Я так понимаю, что обороты не особо большие. ДК - это чисто высокооборотистая тема.   Просто мысли вслух... Если вилка будет сильно расходится, то всегда будет точка пересечения любого луча с кривой BW. Причем точка будет на меньших оборотах по сравнению с обычным престресс-модальником. Какие дальше действия ? С чем тут нужно дополнительно бороться по сравнению с тем, с чем нужно было бы бороться по результатам обычного престресс-модальника Если вилка не сильно расходится (от кривой обычного престресс-мобальника), то смысла стоить строгую ДК большого нет.    
    • boomeeeer
      На сименсе параметрами настраивается многое, что связано с менеджером иснструмента. Вызов возможен как по номерам, так и по именам инструментов. Через равно вроде как по имени вызывается (Т="имя"), а по номеру можно было, например, Т1 M6, а первый корректор подтянулся бы автоматически. И не всегда М6 на смену. Производитель мог нагородить чего угодно. Параметры надо смотреть
×
×
  • Создать...