Редактирование: Моделирование взрыва в SPH постановке
Внимание! Вы не авторизовались на сайте. Ваш IP-адрес будет публично видимым, если вы будете вносить любые правки. Если вы войдёте или создадите учётную запись, правки вместо этого будут связаны с вашим именем пользователя, а также у вас появятся другие преимущества.
Правка может быть отменена. Пожалуйста, просмотрите сравнение версий, чтобы убедиться, что это именно те изменения, которые вас интересуют, и нажмите «Записать страницу», чтобы изменения вступили в силу.
Текущая версия | Ваш текст | ||
Строка 10: | Строка 10: | ||
== Постановка задачи == | == Постановка задачи == | ||
− | + | ||
Дана геометрия реальной осадной пищали 17 века. Упростить геометрию задачи, провести моделирование взрыва пороха в стволе орудия, найти скорость вылета ядра. | Дана геометрия реальной осадной пищали 17 века. Упростить геометрию задачи, провести моделирование взрыва пороха в стволе орудия, найти скорость вылета ядра. | ||
Строка 38: | Строка 38: | ||
Уравнением состоянием было выбрано уравнение Джонсона-Вилкинса-Ли (Jones-Wilkins-Lee). Используется для точного описания поведения давления-объём-энергия продуктов детонации взрывчатых веществ в приложениях связанных с ускорением твердых тел. Значения действительны только для больших расходов. | Уравнением состоянием было выбрано уравнение Джонсона-Вилкинса-Ли (Jones-Wilkins-Lee). Используется для точного описания поведения давления-объём-энергия продуктов детонации взрывчатых веществ в приложениях связанных с ускорением твердых тел. Значения действительны только для больших расходов. | ||
::<math> | ::<math> | ||
− | P = A(1- | + | P = A(1-ω/(R1V))e^(-R1V) + B(1-ω/(R2V))e^(-R2V) + ωE/V |
</math> | </math> | ||
где A и B -линейные коэффициенты [Па]; R1, R2, ω - нелинейные коэффициенты [безразмерные]; P - давление [Па]; E - энергия [Па-м^3/v^3] | где A и B -линейные коэффициенты [Па]; R1, R2, ω - нелинейные коэффициенты [безразмерные]; P - давление [Па]; E - энергия [Па-м^3/v^3] | ||
− | В качестве взрывчатого материала была выбрана модификация пороха LX-04-1 с параметрами приведёнными в таблице ниже. Создаем | + | В качестве взрывчатого материала была выбрана модификация пороха LX-04-1 с параметрами приведёнными в таблице ниже. Создаем |
{| class="wikitable" | {| class="wikitable" | ||
|- | |- | ||
− | ! A !! B !! | + | ! A !! B !! R1 !! R2 !! ω !! E0 !! ρ !! υ |
|- | |- | ||
| 8,364e11 || 1,298e10 || 4,62 || 1,25 || 0,42 || 9,5e9 || 1865 || 8470 | | 8,364e11 || 1,298e10 || 4,62 || 1,25 || 0,42 || 9,5e9 || 1865 || 8470 | ||
|} | |} | ||
− | |||
− | |||
− | |||
== Решение == | == Решение == | ||
Расчёт проводится в решателе LS-DYNA, многоцелевом конечно-элементный комплексе, предназначенный для анализа высоконелинейных и быстротекущих процессов в задачах механики твердого и жидкого тела. LS-DYNA представляет возможность эффективного численного моделирования высоконелинейных термомеханических процессов. | Расчёт проводится в решателе LS-DYNA, многоцелевом конечно-элементный комплексе, предназначенный для анализа высоконелинейных и быстротекущих процессов в задачах механики твердого и жидкого тела. LS-DYNA представляет возможность эффективного численного моделирования высоконелинейных термомеханических процессов. | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
== Результаты == | == Результаты == | ||
Строка 68: | Строка 58: | ||
* В задаче происходят вылеты частиц через поверхность орудия, их процент мал по отношению к количеству частиц оставшихся в области ствола, однако это тоже влияет на решение, поскольку вылетают те частицы, у которых скорость наибольшая. Возможными вариантами решения могут являться: уменьшение размера сетки и более точная настройка контакта частицы-стенки. Дальнейшее уменьшение размера шага по времени не уменьшило количество вылетающих частиц. | * В задаче происходят вылеты частиц через поверхность орудия, их процент мал по отношению к количеству частиц оставшихся в области ствола, однако это тоже влияет на решение, поскольку вылетают те частицы, у которых скорость наибольшая. Возможными вариантами решения могут являться: уменьшение размера сетки и более точная настройка контакта частицы-стенки. Дальнейшее уменьшение размера шага по времени не уменьшило количество вылетающих частиц. | ||
* Скорее всего для данной задачи ALE постановка давала бы более точные решения при тех же затратах, поскольку контакт с границами можно было задать более точно и вылетов за пределы орудия в ALE постановке не было бы. | * Скорее всего для данной задачи ALE постановка давала бы более точные решения при тех же затратах, поскольку контакт с границами можно было задать более точно и вылетов за пределы орудия в ALE постановке не было бы. | ||
− | * Скорость ядра в момент вылета из пушки составляет порядка | + | * Скорость ядра в момент вылета из пушки составляет порядка 280 м/с, что соответствует реальным данным. |
== Список источников == | == Список источников == |