
- •Лекция 17
- •9.2. Конечный элемент. Матрица жесткости конечного элемента
- •17.3. Преобразование матрицы жесткости при переходе от одной системы координат к другой
- •17.4. Примеры построения матрицы жесткости конечного элемента
- •9.5. Формирование и решение системы уравнений мкэ. Определение внутренних усилий в элементах
- •17.6. Пример расчета стержневой системы методом конечных элементов
17.3. Преобразование матрицы жесткости при переходе от одной системы координат к другой
Обычно матрицу жесткости строят в удобной для данного элемента системе координат, которую принято называть локальной или местной. При переходе от отдельных элементов к системе элементов необходимо осуществить переход от локальных систем координат к общей для всех элементов системе координат, которую принято называть глобальной.
Пусть известны перемещения какого-либо жесткого узла k элемента в локальной системе координат, оси которой 1` и 2` повернуты на уголjотносительно осей 1 и 2 глобальной системы координат. Оси 3 и 3`, очевидно, будут совпадать, т.к. в обоих случаях рассматривается одна и та же плоскость (рис.17.15).
Пусть
и
-
перемещения узлаkпо направлению
1 в глобальной и локальной системах
координат соответственно, а
и
-
перемещения по направлению 2 в глобальной
и локальной системах координат
соответственно.
Из рис.17.15 следует:
,
.
Рис.17.15
Очевидно, угловое смещение в обоих
системах осей координат будет одинаковым,
т.е.
.
Полученные соотношения, связывающие
перемещения узла в локальной и в
глобальной системах координат, в
матричной форме будут выглядеть следующим
образом:
.
Матрица
называется матрицей направляющих
косинусов дляk-го узла. Легко
убедиться, что элемент, стоящий в ееi-ой строке иj-ом столбце равен
косинусу угла междуi-ой осью в
локальной системе координат иj-ой
осью в глобальной системе координат.
Для шарнирного элемента связь между перемещениями в локальной и глобальной системах осей будет аналогичной:
.
Для элемента, содержащего несколько узлов, связь между вектором узловых перемещений Uв глобальной системе осей координат и вектором узловых перемещенийU` в локальной системе осей координат осуществляется при помощи квазидиагональной матрицы направляющих косинусовC(e) элемента, составленной из матриц направляющих косинусов входящих в элемент узлов:
.
(17.4)
Например, для элемента, изображенного на рис.9.9, выражение (9.4) будет выглядеть следующим образом:
или
.
Легко убедиться, что это равенство соответствует уравнениям, связывающим перемещения в локальной и глобальной системах осей координат:
,
,
,
,
.
Повторив те же рассуждения для усилий,
действующих на узлы элемента, получим
аналогичную (17.4) зависимость между
векторами усилий, действующих на элемент,
построенных для глобальной системы
координат
и для локальной системы координат
:
(17.5)
Усилия, действующие на узлы элемента и узловые перемещения связаны зависимостью (9.2). Запишем ее также и для локальной системы осей координат:
.
(17.6)
Заменим в (17.6) вектора усилий и перемещений согласно зависимостям (17.4) и (17.5):
и
умножим получившееся равенство слева
на матрицу
,
т.е. на матрицу, обратную матрице
направляющих косинусов элемента:
.
Сопоставляя это выражение с (17.2), получим
зависимость, связывающую матрицы
жесткости, построенные в глобальной и
локальной системах координатных осей:
.
Известно, что у матрицы направляющих
косинусов
обратная матрица совпадает с
транспонированной матрицей. Поэтому,
окончательно получаем:
.
(17.7)
Итак, если матрица жесткости элемента построена в локальной системе осей координат, то по формуле (9.7) можно получить из нее матрицу жесткости элемента в глобальной системе осей координат.