Метод конечных элементов в расчетах стержневых систем. Дьяков И.Ф - 103 стр.

UptoLike

103
2 FMGP(I,J)=FMGP(I,J)-C*FMGP(II,K)
3 CONTINUE
FMGP(II,L)=C
JJ=I+N1-1
P(JJ)=P(JJ)-C*P(N3)
4 CONTINUE
5 P(N3)=P(N3)/FMGP(II,1)
RETURN
END
SUBROUTINE SOL2(N1,N2,LENTA,FMGP,ISTR,ISTL,P,N)
DIMENSION FMGP(ISTR,ISTL),P(N)
N3=N2-N1+1
IF(N2-N)1,2,2
1 N3=N3+1
2 CONTINUE
3 N3=N3-1
IF(N3)7,7,4
4 L=N3+N1-1
DO 6 K=2,LENTA
L=L+1
IF(FMGP(N3,K))5,6,5
5 JJ=N3+N1-1
P(JJ)=P(JJ)-FMGP(N3,K)*P(L)
6 CONTINUE
GO TO 3
7 RETURN
END
4.2.5. Результаты расчёта
Решение уравнений равновесия определяет вектор {Z
0
} узловых пе-
ремещений конструкции в общей системе координат, который выводится
на печать. Вектор {Z
0
} распределяется по КЭ. Согласно выражению (1.7)
вычисляются вектора {S
r
} внутренних узловых сил КЭ в местной системе
координат и выводятся на печать.
Вычисление нормальных напряжений [13] в узлах КЭ выполняется в
соответствии с принятыми положительными направлениями узловых сил
в программах, реализующих МКЭ.
Нормальные напряжения в начале и конце КЭ фермы (началопер-
вый номер узла в матрице индексов, конецвторой номер узла) вычисля-
ются по формуле
F
N
к,н
.
Верхний знак в формуле относится к началу КЭ фермы, а нижний
знакк концу.