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

Моделирование движения резиновой мембраны под действием течения воды


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

Здравствуйте.

Есть задача смоделировать и рассчитать движение ленты (мембраны) под действием течения жидкости. + замерить расходы

post-47674-0-82458500-1449870327_thumb.png

post-47674-0-87010800-1449870332_thumb.png

 

Начальные условия таковы:

1. С незакрытой грани "куба" подсоединен резервуар с водой высотой 1.6м. Т.е. давление на дне как вы понимаете 1.6 м/в. ст.

2. Отверстие регулятор полностью открывается, после чего начинается движение лены

 

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

 

Опыта мало, но прочитав литературы и форумы стало ясно что такая задача решается с помощью 2-way Fluid-Structure Interaction (FSI). 

Стал рассматривать какие же программные комплексы предоставляют такую возможность. И остановился на Ansys и COMSOL Multiphysics.

 

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

 

Поэтому я начал поиск альтернативы, но контекстная реклама не заставила меня долго маяться=) 

И взор пал на Comsol. 

Но и тут неудача следующего характера:

Division by zero.
	- Function: /
Failed to evaluate expression.
	- Expression: nojac(1/sqrt((6*(comp1.fsi.rho^2))*emetric_spatial(comp1.u_fluid-d(x,TIME),comp1.v_fluid-d(y,TIME),comp1.w_fluid-d(z,TIME))))
Failed to evaluate variable.
	- Variable: comp1.fsi.tau_mCross
Failed to evaluate variable.
	- Variable: comp1.fsi.vhNS_quad
Failed to evaluate expression.
	- Expression: nojac((1.4142135623730951^(1-shapeorder(comp1.u_fluid)))*max(comp1.fsi.vhNS_lin-comp1.fsi.vhNS_quad,0))
Failed to evaluate variable.
	- Variable: comp1.fsi.vhNS
Failed to evaluate temporary symbolic derivative variable.
	- Variable: comp1.fsi.crosswindns@VDN${test@15}@VDN$comp1.pz
Failed to evaluate expression.
	- Expression: ((comp1.fsi.vhNS*comp1.fsi.rho)*(-(comp1.fsi.betaT*comp1.fsi.gijxz)))*dvol_spatial
Failed to evaluate Jacobian of expression.
	- Expression: ((comp1.fsi.vhNS*comp1.fsi.rho)*((((-comp1.fsi.gijxx)*(comp1.fsi.betaT*comp1.px))-(comp1.fsi.gijxy*(comp1.fsi.betaT*comp1.py)))-(comp1.fsi.gijxz*(comp1.fsi.betaT*comp1.pz))))*dvol_spatial
Failed to evaluate Jacobian of expression.
	- Expression: (comp1.fsi.crosswindns)*(dvol_spatial)

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

 

Поэтому прошу помочь люди добрые!

 

Рад буду услышать любые ответы

 

P.s. прикладываю файл Comsol

model.rar

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


тема перемещена.

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

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

тема перемещена.

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

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

Может кто поможет решить данную задачу в Ansys? Или лучше создать новую тему с описанием моих нерадивых действий в Ansys, которые приводят меня к ошибке?

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

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

Может кто поможет решить данную задачу в Ansys? Или лучше создать новую тему с описанием моих нерадивых действий в Ansys, которые приводят меня к ошибке?

по ансис лучше отдельную тему

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

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

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

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

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

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

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

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

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

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

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




  • Сообщения

    • gudstartup
      как перевести шпиндель снова первоначальный режим управления скоростью когда ему пофиг число импульсов на оборот? резьбу пока резать не будут а в подпрограмму смены можно что нить добавить или в реализацию m19  чтобы чпу из режима контроля позиционного кодера вышло?
    • gudstartup
    • mactepnew
      Вот я тоже так думал, пока не столкнулся с обратным) Как еще столкнусь, пришлю скриншот. Хаос еще в том, что порядок в дереве построения весьма условный. операция, которая стоит перед групой, имеет бОльший timestamp(пардон, не знаю как это в переводе), чем эта группа, и выполняется позже. При этом все операции в группе имеют timestamp на пару сотен меньше и выполняются раньше, чем две группы, стоящие перед ней. Логика?  Если в эскизе все элементы определены, то можно поставить еще один размер (хочешь переопределить эскиз - да пожалста!), изменить его - и NX самостоятельно выберет, какой размер или ограничение отпустить и весь эскиз летит к чертям. Это ли не хаос?   Назовете хоть один CAD, которым можно пользоваться так? Любой из встреченных мной имеет много неочевидных нюансов, которые конечно не проявляются при выполнении университетского или базового курса. В Creo я на протяжении 5 лет находил всё новые фишки, которые не прописаны в документации. Огромная благодарность Ruslan за пояснения тонкостей. Даже Компас и Solidworks, которые считаются самыми дружественными, оказалось непросто обуздать, когда проходил собеседования и делал тестовые. А ведь университетский курс знал на отлично!   Возвращаясь к NX - вот сегодняшний пример. Утром деталь была оранжевая, единая. Изменил пару уклонов и толщин, стала желтая с элементами оранжевого и частично твердотельной, частично листовой. Я даже не знаю, как сформулировать вопрос в поддержку. 
    • Mikhail_Podgornyj
      Добрый день, кто может поделиться SKF Bearing V1.6 (ANSYS App)? Заранее благодарен!
    • kkk
      Короче, как сделать автоматом не нашел. Раскрасил нужные грани нужными цветами. Затем на вкладке с цветом модели можно выбирать нужный цвет, тогда в модели выбираются грани, покрашенные этим цветом. Можно еще сохранить выбор этих граней, но смысла не вижу, т.к. все легко выбирается по цвету. Затем либо в инструменте Измерить смотрим площадь и как-то ее используем, либо можно макросом (что выкладывал Kir95 ) внести площадь в свойства (на вкладку Конфигурации). Можно еще внести выбранные грани в датчики, но как использовать их в свойствах не нашел.    Макрос прикрепляю.     Площадь выбранных граней в свойства.swp
    • IgnatvKudrin
      Для нашего производства требуются постоянно инструментальные пружины для штампов. Сейчас заказы обрабатываются довольно долго - до 2-х месяцев. А у нас сроки на штампы 1-2 месяца с испытанием, доработкой и т.п. Хотелось бы найти поставщика, у которого уже есть склад с данными пружинами или кто может доставить за 1 месяц.   Пружины стандартные по ISO, но на всякий случай приложил каталог. Вот список пружин, которые нужны очень срочно:   LR 32 x 44 - 1 шт. LR 25 x 51 - 2 шт. LB 63 x 89 - 7 шт. LB 32 x 102 - 10 шт. LB 32 x 89 - 20 шт. LB 32 x 76 - 12 шт. LB 32 x 51 - 6 шт. LB 25 x 127 - 9 шт.   Готов выслушать все предложения. Возможно даже есть те, кто сможет изготовить подобные пружины с такими же параметрами. 20231101.pdf
    • bubblyk
      Знатоки форума, подскажите, как смоделировать простую (на первый взгляд) задачу термомеханики в упругой постановке.   Есть брусок правильный формы ака кирпич, в котором генерируется НЕРАВНОМЕРНОЕ удельное тепловыделение (HGEN). Задняя и передняя поверхности конвективно охлаждаются с известными (но различными) коэффициентами теплоотдачи и температурами жидкости/газа. Свойства кирпича (теплофизика + механика) принимаются зависящими от температуры. Постановка - стационарная Собственно вопрос : как грамотно смоделировать температурное поле в кирпиче с учетом температурного расширения?   Исходный код приложил.     Надеюсь на коллективный разум и опыт! EPTH2D.mac p.s.   Если есть опыт решения этой задачи в Ansys WB, тоже можно сюда приложить командные вставки.
    • OutXL
      Это было давно использовал NC-Link  и данную инструкцию. The Maho.docx
    • Mikhail_Podgornyj
      Добрый день, кто может поделиться SKF Bearing V1.6 (ANSYS App)? Заранее благодарен!
    • alexqwq10
      В данных корпусах ничего лишнего быть не должно. Максимум шпильки заземления, в количестве предусмотренном для конкретного корпуса. Рассматривал очень много подобных корпусов разных производителей и данная проблема актуальна для подавляющего большинства.
×
×
  • Создать...