Численный анализ разлета осколочной оболочки с построением пространственного распределения осколков для АО НПО Базальт
Основная цель расчета состояла в определении пространственного распределения осколков по 5-градусным угловым зонам, необходимого для оценки могущества осколочного действия 82-мм изделия заказчика, на основе трехмерного численного моделирования с применением модели естественного дробления корпусов изделия из сталистого чугуна и специальной стали.
Данная работа включала следующие этапы:
- Создание геометрической модели изделия на основе открытых данных с упрощением конструкции (defeaturing).
- Построение сеточной модели (дискретизация расчетной области) для изделия.
- Задание физических моделей для созданной расчетной области.
- Проведение серии расчетов.
- Обработка и анализ результатов.
Решение
В качестве расчетной области рассматривалась внутренняя полость изделия, заполненная ВВ, и окружающая область воздуха радиусом в 2 раза больше калибра изделия. Расчет в газовой области велся в переменных Эйлера, для дробящегося корпуса изделия – в переменных Лагранжа. .Стыковка сеток осуществлялась с помощью подвижного граничного условия в среде LS-DYNA Исходной геометрией был 2D-чертеж, построенный в программе «Компас». Через нейтральный формат IGES данный чертеж был импортирован в контуры 2D-геометрии в препроцессоре ANSYS LS-DYNA, см.рис.1.
После необходимой очистки геометрии (defeaturing), упрощения конструктивных особенностей, не влияющих на точность расчета осколочного действия (конструкция взрывателя, оперения), в ANSYS Workbench передается трехмерная твердотельная модель БЧ изделия, представленная на рис.2.
Внешний вид конечно-элементной модели в переменных Лагранжа, выполненной в ANSYS и переданной в препроцеcсор LS-Prepost, представлен на рис.3.
Основной проблемой при моделировании было отслеживание каждого из фрагментов естественного дробления, что при общем их числе около 10000 потребовало бы дискретизации минимум на 1 млн по осколочной оболочке. С помощью построения параметризованной регулярной сетки в интересующих зонах удалось снизить общее число ячеек в расчете до 300 000, а время расчета одного варианта до 3 часов на персональной ЭВМ с процессором Intel Core i5.
По результатам моделирования были построены распределения масс и скоростей по 5-градусным угловым зонам, распределение масс представлено на рис.5. Дробление оболочки из высокоосколочной стали дает примерно в 2 раза большее число осколков, чем для сталистого чугуна, что существенно повысит эффективность действия изделия.