ВУЗ:
Составители:
93
WRITE (8,19) (MSVOB(J),J=4,6)
18 FORMAT (I7,I7,2I4)
19 FORMAT(I14,2I4)
WRITE(2)(MSVOB(J),J=1,6)
20 CONTINUE
REWIND 2
21 READ(5,31)IVID
31 FORMAT(I5)
WRITE(8,32)IVID
32 FORMAT(/5X,'Число различных сечений стержней:',I3/)
WRITE(8,33)
33 FORMAT(5X,'Геометрические характеристики сечений КЭ'/
13X,'Сечение',4X,'F',8X,'Iz',8X,'Wz')
DO 35 I=1,IVID
READ(5,34)I1,(GXS(I,J),J=1,3)
34 FORMAT(I5,3F10.2)
35 WRITE(8,36)I1,(GXS(I,J),J=1,3)
36 FORMAT(I7,3F10.2)
DO 38 I=1,KELEM
READ(5,37)IE,IT,(T(J),J=1,3),NMAT
37 FORMAT(6I5)
IV=T(3)
F=GXS(IV,1)
ZI=GXS(IV,2)
WZ=GXS(IV,3)
WRITE(4)IE,IT,(T(J),J=1,2),NMAT,F,ZI,WZ
38 CONTINUE
REWIND 4
RETURN
END
4.2.2. Формирование матрицы жёсткости конструкции
При реализации алгоритма МКЭ формирование матрицы жёсткости
конструкции непосредственно по выражению (1.4) крайне неэффективно,
так как матрица соответствий состоит из единиц и большого количества
нулей, что приводит к нерациональному использованию памяти ЭВМ.
В программном комплексе [5] не формируется матрица жёсткости
конструкции [K
0
], а формируется сразу матрица коэффициентов при неиз-
вестных системы уравнений [K
0
*
]. Формирование матрицы [K
0
*
] выполня-
ется с помощью топологических матриц: матрицы типов узлов [Т
у
] и с
помощью матрицы связи узлов [С].
Матрица типов узлов [Т
у
] – булева матрица, число строк у которой
равно числу узлов расчётной схемы конструкции, а число столбцов равно
максимальному числу степеней свободы узлов. Если расчётная схема кон-
струкции предусматривает наличие в i-м узле перемещения, соответст-
Страницы
- « первая
- ‹ предыдущая
- …
- 91
- 92
- 93
- 94
- 95
- …
- следующая ›
- последняя »