Настройка Contact_Interior в LS-Dyna

[AndreySh 0]

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

Рассчитываю высокоскоростное взаимодействие составного ударника (оболочка металл MAT_PLASTIC_KINEMATIC + упругий материал в середине MAT_ELASTIC) с преградой. Осесимметричная задача, 2Д. Скорость 1200 м/с. Все идет нормально: преграда разрушается, ударник деформируется.

Однако в упругом материале возникают значительные возмущения сетки и расчет вылетает с ошибкой "Error 40509 (SOL+509) negative volume in solid element". Ниже картинка проблемного места. 

Поиск решения проблемы привел к необходимости добавления нового контакта для упругого тела - CONTACT_INTERIOR, который не дает узлам сетки приближаться друг к другу менее определенного расстояния. Применил данный контакт к упругому телу, но результата нет. (В самой карточке контакта есть единственный параметр - выбор номера PART ID - картинку приложил).

Подскажите, пожалуйста, что и где ещё нужно задать, чтобы данный контакт заработал корректно?

Или может есть другой способ решения проблемы "negative in solid element"?

Заранее спасибо.

 

P.s.

Увеличение значений HOURGLASS и изменение его варианта никак на данную проблему не влияют.

 

[ANT0N1DZE 178]

Ремешинг пробовали?

[AndreySh 0]

Да. Сетку упругого тела уменьшал в 2, 4 и даже 5 раз. Результат тот же.

 

[Victoria 104]

Скрытый текст

 

Имелось ввиду adaptive remeshing в самой программе (в процессе вычислений)

Оператор *CONTROL_ADAPTIVE

Изменено 1 ноября 2019 пользователем Victoria

[AndreySh 0]

@Victoria

Теперь понятно. Спасибо за разъяснение.

Нет, *CONTROL_ADAPTIVE  не пробовал. Интересный вариант - посмотрю в эту сторону.

 

С другой стороны хотелось бы понять как пользоваться именно CONTACT_INTERIOR, т.к. это видится, как более простой способ.

[Victoria 104]

При сильной деформации элементов CONTACT_INTERIOR вряд-ли поможет.

А вот адаптивная перестройка элементов дает шанс на успех.

И реализуется для 2D задач достаточно просто

Изменено 2 ноября 2019 пользователем Victoria

[AndreySh 0]

@Victoria

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

Поделитесь ссылками, где можно ознакомиться с реализацией "достаточно просто"?

 

4 часа назад, Victoria сказал:

При сильной деформации элементов CONTACT_INTERIOR вряд-ли поможет.

Он идеально подходит для подобных задач, т.к.:

"Назначение: задавать внутренний контакт для элементов-брусков из пеноматериала. Когда пеноматериалы сжимаются при высоком давлении, объемные элементы, используемые для дискретизации этих материалов, нередко инвертируют, порождая отрицательные объемы элемента и ошибочное прекращение счета. Чтобы предотвратить инвертирование элементов, имеет смысл рассматривать контакты внутри пеноматериала между внутренними поверхностями, образующими грани объемных элементов."

[Victoria 104]

=============================================

*CONTROL_ADAPTIVE
         2      0.05         8         2
*PART
Part          1 for Mat         1 and Elem Type         1
         1         1         1         1         0         0         0
*PART
Part          2 for Mat         2 and Elem Type         1
         2         1         2         2         0         0         2

=============================================

Это все, что нужно для запуска адаптивной перестройки сетки для 2-х частей

=============================================

Примеры см ниже

=

penetration2d-adaptive.k

pen-adapt.k

metal_forming.k

=

Из Вашего описания контакта следует, что он предназначен для ЭЛЕМЕНТОВ-БРУСКОВ (объемных элементов)

А где-же они в этой задаче?

Изменено 3 ноября 2019 пользователем Victoria

[AndreySh 0]

Спасибо за ответ и отдельное большое спасибо за примеры. Теперь все стало намного более понятно.

 

4 часа назад, Victoria сказал:

А где-же они в этой задаче?

Задача представлена двумерная и, естественно, объемных элементов нет. В тоже время я не нашел более подробного описания данного контакта и, возможно, есть реализация для 2Д случая.

 

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

 

И все же, если у кого-то есть опыт использования Contact_Interior для 2Д или 3Д варианта поделитесь пожалуйста.

[AlexKaz 1 822]

@AndreySh , IMHO, врубите hourglass, а то выше чётко видно что элементы расползаются. А также демпфирование и вязкость.

Изменено 5 ноября 2019 пользователем AlexKaz

[Victoria 104]

В 01.11.2019 в 14:05, AndreySh сказал:

P.s.

Увеличение значений HOURGLASS и изменение его варианта никак на данную проблему не влияют.

 

[AlexKaz 1 822]

В *part может быть отключено.

[fzbm 44]

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

[AndreySh 0]

Странно не приходили сообщения о новых ответах.

 

В 05.11.2019 в 07:09, AlexKaz сказал:

В *part может быть отключено.

Данный элемент включен в *part.

 

В 11.11.2019 в 10:55, fzbm сказал:

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

Забыл упомянуть. При соударении со скоростью 800 м/с подобные артефакты отсутствуют.

[fzbm 44]

20 часов назад, AndreySh сказал:

Забыл упомянуть. При соударении со скоростью 800 м/с подобные артефакты отсутствуют.

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