Падение клина в жидкость [2D]

[Dem_SMK 0]

Доброго времени суток.

 

Задача следующая: абсолютно жесткий клин проникает в покоящуюся жидкость с постоянной скоростью. Сетка клина и жидкости разбита 1 к 1, используется 11 формулировка (ALE мульти-материал). Материал жидкости - Fluid_Elastic, вместо воздуха - mat_vacuum. Клин разбит shell-элементами. По ГУ: левая кромка по симметрии, остальные три кромки ГУ с не отражением акустических волн. Клин может двигаться только вертикально. Вроде бы все хорошо, но получается два варианта:

 

1. При задании модуля объемной упругости К=2,2*10^7 Па и коэф. штрафа PFAC = 0.1, давление выходит на уровень 0,45 МПа.

2. При задании модуля объемной упругости К=2,2*10^9 Па и коэф. штрафа PFAC = 0.001, давление выходит на уровень 0,75 МПа. Но после полного погружения клина уходит в ноль.

Аналитика дает около 0,7 МПа. Основной вопрос почему давление падает до 0 во втором случае? Может кто-то наведет на правильный путь мыслей. Где-то промах с подбором настроек контакта или в постановке задачи?

Изменено 11 марта 2016 пользователем Dem_SMK

[vl 327]

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

вот поэтому и 0

[Dem_SMK 0]

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