Контакты frictional выдают ошибку при статичесом расчете

[rasta89 113]

Всем привет.

Предлагаю к рассмотрению мою проблему с контактами.

Задача такая, есть устройство с резиной внутри трубы, резина (3 шт, с кольцами между) сжимается под давлением и заполняет межтрубное пространство, цель - увидеть как деформируется резина, в идеале как она затекает в пространстве между трубами. Решил начать статического расчета, заданы фиксации, между резиной и кольцами задано frictional контакты, с одной стороны на набор резинок давит усилие. Резине задан свой материал.

 

труба зафиксирована с двух сторон

нижний конус зафиксирован чтобы резина с кольцами не ушли вниз

верхний конус, резина и кольца между резинками -frictional

труба в которой все это стоит зафиксирована с двух сторон

 

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

 

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

выложил проект, если кому интересно посмотрите

https://disk.yandex.ru/d/kWW0uScLdYysuA

[rasta89 113]

3 часа назад, karachun сказал:

точки в студию, если речь о Curve Fitting

речь о том что ни смотря ни на что расчет все равно дальше ваших 8 мм не продвигается

[karachun 2 038]

 

 

@rasta89 Что это за хрень? Почему такие большие величины у Strain? Если мы задаем относительные деформации в мм/мм то единица соответствует 100%. А график на сайте заканчивается на 200%. Почему тогда Strain в таблице равен 50 тысячам?

 

Почему Incompressibility Parameter не задан? Если он нулевой то материал абсолютно несжимаемый а это может вызвать численные проблемы. Лучше взять хотя бы значение из моей резины.

 

Вот как так можно тупить? Мы же все это обговариваем уже два месяца и тут такие банальные ошибки.

 

Короче все, я рассказал все что хотел рассказать, и дал ссылки на те справочные материалы которые знаю. пусть вам дальше @Jesse, @ak762 и @none помогают. А с меня хватит.

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

А еще говорите что не самый глупый человек на форуме. Вот ей богу. даже @Kotikov бы справился быстрее.

[rasta89 113]

2 часа назад, karachun сказал:

 

 

@rasta89 Что это за хрень? Почему такие большие величины у Strain? Если мы задаем относительные деформации в мм/мм то единица соответствует 100%. А график на сайте заканчивается на 200%. Почему тогда Strain в таблице равен 50 тысячам?

 

Почему Incompressibility Parameter не задан? Если он нулевой то материал абсолютно несжимаемый а это может вызвать численные проблемы. Лучше взять хотя бы значение из моей резины.

 

Вот как так можно тупить? Мы же все это обговариваем уже два месяца и тут такие банальные ошибки.

 

Короче все, я рассказал все что хотел рассказать, и дал ссылки на те справочные материалы которые знаю. пусть вам дальше @Jesse, @ak762 и @none помогают. А с меня хватит.

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

А еще говорите что не самый глупый человек на форуме. Вот ей богу. даже @Kotikov бы справился быстрее.

ну, с этими скринами я конечно профанулся, не заметил что импортировалось без разделительной запятой, про Incompressibility Parameter тоже не доглядел

но даже учитывая эти параметры правильными, все равно расчет пока что рушится, перебрал все что можно, контакты, сетки, воздействия

[karachun 2 038]

@rasta89 Расчету это не поможет но надо еще точку (0; 0) добавить.

[rasta89 113]

Пока лучше результата не вышло, происходит излом сетки, сейчас она адаптивная, после излома расчет рушится, 14мм

[ak762 307]

On 8/29/2021 at 8:37 PM, karachun said:

Расчету это не поможет но надо еще точку (0; 0) добавить.

а в модели материала Mooney-Rivlin/Ogden разве график не может представлять компрессию в минусовой области т.е может надо не только с нуля начинать?

[karachun 2 038]

@ak762 Можно и нужно но данных нет. Предоставьте данные для сжатия тогда и начнем.

[rasta89 113]

 

[soklakov 1 475]

красненько - там контакт настроен?

желтенько - как будто локальная потеря устойчивости.  чтобы отмесТи эту версию, можно потестить non-linear stabilization... или поменять задачку, чтобы исключить баклинг.

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

[rasta89 113]

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

Там где вы отмечали красным там стоит фрикционный контакт, он там есть.

 

С контактами перетестил тоже кучу вариантов от стандартных до вариаций: formulation - pure penalty, detection method - nodal-normal to target, выставлял normal stufines, behaviour -assymmetric/symmetric и т д. Я наверно все перепробовал, измельчал сетку до 0,4 мм, сетки перетестил все, перестроение сетки тоже не дает особого результата.

 

Задался даже задачей подогнать тупо резину (константы, законы сжатия) чтобы у нее уровень сжатия был максимальный, хотя бы чтобы увидеть визуально ее сжатие (без соответствия там ее настоящим свойствам), не вышло. Сделал элементарный цилиндр и сжимал его, и вот однажды при помощи polynormal он начал сжиматься на 80-90% удача подумал я, спроецировал все на необходимую модель, и снова нет, сетка разрушается.

 

Запросил у поставщика максимально возможные данные по резине, буду ждать ответ.

[AlexArt 278]

Прошел месяц... осевую постановку вы так и не осилили. Перестроение сетки не вижу.

[rasta89 113]

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

Прошел месяц... осевую постановку вы так и не осилили. Перестроение сетки не вижу.

Картина не меняется, позже выложу с перестроением.

Осевую постановку - в смысле осесимметричная постановка в 2Д? Как уж не осилил, осилил и тестил.

[soklakov 1 475]

 

опять эта шляпа)

гляньте, кстати, напряжения и невязки, а то эти поля перемещенеий совсем неинформативны.

 

 

1 час назад, rasta89 сказал:

Как уж не осилил, осилил и тестил.

а чего картинки не с нее?

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

[rasta89 113]

7 часов назад, soklakov сказал:

 

опять эта шляпа)

гляньте, кстати, напряжения и невязки, а то эти поля перемещенеий совсем неинформативны.

 

 

а чего картинки не с нее?

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

 

Цитата

а чего картинки не с нее?

Вроде версия сохранилось, загружу скрины, проблемы те же. Решил делать сразу в "трушном" виде. Как-то было такое что добился нужной сходимости в плоской задаче а при переносе уже на рабочую модель проблемы не уходили. Поэтому делал все в оригинале как говорится.

 

Цитата

опять эта шляпа)

Она и не решается никак(

Кстати она (сетка в месте излома) ведёт себя по-разному с разным законом сжатия и константами. Такое ощущение что дело в материале а не в контактах, но опять таки, хотя бы одна из 100 попыток должна была дать результат

Изменено 14 сентября 2021 пользователем rasta89

[soklakov 1 475]

10 часов назад, rasta89 сказал:

Решил делать сразу в "трушном" виде.

а зря

10 часов назад, rasta89 сказал:

Такое ощущение что дело в материале а не в контактах, но опять таки, хотя бы одна из 100 попыток должна была дать результат

нет. хоть миллион попыток, если все они бестолоковые, не дадут результата.

решить задачу случайным перебором настроек- не вариант.

[rasta89 113]

5 часов назад, soklakov сказал:

а зря

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

 

Пока я разбираюсь отклонимся немного от темы.

Такой вопрос, можно ли решить задачу по герметичности в данном случае не прибегая к MAPDL?

Например используя предстресс полученных спрессованных резинок.

Например это было бы так, спрессовали резинки, они закрыли промежуток, затем на спрессованные резинки даем давление, и в результатах смотрим на контакте резины с железом в contact status, если статус контакта «closed», без зазора, значит герметично.

[soklakov 1 475]

4 минуты назад, rasta89 сказал:

Например это было бы так, спрессовали резинки, они закрыли промежуток, затем на спрессованные резинки даем давление, и в результатах смотрим на контакте резины с железом в contact status, если статус контакта «closed», без зазора, значит герметично.

только смотреть надо давление в контакте, а не статус. если оно выше давления в жидкости - то герметично. лучше с запасом 2.

[karachun 2 038]

1 час назад, rasta89 сказал:

не прибегая к MAPDL?

Вот по этой ссылке есть рабочий пример с небольшой концентрацией APDL, порядка 10-15%.

 

13.09.2021 в 22:58, soklakov сказал:

опять эта шляпа)

Предположу что решение может улучшить включение mixed u-P Formulation?

Ну и другие опции, наподобие Reduced Integration. Что именно нужно я так и не скажу. у меня в Настране оно работает как-то само, безо всех этих дополнительных костылей, ну или костыли уже зашиты в решатель)

[soklakov 1 475]

8 минут назад, karachun сказал:

Предположу что решение может улучшить включение mixed u-P Formulation?

по идее она автоматически включилась из-за гиперэластика. можно проверить в аутпуте

[karachun 2 038]

35 минут назад, karachun сказал:

Вот по этой ссылке есть рабочий пример с небольшой концентрацией APDL, порядка 10-15%.

А ссылку не приложил.

https://cccp3d.ru/topic/86970-решение-контактной-задачи-с-fluid-pressure-penetration-в-воркбенче/