Метод Ньюмарка¶

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

¶История и предпосылки создания

В середине XX века развитие строительной механики и сейсмостойкого строительства потребовало точных методов расчёта несущих конструкций на динамические воздействия. Классические аналитические решения были применимы лишь для простейших систем с одной степенью свободы. Для реальных сооружений с множеством степеней свободы требовались численные подходы. Ньюмарк, работавший в Иллинойсском университете в Урбана-Шампейне, разработал метод, который позволял эффективно интегрировать уравнения движения, сохраняя устойчивость при больших временных шагах. Первоначально метод был ориентирован на расчёт железобетонных конструкций при землетрясениях, но впоследствии нашёл применение в машиностроении, авиастроении и других отраслях.

¶Математическая основа

Метод Ньюмарка основан на аппроксимации перемещений и скоростей на временном шаге с помощью разложения в ряд Тейлора с учётом ускорений. Уравнение движения системы с несколькими степенями свободы в матричной форме имеет вид:

\[ M \ddot{u}(t) + C \dot{u}(t) + K u(t) = F(t) \]

где \(M\), \(C\), \(K\) — матрицы масс, демпфирования и жёсткости соответственно; \(\ddot{u}\), \(\dot{u}\), \(u\) — векторы ускорений, скоростей и перемещений; \(F(t)\) — вектор внешних сил.

Основные предположения метода на шаге интегрирования \([t_n, t_{n+1}]\):

• Ускорение изменяется по линейному закону (или по параболическому — в зависимости от выбора параметров).
• Скорость и перемещение на следующем шаге выражаются через ускорения на текущем и следующем шагах.

Формулы Ньюмарка:

\[ \dot{u}_{n+1} = \dot{u}_n + (1 - \gamma) \Delta t \ddot{u}_n + \gamma \Delta t \ddot{u}_{n+1} \]

\[ u_{n+1} = u_n + \Delta t \dot{u}_n + \left( \frac{1}{2} - \beta \right) \Delta t^2 \ddot{u}_n + \beta \Delta t^2 \ddot{u}_{n+1} \]

где \(\Delta t\) — шаг по времени; \(\gamma\) и \(\beta\) — параметры метода, определяющие его устойчивость и точность.

¶Параметры метода

Выбор параметров \(\gamma\) и \(\beta\) критически влияет на свойства алгоритма:

• Безусловная устойчивость (метод не требует малого шага по времени для устойчивости) достигается при \(\gamma \geq 0.5\) и \(\beta \geq \frac{1}{4}(\gamma + 0.5)^2\).
• Условная устойчивость — при других значениях.
• Точность — при \(\gamma = 0.5\) метод сохраняет энергию системы (отсутствует численное демпфирование), при \(\gamma > 0.5\) вносится искусственное демпфирование.

Наиболее распространённые варианты:

¶Алгоритм реализации

На практике метод Ньюмарка реализуется как неявная схема: на каждом шаге необходимо решать систему линейных алгебраических уравнений относительно \(\ddot{u}_{n+1}\). Алгоритм для шага \(n \to n+1\):

1. Вычислить эффективную матрицу жёсткости:

\[ \hat{K} = K + \frac{\gamma}{\beta \Delta t} C + \frac{1}{\beta \Delta t^2} M \]

1. Вычислить эффективный вектор нагрузки:

\[ \hat{F}_{n+1} = F_{n+1} + M \left( \frac{1}{\beta \Delta t^2} u_n + \frac{1}{\beta \Delta t} \dot{u}_n + \left( \frac{1}{2\beta} - 1 \right) \ddot{u}_n \right) + C \left( \frac{\gamma}{\beta \Delta t} u_n + \left( \frac{\gamma}{\beta} - 1 \right) \dot{u}_n + \Delta t \left( \frac{\gamma}{2\beta} - 1 \right) \ddot{u}_n \right) \]

1. Решить систему \(\hat{K} \ddot{u}_{n+1} = \hat{F}_{n+1}\).
2. Обновить скорости и перемещения по формулам Ньюмарка.

Для больших систем (тысячи и миллионы степеней свободы) решение системы на каждом шаге требует значительных вычислительных ресурсов, поэтому часто применяются итерационные решатели или методы понижения размерности.

¶Применение

Метод Ньюмарка является стандартным инструментом в программных комплексах конечно-элементного анализа (FEA), таких как ANSYS, ABAQUS, LS-DYNA, NASTRAN, а также в специализированных пакетах для сейсмического расчёта (SAP2000, ETABS, SCAD). Основные области применения:

• Сейсмостойкое строительство — расчёт зданий, мостов, плотин на землетрясения. Метод позволяет моделировать нелинейное поведение материалов (пластичность, трещинообразование) при циклических нагрузках.
• Машиностроение — анализ ударных нагрузок, вибраций роторов, динамики механизмов.
• Авиационная и ракетно-космическая техника — расчёт конструкций на акустические и ударные воздействия.
• Геотехника — моделирование динамики грунтов и взаимодействия с фундаментами.

¶Достоинства и недостатки

¶Достоинства

• Безусловная устойчивость при правильном выборе параметров — позволяет использовать относительно большие шаги по времени.
• Высокая точность для низкочастотных составляющих отклика.
• Возможность введения искусственного демпфирования для подавления высокочастотных шумов.
• Хорошо изучен и реализован во всех основных FEA-пакетах.

¶Недостатки

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

¶Сравнение с другими методами

Метод Ньюмарка относится к семейству методов прямого интегрирования. Его альтернативы:

• Явные методы (например, метод центральных разностей) — не требуют решения системы на шаге, но условно устойчивы (шаг должен быть меньше критического, определяемого наименьшим периодом системы). Эффективны для задач с короткими временными интервалами и высокочастотными процессами (удары, взрывы).
• Методы Ньюмарка с модификациями (например, метод Хилбера — Хьюза — Тейлора, HHT) — позволяют улучшить демпфирование высоких частот при сохранении точности.
• Методы спектрального разложения (модальный анализ) — эффективны для линейных задач с малым числом мод, но не применимы к нелинейным системам.

¶Интересные факты

• Натан Ньюмарк (1910–1981) был одним из пионеров сейсмостойкого строительства. Его метод до сих пор остаётся одним из самых популярных в инженерной практике.
• Параметр \(\beta\) в методе Ньюмарка часто называют «параметром Ньюмарка», а \(\gamma\) — «параметром демпфирования».
• В 1970-х годах метод был расширен на задачи с нелинейным демпфированием и пластичностью, что позволило моделировать разрушение конструкций при сильных землетрясениях.

¶Источники

• Newmark, N. M. (1959). A Method of Computation for Structural Dynamics. Journal of the Engineering Mechanics Division, ASCE, 85(EM3), 67–94.
• Bathe, K. J. (2006). Finite Element Procedures. Prentice Hall.
• Chopra, A. K. (2012). Dynamics of Structures: Theory and Applications to Earthquake Engineering. Prentice Hall.
• Clough, R. W., Penzien, J. (1993). Dynamics of Structures. McGraw-Hill.
• Zienkiewicz, O. C., Taylor, R. L. (2005). The Finite Element Method for Solid and Structural Mechanics. Elsevier.