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

UptoLike

ВУЗ: 

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

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

Рубрика: 

95
COMMON/MG/FMG(6,6),FMGO(6,6),MAT(20)
COMMON/XAR/IE,IT,T(3),F,ZI,NMAT,WZ,NY
DIMENSION FMGP(N,LENTA)
INTEGER T
REAL MAT
KFERM=0
KSTR=0
LLENTA=1
DO 1 I=1,N
DO 1 J=1,LENTA
1 FMGP(I,J)=0.
WRITE (8,2)
2 FORMAT(/1X,'Элемент Tип Mатрица индексов',6X,
1'F',9X,'Jz',5X,'NMAT')
DO 10 K=1,KELEM
READ(4)IE,IT,(T(J),J=1,2),NMAT,F,ZI,WZ
GO TO(101,102,102),IT
101 CALL MGSTF
KFERM=KFERM+1
GOTO 4
102 CALL MGSTR
KSTR=KSTR+1
4 DO 9 L=1,NY
II=T(L)
K1=(II-1)*NS
DO 8 M=1,NS
K1=K1+1
K2=(L-1)*NS+M
DO 7 LL=1,NY
K3=(T(LL)-1)*NS
DO 6 MM=1,NS
K4=(LL-1)*NS+MM
K5=K3+MM+1-K1
IF(K5)6,6,5
5 FMGP(K1,K5)=FMGP(K1,K5)+FMGO(K2,K4)
IF(LLENTA-K5)200,6,6
200 LLENTA=K5
6 CONTINUE
7 CONTINUE
8 CONTINUE
9 CONTINUE
10 CONTINUE
WRITE (8,12) LLENTA,KFERM,KSTR
12 FORMAT(/3X,'ПPOBEPKA',/
13X,'Лента MЖK: ',I5,/
23X,'КЭ фермы: ',I5,/
33X,'Балочных КЭ: ',I5)
RETURN
END
4.2.3. Граничные условия
Граничные условия в МКЭ автоматически учитываются вектором уз-
ловой нагрузки {P
0
}. В том случае, когда граница свободна, элементы
вектором нагрузки, соответствующие расположенным на границе узлам,
равны нулю.

Страницы