cfs post вектор напряжений

[dolivanov 7]

43 минуты назад, soklakov сказал:

Цепляем бляшку на контакт к стенке. Потом, когда-нибудь, не сразу, можно сделать контакт CZM. Но пока просто смотрим напряжения в контакте. Когда превышены усилия адгезии в большой области - бляшка сорвалась.

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

Вы без касательных-то посчитали?

Я считаю без бляшки, пока в постановке с жесткими стенками. Со стенкой и бляшкой сложно всё, так как непонятно, где ее свойства брать, надо КТ томограмму анализировать, это отдельная задача. Или натурные эксперименты проводить, но материала и нет, и он будет трупный, свойства будут неидивидуальными. 

 

Что значит без касательных посчитали?) С ними так же, как и с давлением. 

 

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

 

47 минут назад, ДОБРЯК сказал:

Не получится. 

Напряжение в точке (в узле) это не 3-х мерный вектор, как сила или перемещение, а тензор 2-го ранга. 

Матрица 3х3. Шесть независимых величин. Шесть, а не три. )

Я это понимаю, но я речь веду о напряжении в точке, которая находится на площадке с нормалью. По-другому - есть поверхность, я на этой поверхности хочу оценить вектор напряжения, поверхность в своих точках имеет вектор нормали, значит площадки в этих точках уже определены и на них можно считать вектор напряжений. Не так? 

Даже задачи есть для студентов, типа "определить вектор напряжения в точке на площадке с нормалью n". 

 

Изменено 9 сентября 2019 пользователем dolivanov

[piden 2 874]

5 hours ago, dolivanov said:

делаю сумму векторов давления и Wall Shear (переменная PressurePlusWSS)  - не работает. Не вычисляет и не отрисовывает.

 

Делал в Ансис 19.2.

 

Вот архив с проектом - Fsapr_Tube.zip

Таким методом у меня отрисовывает:

 

Если задавать через pressVec X + Wall Shear X ... - то тоже не рисует. Скорее всего User Variables должны выражаться через некие базовые величины, а не через другие User Variables, иначе придется столько всяких пользовательских переменных повычислять перед отображением!)

 

Могу скинуть .cst, но думаю, разберетесь.

 

Нет времени вникать во все, что было выше, но вот пару советов:

1. @soklakov прав на счет интеграла -- если вам нужен результирующий вектор, то берите Force. Там уже учтен интеграл по давлению + касательная составляющая от wall shear. Фактически, напр. Force X == areaInt(Pressure * Normal X + Wall Shear X).

2. Не используйте ave() -- это алгебраическое среднее, без учета площади каждой элементарной площадки. Используйте areaAve() для получения "правильного" среднего.

3. Пока писал про Force, подумал, насколько правомерно использовать Total Pressure вместо Pressure. Скорость на стенке должна быть = 0, но численно не всегда-то так получается . И тогда в случае использования Total Pressure у нас динамическая составляющая (касательная) вдруг учитывается как нормальная.

Сравните эпюры Total Pressure и Pressure на стенке в районе inlet'a, где профиль скорости потока еще не успел сформироваться.

0. Не тратьте время и нервы на @ДОХЛЯК'а. 

 

Изменено 9 сентября 2019 пользователем piden

[dolivanov 7]

14 минуты назад, piden сказал:

Таким методом у меня отрисовывает:

 

Если задавать через pressVec X + Wall Shear X ... - то тоже не рисует. Скорее всего User Variables должны выражаться через некие базовые величины, а не через другие User Variables, иначе придется столько всяких пользовательских переменных повычислять перед отображением!)

 

Могу скинуть .cst, но думаю, разберетесь.

 

Нет времени вникать во все, что было выше, но вот пару советов:

1. @soklakov прав на счет интеграла -- если вам нужен результирующий вектор, то берите Force. Там уже учтен интеграл по давлению + касательная составляющая от wall shear. Фактически, напр. Force X == areaInt(Pressure * Normal X + Wall Shear X).

2. Не используйте ave() -- это алгебраическое среднее, без учета площади каждой элементарной площадки. Используйте areaAve() для получения "правильного" среднего.

3. Пока писал про Force, подумал, насколько правомерно использовать Total Pressure вместо Pressure. Скорость на стенке должна быть = 0, но численно не всегда-то так получается . И тогда в случае использования Total Pressure у нас динамическая составляющая (касательная) вдруг учитывается как нормальная.

Сравните эпюры Total Pressure и Pressure на стенке в районе inlet'a, где профиль скорости потока еще не успел сформироваться.

0. Не тратьте время и нервы на @ДОХЛЯК'а. 

 

Спасибо.

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

Тоже подумал уже дома, что variables не могут выражаться через новые переменные.

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

Давление сравнивал вроде, посмотрю ещё завтра на работе. 

Про виды средних почитаю, спасибо.

[ДОБРЯК 602]

8 часов назад, dolivanov сказал:

По-другому - есть поверхность, я на этой поверхности хочу оценить вектор напряжения

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

Изменено 10 сентября 2019 пользователем ДОБРЯК

[piden 2 874]

5 hours ago, ДОБРЯК said:

И что вы сможете оценить в результате?

А я думал, ты узнаешь свое любимое векторное давление 

Он нагрузку хочет оценить. "Ты даже этого понять не можешь"...

[ДОБРЯК 602]

24 минуты назад, piden сказал:

А я думал, ты узнаешь свое любимое векторное давление 

Он нагрузку хочет оценить. "Ты даже этого понять не можешь"...

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

Теперь понятно, что ты не понимаешь, что механическое напряжение тензор 2-го ранга. )

И если бы было 6 величин, ты бы все 6 сложил, чтобы нагрузку оценить. ))

 

Изменено 10 сентября 2019 пользователем ДОБРЯК

[piden 2 874]

2 minutes ago, ДОБРЯК said:

Но там еще есть и касательное напряжение для этой же точки.

Расскажи, как он в CFD касательные напряжения получил. Наверняка это увлекательная история.

[ДОБРЯК 602]

9 минут назад, piden сказал:

Расскажи, как он в CFD касательные напряжения получил. Наверняка это увлекательная история.

Так ты у него и спроси. Почему ты сразу не спросил? 

Объяснял что-то ), а потом выясняется что ты не понимаешь откуда ТС касательные напряжения получил. )

[dolivanov 7]

Вот, что мне нужно оценить.

[dolivanov 7]

Кстати, в гидродинамике для вязкой жидкости, как мне известно, тензор напряжений имеет 4 независимые компоненты, а не 6.