Расчёт ферм методом конечных элементов. Учебное пособие по дисциплине Строительная механика
.pdf
примыкает к узлу с меньшим номером. В каждой строке (а строка соответствует узлу, совпадающему с ней по номеру) значащие числа характеризуют номера элементов, соединяемых в узле, причем 1 указывает на то, что к узлу подходит начало стержня, а -1, - что к узлу примыкает его конец.
Построенная таким образом матрица, является первой матрицей инциденций ориентированного графа и определяет топологию стержневой системы.
В качестве иллюстрации приведем пример составления структурной матрицы для фермы, показанной на рис. 7.
|
1 |
0 |
0 |
0 |
1 |
0 |
0 |
|
||
|
|
|
1 |
1 |
0 |
0 |
0 |
1 |
0 |
|
|
|
C |
0 |
1 1 |
|
|
|
1 |
. |
|
S |
0 |
0 |
0 |
|||||||
|
0 |
0 |
1 1 |
0 |
0 |
0 |
|
|||
|
0 |
0 |
0 |
1 1 |
1 |
1 |
|
|||
С помощью структурной матрицы SC можно найти проекции длин стержней системы на координатные оси (рис. 8).
Y
Yj |
|
|
j |
|
n |
|
|
|
|
ny |
|
|
|
n n |
|
Yi |
|
|
|
i |
|
|
|
|
nx |
|
|
|
|
|
0 |
Xi |
X |
|
Xj |
Рис. 8. проекции длин стержней n системы на координатные оси
Для этого сформируем матрицу-столбец координат узлов фермы
10
|
|
|
|
|
|
|
|
|
|
С |
|
|
|
1 |
|
|
|
С2 |
|
|
|
|
|
С |
* |
|
|
|
|
* |
(1) |
|
|
|
|
|
|
|
|
|
|
Сm |
|
элементами которой являются матрицы-столбцы, определяющие координаты X,Y узлов фермы:
|
|
|
X i |
|
|
|
|
X j |
|
|
|
|
X m |
|
||
|
|
|
|
|
|
|||||||||||
Сi |
|
|
|
, |
С j |
|
|
|
, ……… Сm |
|
|
|
. |
|||
|
Yi |
|
Yj |
|
Ym |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Вектор проекции , элементы которого определяют проекции длины элемента на оси X,Y, можно представить как
|
|
|
|
|
|
|
|
||
|
SCT C . |
(2) |
||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
2 |
|
|
. |
|
|
. |
|
|
(3) |
|
|
|
|
|
. |
|
|
|
|
|
S |
|
Примечание. Здесь, и далее везде, индекс "Т" означает, что матрица является транспонированной.
Для n-го стержня фермы элемент матрицы-столбца:
|
|
|
|
nx |
|
|
|
cos |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
n |
|
|
|
|
|
sin |
|
n |
n |
|
n |
n |
; |
(4) |
||||
|
|
|
ny |
|
|
|
n |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
11
|
|
|
|
cos n |
|
|
1 |
|
|
|
|
|
|
|
|
||||
|
|
n |
|
sin n |
|
|
|
n |
(5) |
|
|
n |
|||||||
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|||
где |
nx и |
ny |
- |
проекции |
этого стержня в глобальной системе |
||||
координат; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
- матрица-столбец направляющих косинусов рассматриваемого |
||||||||
стержня.
Длины отдельных элементов определяются выражением:
|
|
|
n |
|
|
|
|
|
|
|
|
|
nT n . |
(6) |
|||
Сформируем матрицу направляющих косинусов для всей |
||||||||
стержневой системы |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
. |
|
(7) |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S |
|
Заметим, что матрица направляющих косинусов - это всегда квадратная матрица размера S x S.
Перейдем к рассмотрению условий равновесия отдельного n- ного стержня системы, предполагая, что сосредоточенные силы приложены только в узлах. Начало координат совместим с началом стержня; глобальную систему координат - с горизонтальной и вертикальной осями (Х,У), а местную (x,y) - с осью стержня (рис. 9).
12
Y
y 


n
i 


n
Nнn
x
j
Nкn

n
X
Y
Yнn 
i


Yкn
j
Xкn

n
n
n
Xнn X
Рис. 9. Внутренние усилия в стержнях в локальных и глобальных координатах
По длине стержня нормальная сила постоянна, т.е.:
|
Nнn NKn Nn , |
где |
Nнn , NKn - нормальные силы в начале и в конце n-го стержня, |
соответственно.
Связь между усилиями в краевых сечениях в местной ( Nнn , глобальной ( X нn ,Yнn , X Kn ,YKn ) системах координат устанавливается матричными отношениями:
|
|
|
|
Nn ; |
|
|
|
|
|
|
Nn ; |
|
|
X нn |
|
n |
|
|
|
X Kn |
n |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Nn T X нn ; |
|
|
|
Nn Tn |
X Kn ; |
(8) |
|||||||
|
|
X нn |
|
|
|
|
|
X Kn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
X нn |
|
Yнn |
, |
|
X Kn |
|
YKn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где X нn и |
|
X Kn |
- матрицы-столбцы. Они |
определяют проекции усилия |
|||||||||
Nn в начале и в конце n-го стержня, отнесенные к глобальной системе
координат, а матрица-столбец n осуществляет переход от местной системы координат к глобальной.
Рассмотрим условия равновесия узлов системы (рис. 10).
13
Y2i
X2i
2 
Y |
|
|
|
|
Yni |
|
Fyi |
Xni |
|
n |
|
|
|
|
|
|
Fxi |
i |
|
X |
|
1 |
X1i |
|
|
|
|
|
Y1i |
Рис. 10. Загружение стержней фермы примыкающих к узлу.
Внешние нагрузки, действующие в i- м узле характеризуются матрицей-столбцом:
|
|
FXi |
|
Fi |
FYi |
, |
|
|
|
||
где FXi и FYi - проекции внешних сил, |
действующих в i-м узле на |
||
координатные оси.
Уравнения равновесия i-ro узла в матричной форме имеют вид: |
|
|
|
Fi |
X n , |
где знак суммы распространяется на все примыкающие к i-му узлу стержни.
Таких уравнений нужно составить столько, сколько в системе
имеется узлов.
Совокупность векторов Fi образует матрицу-столбец внешних нагрузок
14
|
|
|
|
|
|
|
|
||
|
|
F |
|
|
|
|
1 |
|
|
|
F2 |
|
|
|
|
|
|
||
F |
... |
|
. |
|
|
|
|
|
|
|
|
Fm |
|
|
Примечание. Вектор внешних сил вестные опорные реакции.
Связь между внешними узловыми усилиями, в примыкающих к узлам матричным уравнением
|
|
|
|
F |
|
S |
C X . |
(9)
F содержит также и неиз-
нагрузками и внутренними стержнях, устанавливается
(10)
Следует заметить, что каждая строка матрицы SC объединяет стержни, сходящиеся в соответствующем узле:
|
|
|
|
|
|
|
|
|
|
|
|
X |
1 |
|
|
|
|
|
|
|
X |
2 |
|
|
|
|
|
||
X |
... |
, |
||
|
|
|
|
|
|
|
X m |
|
|
где X - матрица-столбец проекций внутренних сил в стержнях, сходящихся, соответственно, в узлах 1,2………..m.
Совместное рассмотрение выражений (8), (10) дает
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
S |
|
|
N |
S |
N ; |
|
S |
|
S |
C |
|
; |
|
|
|
|
(11) |
|
|
|
N1 |
|
|
|
|
|
|
N 2 |
|
|
|
N |
, |
(12) |
||||
... |
||||||
|
|
|
N S |
|
|
|
15
- матрица-столбец, элементами которого являются внутренние усилия в стержнях фермы.
Для определения неизвестных внутренних усилий в Ni стержнях, из вектора F следует исключить элементы, соответствующие опорным связям системы (т.е. неизвестные реакции в опорах) и сформировать вектор Q, представляющий собой вектор только активных угловых нагрузок. Эту процедуру можно осуществить, умножив обе части
уравнения (8) на прямоугольную матрицу S0 , которая получается из единичной матрицы порядка 2m путем удаления строк, соответствующих опорным связям.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
S |
0 F |
S |
|
|
S |
|
|
|
|
N |
|
||||||||||||
Итак |
|
o |
C |
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
или |
Q |
|
S |
p N ; |
(13) |
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
Q |
|
S |
o F ; |
|
|
|
|
|
(14) |
|
||||||||||||
|
|
|
|
|
|
o |
|
C |
|
|
|
|
|
|
|||||||||||
|
S p |
S |
S |
|
|
|
|
||||||||||||||||||
|
|
. |
|
|
(15) |
||||||||||||||||||||
В качестве иллюстрации приведем составление матрицы S0 для фермы, показанной на рис. 7:
|
|
|
|
0 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
|
|
|
|
|
0 |
0 |
0 |
1 |
0 |
0 |
0 |
0 |
0 |
0 |
|
|
|
|
|
0 0 0 0 1 0 0 0 0 0 |
|
|||||||||
S0 |
. |
|||||||||||||
|
|
|
|
0 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
0 |
0 |
|
|
|
|
|
0 |
0 |
0 |
0 |
0 |
0 |
1 |
0 |
0 |
0 |
|
|
|
|
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
1 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Здесь из десяти строк (2m = 10) единичной матрицы вычеркнуты две первые строки, соответствующие опорным связям по осям X и Y в узле 1, и строки, соответствующие опорным связям по осям X и Y в узле I, и строки 8 и 10 , соответствующие связям по оси Y в узлах 4 и 5.
Перейдем к рассмотрению геометрических уравнений, устанавливающих связь между деформациями и перемещениями.
16
Условия неразрывности требуют, чтобы перемещения всех концов элементов (стержней), примыкающих к i-му узлу, а также их проекции на координатные оси, были равны, т. е.:
Ui1 Ui 2 .... Ui
Vi1 Vi 2 .... Vi
где Ui1 и Vi1 - проекции перемещения соответствующего конца n -го стержня на координатные оси.
Установим связь между удлинениями отдельных элементов и
перемещениями их концов (рис. 11): |
|
n (U j Ui ) cos n (V j Vi )sin n . |
(17) |
Y

n

n i
i
Ui
|
|
j |
j |
Uj |
|
|
Vj |
|
|
|
|
|
|
n |
|
n |
|
|
|
|
Vi
0
X
Рис.11. Перемещение и удлинение стержня на плоскости
Получение зависимости (17), для большей наглядности, можно подкрепить графическим решением (рис. 12).

n




n
(Uj
n 
(Vj - Vi) - Ui)
Рис.12.Геометрическая интерпретация удлинения стержня
17
Представим выражение (17) в матричной форме:
|
|
|
|
|
|
|
|
|
|
||
n nT U ; |
|
|
|
(18) |
|||||||
|
|
|
|
|
|
|
|
U i |
|
|
|
|
|
|
|
U j |
|
||||||
|
|
||||||||||
nT |
|
cos n |
sin n |
|
; U |
V |
|
V |
|
. |
|
|
|
|
j |
||||||||
|
|
|
|
|
|
|
|
i |
|
|
|
В соответствии с выражениями типа (18), вектор стержневой системы, имеет вид:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T |
S |
CT U ; |
|
|
|
|
|
|
|
|
|
(19) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
U 2 |
|
|
|
|
|
|
|
|
|
|
... |
|
|
|
|
U i |
|
|
||
U |
|
|
; |
|
|
|
U i |
|
. |
(20) |
|
U i |
|
|
|
V |
|||||
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
i |
|
|
|
... |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U m |
|
|
|
|
|
|
|
|
|
Рассмотрим теперь физические уравнения, устанавливающие связь между удлинениями стержней n и продольными усилиями;
n Nn n /(En An ) . |
(21) |
Для всей системы стержней эта связь имеет вид:
|
|
|
|
|
|
|
|
||||
|
|
|
|
L |
n |
N ; |
(22) |
||||
|
|
|
1 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||||
|
|
|
E1 A1 |
|
|
|
|
||||
|
|
|
|
|
2 |
|
|
|
|
||
|
|
|
|
|
E2 A2 |
|
|
|
|
||
|
|
S |
|
. |
|
(23) |
|||||
L |
|||||||||||
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
|
|
|
|
|
|
ES AS |
|
18
- матрица податливости стержней фермы.
Перейдем теперь к получению зависимостей, позволяющих определять внутренние усилия в стержнях узлов системы.
Совместное рассмотрение выражений (19) и (22), дает
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
||
N |
L |
|
T |
S |
T |
U |
|
|
|||||||||
|
. |
(24) |
|||||||||||||||
|
|
|
|
|
S |
|
|
c |
|
||||||||
После подстановки формулы (24) в уравнение (13), получим: |
|||||||||||||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
||
Q |
S |
P |
|
L |
T |
S |
T U |
. |
(25) |
||||||||
|
|
|
|
S |
|
|
|
|
c |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Удалим из |
вектора |
|
U |
заведомо |
нулевые перемещения |
||||||||||||
(соответствующие опорным связям) с помощью соотношения. |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
U |
S0T V , |
|
|
|
|
|
|
|
|
(26) |
|||||||
где V - вектор содержащий только ненулевые перемещения. Подставим выражение (26) в уравнение (25) и выразим:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V |
( |
S |
P |
|
L |
S 1 |
|
S |
PT ) 1 Q |
(27) |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
или |
V |
|
L |
Q ; |
|
|
|
|
|
|
|
|
|
|
(28) |
|||||
|
|
|
|
( |
|
|
|
|
1 |
|
|
T ) 1 |
|
|||||||
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
L |
S |
P |
L |
S |
(29) |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
S |
|
|
|
P |
||||
где L - называется матрицей податливости системы.
Каждый ее элемент определяет перемещения по направлению, определяемому соответствующей строкой матрицы от единичной силы, приложенной в точке, определяемой соответствующим столбцом матрицы.
Усилия в стержнях фермы могут быть определены из совместного рассмотрения выражений (24), (26) и (15):
|
|
|
|
|
|
||
N |
L |
1 |
|
S |
T V |
. |
(30) |
|
S |
|
P |
||||
И, наконец, рассмотрим алгоритм вычисления опорных реакций шарнирно-стержневой системы. Для этого из вектора внешних
19
