№6 Сетка сошлась, грунт разошелся

Я ранее уже рассказывала, что «Каменика» пришла в геотехнику из машиностроения. И когда мы начали знакомиться с практикой расчетов в этой сфере, то наткнулись на вещь, которая нас, прямо скажем, изрядно удивила.

Конечно-элементная сетка

В машиностроении исследование сеточной сходимости — это святое. Прям база базовое. Первое, чему учат молодого расчетчика: не можешь доказать, что результат не зависит от размера сетки — не имеешь права им пользоваться. Это как в лаборатории: не откалибровал прибор — не пиши протокол. И второе, чему учат: всегда делай локальное сгущение в зонах ожидаемых высоких (концентраторов) напряжений — у отверстий, входящих углов, мест приложения нагрузок.

Придя в геотехнику, мы увидели, что многие обходятся без этого. Причем часто по объективным причинам: геологические модели огромны, расчет может идти днями, а в коммерческой работе за каждый час приходится платить. В итоге сетка строится «на глазок», и если картинка в PLAXIS или MIDAS выглядит правдоподобно — расчет считается завершенным.

Проблема в том, что это не просто «погрешность». Это слепое пятно. В одной из статей про типичные ошибки проектирования армогрунтовых стен с использованием МКЭ был сделан убийственный вывод: многие инженеры учатся «нажимать кнопки» по туториалу, не понимая физики и математики метода. А старый принцип GIGO (Garbage In — Garbage Out) никто не отменял.

О чем говорит метод

Теорема сходимости МКЭ обещает нам точное решение, но только при бесконечном измельчении сетки, что физически невозможно (хотя и здесь всё не так гладко). Слишком мелкая сетка — это экспоненциальный рост времени счета. Слишком крупная — решение уравнения может не «попасть в конечный элемент». Именно поэтому исследование сходимости — являет собой задачу про поиск баланса между потребной точностью и вычислительными затратами.

Простой пример

Представьте, что вы считаете, ну пусть будет прочность двутавровой балки, которая испытвает одновременно изгиб, кручение и распределенную нагрузку. На очень крупной сетке конструкция может показать ложный  двухкратный запас проности. Начинаете измельчать — и только на третьем шаге видите, что коэффициент запаса падает до 1.1. 

Еще более критична ситуация с концентраторами напряжений — отверстиями, углами, локальными облегчениями. Грубая сетка их просто не видит. Это полная потеря информации. Локальное сгущение сетки в таких зонах — необходимость, позволяющая уловить градиенты напряжений.

Что в сухом остатке

1. Без исследования сходимости ваш расчет — это просто цветная картинка. Вы не знаете, насколько она врет. Может, на 2%, а может, на 40% в опасную сторону.

2. Локальное сгущение — это стандарт. Если его нет, вы просто не видите реальных градиентов напряжений в опасной зоне.

Нас, как производителей приборов, это задевает напрямую. Поэтому при проектировании новых едениц изделий, всегда обращаем внимания на методику оценки предварительной работоспособности конструкции.