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

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


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

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

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

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

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

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 пользователей

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




  • Сообщения

    • Рома калужский
      это код подключения шпинделя как ось "С" но при этом он не перемещает револьверку, т.е. не работает по программе
    • Горыныч
      Еще совет: строить нужно только половину, если модель симметричная, тогда Солид меньше "косячит". И соглашусь с коллегами, идеальный вариант для гладкого построения - два сечения и две направляющие кривые, т.е. 4 управляющие кривые (кромки). При необходимости можно добавлять не только промежуточное сечение, но и направление. Не нужно пятаться построить всю поверхность одним фичером.  Пусть поверхности будут лоскутами, но соединяются по касанию или кривизне. 
    • ДОБРЯК
      Откуда появится смещение центра масс колеса, если колесо сбалансировано.  Жесткий сбалансированный вал, это условие данной задачи. На 4000 оборотах эти силы будут минимальны. Это разговор ни о чем в данной задаче. :=) Но даже если вы абсолютно точно определите эти частоты, с учетом всех сил. Что дальше с этими частотами делать? На основании чего делать вывод о работоспособности изделия? Это в учебной работе получили число, работу приняли и сдали в архив. Какая практическая ценность всех этих уточнений в данном случае?:=)  
    • Борман
      Если еще не поздно... я и сам погрузился бы. Физика такая, для вала с одним колесом... F - это понятно, центробежное смещение центра масс колеса M - момент от центробежных сил, который не появится, если колесо не повернется при изгибе вала. Обрати внимание, что BW-момент выворачивает ось в сторону действия силы F, т.е. при обратной прецессии они действуют сообща, и в этом опасность.   Отсюда будут разные собственные частоты при прямой и обратной прецессиях.  Теперь осталось понять, что такое прямая и обратная прецессии :), и почему момент при них направлен в разные стороны. По первой части вопроса см. мою картинку в этой теме от 27 августа 2013. По второй части вопроса - нужно включать воображение и усердно рисовать на бумажке. Мне пока сложно.    Физика именно такая. Нужно как то оценить величину этого момента и приложить его в софте, который не умеет ДК. А с тестовыми значениями можно и так порешать, и на выходе должна быть вилка. Т.е. F - задает типа "биссектрису" вилки (престресс-модальник), а момент - угол вилки.  
    • 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/  Ну и классику жанра :) 
×
×
  • Создать...