Матрица жёсткости

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

История и развитие понятия

Понятие матрицы жёсткости возникло в середине XX века в связи с развитием метода конечных элементов (МКЭ). Основоположниками метода считаются Р. Курант (1943), Дж. Аргирис (1954) и О. Зенкевич (1960-е). Аргирис впервые систематически применил матричную алгебру для расчёта стержневых систем, введя понятие «матрицы жёсткости элемента». В 1956 году М. Тёрнер, Р. Клаф, Г. Мартин и Л. Топп опубликовали работу, где матрица жёсткости была использована для треугольных и прямоугольных конечных элементов в плоской задаче теории упругости. С тех пор понятие стало стандартным инструментом в инженерных расчётах.

Физический смысл

Матрица жёсткости K размером n×n (где n — число степеней свободы системы) связывает вектор узловых сил F и вектор узловых перемещений U уравнением:

F = K · U

Каждый элемент kᵢⱼ матрицы — это сила, возникающая в i-й степени свободы при единичном перемещении по j-й степени свободе, при условии, что все остальные степени свободы зафиксированы. Диагональные элементы kᵢᵢ всегда положительны и представляют собой жёсткость системы по соответствующей степени свободы. Матрица жёсткости является симметричной (kᵢⱼ = kⱼᵢ), что следует из теоремы взаимности Максвелла — Бетти.

Типы матриц жёсткости

1. Матрица жёсткости конечного элемента

Относится к отдельному элементу (стержню, балке, пластине, оболочке, объёмному элементу). Формируется на основе аппроксимирующих функций формы и физических свойств материала. Для стержневого элемента, работающего на растяжение-сжатие, матрица жёсткости имеет размер 2×2: Kₑ = (EA/L) · [[1, -1], [-1, 1]], где E — модуль упругости, A — площадь поперечного сечения, L — длина элемента.

2. Глобальная матрица жёсткости

Собирается из матриц жёсткости отдельных элементов путём ансамблирования (суммирования) с учётом топологии соединений. Учитывает граничные условия (закрепления). После сборки система уравнений K·U = F решается относительно неизвестных перемещений.

3. Секционная матрица жёсткости

В строительной механике — матрица, описывающая упругие свойства поперечного сечения стержня. Для изгибаемого элемента она включает изгибную жёсткость EI, для кручения — крутильную жёсткость GIₜ.

4. Матрица жёсткости основания

В геотехнике — матрица, моделирующая упругое основание (например, по модели Винклера или Пастернака). Учитывает распределённую реакцию грунта.

Свойства матрицы жёсткости

• Симметричность: Kᵢⱼ = Kⱼᵢ (следствие принципа взаимности).
• Положительная определённость: для недеформированной системы без свободных движений (т.е. при наличии достаточных закреплений) матрица жёсткости положительно определена, что гарантирует единственность решения.
• Разрежённость: в глобальной матрице жёсткости большинство элементов равны нулю (ненулевые только для степеней свободы, связанных общими узлами). Это свойство используется для эффективного хранения и решения систем.
• Вырожденность при отсутствии закреплений: если система не закреплена (плавает), матрица жёсткости вырождена (имеет нулевые собственные значения), что соответствует возможности жёсткого смещения как целого.

Применение

Метод конечных элементов (МКЭ)

Матрица жёсткости — центральный элемент МКЭ. Используется в расчётах на прочность, устойчивость, колебания. Позволяет моделировать:

• стержневые системы (фермы, рамы);
• пластины и оболочки (корпуса, крылья);
• объёмные тела (детали машин, фундаменты);
• композитные и анизотропные материалы.

Динамический анализ

В задачах собственных колебаний решается обобщённая проблема собственных значений: (K — ω²M)·U = 0, где M — матрица масс. Собственные значения ω² дают квадраты собственных частот, а собственные векторы — формы колебаний.

Оптимизация конструкций

Матрица жёсткости используется в задачах топологической и параметрической оптимизации для поиска минимальной массы при заданных ограничениях по жёсткости.

Адаптивное перестроение сетки

В МКЭ с адаптивным измельчением сетки матрица жёсткости пересчитывается на каждом шаге.

Примеры

Пример 1: Пружина

Для пружины жёсткостью k матрица жёсткости 2×2: K = k·[[1, -1], [-1, 1]]. При закреплении левого конца (u₁ = 0) система сводится к одному уравнению: F₂ = k·u₂.

Пример 2: Ферма из двух стержней

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

Критика и ограничения

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

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

• Матрица жёсткости является аналогом тензора жёсткости в теории упругости в дискретной форме.
• В методе граничных элементов (МГЭ) аналогом является матрица влияния.
• Для систем с демпфированием вводится матрица демпфирования, которая может быть пропорциональной матрице жёсткости (модель Рэлея).
• В квантовой механике аналогом матрицы жёсткости является матрица силовых констант в молекулярной динамике.

Источники

1. Зенкевич О. К. Метод конечных элементов в технике. — М.: Мир, 1975.
2. Аргирис Дж. Современные достижения в методах расчёта конструкций с применением матриц. — М.: Стройиздат, 1968.
3. Бате К., Вилсон Е. Численные методы анализа и метод конечных элементов. — М.: Стройиздат, 1982.
4. Тимошенко С. П., Гудьер Дж. Теория упругости. — М.: Наука, 1975.
5. Розин Л. А. Метод конечных элементов для упругих систем. — М.: Стройиздат, 1971.