Отладка осесимметричной задачи при использовании SPH

[Kracozebr 4]

Разрабатываю осесимметричную задачу пробития стальной пластины стальным пробойником. В control_SPH параметр IDIM ставлю равным -2 (что означает формулировку для осесимметричной задачи). При запуске на расчет выдает ошибку Error 30282 <INI+282> SPH initianilization SPH element number 57395 has a negative mass.

Захожу в keyword manager и смотрю массу узла, она положительная все вроде бы нормально. Удалил узел, в котором масса отрицательная запускаю на расчет, уже в другом узле масса отрицательная.

Кто знает в в каком направлении копать?

Что можно почитать по этому поводу? Вообще материала по SPH в  LS-Dyna не очень много.

[fzbm 44]

А пробовали решать в задачу в полной постановке?

[Kracozebr 4]

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

[fzbm 44]

Да спасибо. Как то я делал ALE 2D там тоже важно правильно расположить относительно осей (вроде в положительном диапазоне...)