Расчетный сервис не должен быть черным ящиком: инженер обязан знать, каким методом получены числа, в каких допущениях они верны и как их проверить. Эта статья — о том, как устроен решатель модуля балок: два вычислительных пути, метод начальных параметров для перемещений, принятая модель и ее границы, встроенные механизмы контроля точности — и сверка сквозного примера раздела с ручным счетом.
Два решателя
Модуль выбирает вычислительный путь автоматически, по кинематическому анализу схемы; текущий статус всегда виден в статус-баре («Статически определимая»).
- Аналитический решатель — для статически определимых балок без врезанных шарниров. Реакции находятся из уравнений равновесия ΣF_x = 0, ΣF_y = 0, ΣM = 0; внутренние усилия — методом сечений: балка разбивается на участки по опорам, точкам приложения сил и моментов, границам распределенных нагрузок, и на каждом участке выписываются функции Q(x), M(x), N(x) в замкнутом виде. Именно эти выкладки разворачивает пошаговый отчет.
- Конечно-элементный решатель — для статически неопределимых схем и схем с врезанными шарнирами, где уравнений равновесия недостаточно. Балка моделируется стержневыми конечными элементами с кубическими (эрмитовыми) функциями формы; узловые перемещения точны, а для корректной формы эпюр между узлами сетка автоматически сгущается — внутренние усилия при этом снимаются в произвольных сечениях элемента.
Оба пути дают один и тот же состав результатов: реакции, эпюры усилий и перемещений, проверку и подбор сечения. Различие одно: замкнутый ручной вывод формул существует только для определимых схем — у МКЭ промежуточные выкладки к компактному решению не сводятся.
Модель и допущения
Решатель реализует техническую теорию изгиба (модель Эйлера — Бернулли) в классическом наборе гипотез:
- плоский изгиб: схема, нагрузки и перемещения лежат в одной плоскости;
- гипотеза плоских сечений: сечения остаются плоскими и нормальными к изогнутой оси;
- линейная упругость материала и малые перемещения;
- деформации сдвига не учитываются (модель Тимошенко не применяется).
Практическое следствие последнего пункта: сдвиговая добавка к прогибу имеет порядок (h/L)² и для балок обычных пропорций (L/h ≥ 10) составляет доли процента — ей пренебрегают и ручные методики сопротивления материалов, которые модуль воспроизводит. Для коротких высоких балок и составных сечений с податливой стенкой недооценка прогиба растет — это граница применимости модели.
Перемещения: метод начальных параметров
Эпюры прогибов y(x) и углов поворота θ(x) аналитический решатель строит методом начальных параметров: интегрированием уравнения изогнутой оси EI·y″(x) = M(x) с универсальными уравнениями для каждого типа нагрузки. Результат — аналитические полиномы EIθ(x) и EIw(x) по участкам; начальные параметры находятся из граничных условий на опорах (нулевой прогиб на шарнирах, нулевые прогиб и угол в заделке). Экстремумы перемещений отыскиваются по самим функциям, а не по сетке — в отчет попадают точные координаты.
Из этой же механики следует символьный режим: пока жесткость не задана, полином EIw(x) полностью определен нагрузками и геометрией — модуль строит универсальную безразмерную эпюру, а физический прогиб получается делением w = EIw/EI после ввода жесткости (см. способы задания сечения).
Контроль точности
- Эталонный набор схем. Аналитический решатель и МКЭ считают один и тот же набор верификационных балок в автоматических тестах проекта; сходимость двух независимых методов — условие выпуска каждой сборки. Сквозной пример раздела взят из этого набора.
- Численная проверка эпюр. Тумблер «Проверка построения эпюр y, θ численным методом» в панели «Расчет» запускает независимое численное интегрирование уравнения изогнутой оси и сверяет его с аналитическими полиномами — контроль выполняется на ваших данных, а не только на эталонах.
- Пошаговый отчет как открытый код решения. Каждое число сводки прослеживается до уравнения с подстановкой — ошибку постановки видно построчно.
Сверка с ручным счетом
Выполним независимый расчет сквозного примера (балка L = 10 м, опоры в 2 и 8 м, сила 10 кН вверх в 1 м, трапеция 10→17 кН/м на участке 3–6 м, момент 15 кН·м в 9 м) и сопоставим с выдачей модуля — тот же формат, в котором публикуют верификационные примеры большие расчетные комплексы.
Равнодействующая трапеции R_q = (10+17)/2·3 = 40.5 кН приложена в центре тяжести эпюры нагрузки x = 3 + 1·(10+2·17)/(10+17) = 4.63 м. Из ΣM_A = 0: R_B = (−10·1 − 40.5·2.63 + 15)/(−6) ≈ 16.92 кН; из ΣM_B = 0: R_A ≈ 13.58 кН; контроль ΣF_y = 13.58 + 16.92 + 10 − 40.5 = 0.
| Величина | Ручной счет | Модуль | Расхождение |
|---|---|---|---|
| Реакция R_A, кН | 13.58 | 13.58 | 0 |
| Реакция R_B, кН | 16.92 | 16.92 | 0 |
| Q_max, кН (участок 2–3 м) | 23.58 | 23.58 | 0 |
| M_max, кН·м | 57.7 (x = 4.93 м) | 57.57 (x = 4.80 м) | 0.2% |
Реакции и поперечная сила совпадают тождественно — они выводятся из одних и тех же уравнений равновесия. Максимум момента аналитически лежит в нуле поперечной силы: 23.58 = 10t + 7t²/6 дает t = 1.93 м от начала трапеции, то есть x = 4.93 м. Значение сводки 57.57 кН·м отвечает функции M(x) в точке x = 4.80 м и отличается от аналитического экстремума на 0.2% — обе величины воспроизводятся подстановкой в функции пошагового отчета, эпюра M в окрестности максимума пологая. Для проверок прочности и подбора сечения расхождение такого порядка значения не имеет.
Теоретическая база методов — в справочнике: внутренние усилия и метод сечений, построение эпюр, прогибы и их нормирование.