avasyukov / fenics-elastic-waves

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

В этом репозитории обитают попытки честно посчитать полную упругую динамику FEniCS-ом. Пока что только в 2D.

Первая реализация примерно полностью основана на https://comet-fenics.readthedocs.io/en/latest/demo/elastodynamics/demo_elastodynamics.py.html

Результат первой реализации выглядит примерно так: pulse2d demo

Это короткий импульс справа области, прямоугольный по времени.

По меткому выражению @amisto текущий результат выглядит как корги из хлебушка. На первый взгляд мило, а потом возникают вопросы. (Цитата неточная.)

Примерный перечень вопросов:

  1. Несимметрично, волны направо и налево разные идут.
  2. Продольная волна, которая большая, выглядит не как одна волна, а как вермишель из нескольких фронтов, как будто переотражение где-то было.
  3. Неясно, где сдвиговая волна - то ли склеилась с волнами Релея, то ли пропала. Может быть, это вопрос подбора параметров материала.
  4. Отражения красивые, но в одном месте странная волна, которую надо проверять - потенциальный численный глюк.

About


Languages

Language:Python 100.0%