Моделирование динамических задач МДТТ с помощью параллельной версии сеточно- характеристического метода Васюков Алексей Викторович, МФТИ, аспирант Научный руководитель: д.ф.-м.н Петров Игорь Борисович МФТИ, 2011
План Кратко о методе Примеры решения задач Прочность конструкций (многослойные преграды, композиты, трубы под давлением) Медицинские задачи (черепно-мозговая травма, экстракция катаракты) Подробно о методе и реализации Если будет время и желание МФТИ, 2011
- уравнения движения - реологические соотношения - приближение малых деформаций - матричная форма записи Математическая модель
Расщепление по направлениям. Гибридная схема Схема Куранта-Изаксона-Рис (1-ый порядок) Схема Лакса-Вендрофа (2-ой порядок) Параметр переключения Условие устойчивости Метод численного решения
Сравнение методов
Возможность расчета областей сложной геометрии. Отсутствие выделенных направлений. Неструктурированные сетки
Граничные условия Спектр матрицы A Граничные условия Свободная граница
Контактные границы Задача контактного разрыва Контактные условия Свободное скольжение Многомерный случай
Деформируемая пуля, v = 100 м/с Защитный слой Тело Пространство между пулей и защитой Защита тела МФТИ, 2011
Защита тела МФТИ, 2011
Защита тела МФТИ, 2011
Защита тела (косой удар) МФТИ, 2011
Защита тела (косой удар) МФТИ, 2011
Рассматриваемая задача Труба заполненная газом или жидкостью под давлением Образование трещины в трубе (надрез) и ее развитие под действием внутреннего давления Деформация трубы в ходе процесса Деформация и разрушение труб (1/4) МФТИ, 2011
Деформация и разрушение труб (2/4) Область надреза трубы МФТИ, 2011
Деформация и разрушение труб (3/4) Давление в трубе вызывает деформации МФТИ, 2011
Деформация и разрушение труб (4/4) Итоговое раскрытие трещины МФТИ, 2011
Задача о соударении с многослойной преградой
Задача о черепно-мозговой травме Сагитальное сечение Трансверсальное сечение Skull (elastic) Brain (viscoelastic) 10 m/s impact МФТИ, 2011
Моделирование методами МКЭ Задача о черепно-мозговой травме
Max Pressure = 100 atm Max Time = 3·10 -6 sec Моделирование сеточно-характеристическим методом Задача о черепно-мозговой травме
Зависимость повреждений от направления удара Задача о черепно-мозговой травме Максимальное сдвиговое напряжение. Сила удара одинакова во всех случаях, меняется только направление удара.
Расчет и томограмма Задача о черепно-мозговой травме
Удаление катаракты Sclera Vitreous body Muscle Crystalline Heat source in laser case or stress source in ultrasonic case Anterior chamber Surgical instruments: Laser needle (heat source) Ultrasonic needle (stress source)
Удаление катаракты Моделирование методами МКЭ
Удаление катаракты Моделирование сеточно-характеристическим методом
Удаление катаракты Области поражения
Сеточно-характеристический метод Динамические задачи Волновые процессы Аккуратный учет границ Заключение МФТИ, 2011 Существующая реализация Гибридные методы Задачи с малыми и с конечными деформациями Параллельная реализация (ускорение в 56 раз на 64 процессорах)
Спасибо за внимание! МФТИ, 2011