В этом репозитории обитают попытки честно посчитать полную упругую динамику FEniCS-ом. Пока что только в 2D.
Первая реализация примерно полностью основана на https://comet-fenics.readthedocs.io/en/latest/demo/elastodynamics/demo_elastodynamics.py.html
Результат первой реализации выглядит примерно так:
Это короткий импульс справа области, прямоугольный по времени.
По меткому выражению @amisto текущий результат выглядит как корги из хлебушка. На первый взгляд мило, а потом возникают вопросы. (Цитата неточная.)
Примерный перечень вопросов:
- Несимметрично, волны направо и налево разные идут.
- Продольная волна, которая большая, выглядит не как одна волна, а как вермишель из нескольких фронтов, как будто переотражение где-то было.
- Неясно, где сдвиговая волна - то ли склеилась с волнами Релея, то ли пропала. Может быть, это вопрос подбора параметров материала.
- Отражения красивые, но в одном месте странная волна, которую надо проверять - потенциальный численный глюк.