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

Ошибка при задании скорости на входе


Кирилл Александрович

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

Кирилл Александрович

Доброго времени суток всем участникам. При попытке задать скорость (изменяющуюся  по закону) на входе выдаётся ошибка. 

Задаются выражения: 

alpha=0.15 

z0=10 [m]

rho=1.225 [kg m^-3]

w0=230 [Pa]

k=(z/z0)^alpha 

u=sqrt(2*w0*k/rho)

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

Подскажите пожалуйста, заранее благодарен.

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


Скиньте кусок модели, куда эта скорость прикладывается, чтобы проще проверять было. И тема чего-то про дублировалась.

 

Изменено пользователем karachun
Ссылка на сообщение
Поделиться на других сайтах
Кирилл Александрович

Почему продублировала тема не знаю (

Кусок модели это как? Фото подойдёт ?

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

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

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

К сожалению на рабочем компьютере нет выхода в интернет ( постараюсь придумать что-нибудь.

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

Модель никак прикрепить не получается(.

Я был бы Вам очень признателен, если Вы поможете решить этот вопрос. 

Поясню суть задачи: имелся группа сооружений, ее необходимо рассчитать на ветровое воздействие. Данные по ветру из ГОСТа . Домен взят как прямоугольник , ГУ inlet на левой боковой грани, на нижней грани ( земля) и гранях которые относятся к  зданиям wall no slip, все остальные opening.  Если задавать скорость постоянную то задача считается 

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

Ну это все понятно, хорошо сейчас попробую разобраться.

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

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

Вроде как вы скорость под корень засунули. Размерность получается не та. Плотность постоянная, превратите динамический напор в скорость и вынесите ее из под корня.

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

У меня CFX ругался на названия выражений k и u, они зарезервированы. Попробуй заменить их на k1 и u1.

Все считается и не вылетает.

Еще вопрос - а должна ли быть k1 под корнем или это множитель скорости?

Снимок16.JPG

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

Почему продублировала тема не знаю

 

Бывает такой глюк при проблемах с интернетом..

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

У меня CFX ругался на названия выражений k и u, они зарезервированы. Попробуй заменить их на k1 и u1.

Все считается и не вылетает.

Еще вопрос - а должна ли быть k1 под корнем или это множитель скорости?

Снимок16.JPG

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

Большое спасибо @karachun за проявленный интерес и помощь!)

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

Не за что, только с файлом примером помощь приходит гораздо быстрее)

У меня в версии 19.0 название выражения "u" занято под скорость в направлении Х, а "k" это, скорее всего, постоянная Больцмана.

У меня CFX даже не позволил создать функцию с таким названием.

 

Снимок.JPG

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

Не за что, только с файлом примером помощь приходит гораздо быстрее)

У меня в версии 19.0 название выражения "u" занято под скорость в направлении Х, а "k" это, скорее всего, постоянная Больцмана.

У меня CFX даже не позволил создать функцию с таким названием.

 

Снимок.JPG

Я тоже заменил , все работает в 15.0.7

В чем была проблема ранее пока затрудняюсь сказать. 

Позже дополню тему. 

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

Модель никак прикрепить не получается(.

 

На форуме можно выкладывать файлы определенных типов.. Лучше всего архив..

 

Если размер большой, то на файлообменник..

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

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

Однако теперь встал вопрос в корректности поставленных ГУ для данной задачи. 

В 14.08.2018 в 15:59, Кирилл Александрович сказал:

Модель никак прикрепить не получается(.

Я был бы Вам очень признателен, если Вы поможете решить этот вопрос. 

Поясню суть задачи: имелся группа сооружений, ее необходимо рассчитать на ветровое воздействие. Данные по ветру из ГОСТа . Домен взят как прямоугольник , ГУ inlet на левой боковой грани, на нижней грани ( земля) и гранях которые относятся к  зданиям wall no slip, все остальные opening.  Если задавать скорость постоянную то задача считается 

 В таком случае аэродинамические коэффициенты не получаются близкими к значениям по СНиП. 

Вариации с ГУ особых изменений не дали.

Собственно вопрос: что делаю не так? 

Может есть рекомендации по заданию ГУ для подобных задач?

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

В таком случае аэродинамические коэффициенты не получаются близкими к значениям по СНиП.

Собственно вопрос: что делаю не так?

В CFD можно наделать кучу ошибок кроме неправильно заданных ГУ, и это нормально. На результат влияет качество и плотность сетки, выбор модели турбулентности, учет или игнорирование физических явлений (сжимаемая/несжимаемая среда, учет теплообмена - скорее всего это не ваш случай, но так в общем надо смотреть и на это). Можно расчетный домен сделать слишком маленьким.

Короче рекомендаций "как сделать правильно вообще" хватит на целую книгу. Скидывайте пример, тогда можно будет предметно обсуждать. Здания я не считал (знаю только, что там нужно использовать модель турбулентности SA), но постараюсь помочь.

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

Ссылка на сообщение
Поделиться на других сайтах
Кирилл Александрович
19 минут назад, karachun сказал:

В CFD можно наделать кучу ошибок кроме неправильно заданных ГУ, и это нормально. На результат влияет качество и плотность сетки, выбор модели турбулентности, учет или игнорирование физических явлений (сжимаемая/несжимаемая среда, учет теплообмена - скорее всего это не ваш случай, но так в общем надо смотреть и на это). Можно расчетный домен сделать слишком маленьким.

Короче рекомендаций "как сделать правильно вообще" хватит на целую книгу. Скидывайте пример, тогда можно будет предметно обсуждать. Здания я не считал (знаю только, что там нужно использовать модель турбулентности SA), но постараюсь помочь.

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

Хм. Тогда это тянет на новую тему. 

Спасибо за ответ. 

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

Тут еще такое дело - в СНиПе или ГОСТе коэффициент это просто одна цифра, какими соображениями руководствовались авторы Вам не известно. Возможно цифра взята в запас, или усреднена для нескольких типов конструкций. Если хотите проверить правильность расчета, то поищите статьи где описан простой эксперимент, возможно продут макет одного отдельно стоящего здания. Есть поля давлений, скоростей. Возможно есть статьи, в которых уже есть сравнение эксперимента и CFD - тогда будет ясно какие подходы использовать и что можно ожидать от расчета. И начинать лучше всего с простых моделей и добившись хорошего совпадения с экспериментом переходить на сложную модель.

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

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

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

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

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

  Разрешено не более 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 на смену. Производитель мог нагородить чего угодно. Параметры надо смотреть
×
×
  • Создать...