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

[supwgeneral 1]

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

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

 

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

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

Изменено 11 декабря 2015 пользователем supwgeneral

[Bonusfrag 104]

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

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

[supwgeneral 1]

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

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

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

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

Изменено 13 декабря 2015 пользователем supwgeneral

[Bonusfrag 104]

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

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

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