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

UptoLike

ВУЗ: 

УлГТУ | Ульяновск

Составители: 

Рубрика: 

98
4 FORMAT(3X,'Узел ст. св. заданного перемещения')
DO 6 J=1,KPER
READ(5,1)NYZ,IPER
WRITE (8,5) NYZ,IPER
5 FORMAT(I6,I20)
I=IPER+(NYZ-1)*NS
FMGP(I,1)=FMGP(I,1)*1.E8
P(I)=FMGP(I,1)*P(I)
6 CONTINUE
7 RETURN
END
C РЕАЛИЗАЦИЯ CBЯЗEЙ C ЗAДAHHOЙ ЖECTKOCTЬЮ
SUBROUTINE PODAT (FMGP)
COMMON /RAZM/KYZO,NS,KELEM,N,LENTA
DIMENSION FMGP(N,LENTA)
READ (5,1) KPOD
1 FORMAT (10I5)
WRITE (8,2) KPOD
2 FORMAT (/3X,'Cвязей с заданной жёсткостью:',I5)
IF (KPOD) 8,8,3
3 WRITE (8,4)
4 FORMAT (3X,'Узел Cвязь Жёсткость')
DO 7 J=1,KPOD
READ (5,5) NYZ,NSV,GST
5 FORMAT (2I5,F10.2)
WRITE (8,6) NYZ,NSV,GST
6 FORMAT (I6,I8,F15.2)
I=(NYZ-1)*NS+NSV
FMGP(I,1)=FMGP(I,1)+GST
7 CONTINUE
8 CONTINUE
RETURN
END
4.2.4. Решение разрешающей системы уравнений
Вследствие симметрии матрицы и её ленточной структуры при реше-
нии линейной системы уравнений можно сэкономить память. Память эко-
номится благодаря тому, что только ненулевую часть матрицы требуется
хранить в ЭВМ. Эта часть ленточной матрицы располагается в памяти в
виде прямоугольного массива. Структура матрицы (рис. 4.6, а) является
следствием регулярной нумерации узлов расчётной схемы конструкции.
Половина ширины ленты матрицы, с учетом элементов главной диа-
гонали, как отмечалось, может быть вычислена по формуле:
LENTA = (R+1)·NS,
где R – наибольшая разность номеров узлов КЭ;
NS – число степеней свободы узла.

Страницы