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

UptoLike

ВУЗ: 

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

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

Рубрика: 

102
KR=(I2-I1+1)/NS
CALL RDDA(NR,KR,1,NS,FMGP,ISTR,ISTL)
GO TO 6
10 IT=1
I1=1
GO TO 9
11 RETURN
END
SUBROUTINE RDDA(NN,KZ,NF,NS,FMGP,ISTR,ISTL)
DIMENSION FMGP(ISTR,ISTL)
NW=NN
N1=NF
DO 1 K=1,KZ
N2=N1+NS-1
READ(11,rec=NW)((FMGP(I,J),I=N1,N2),J=1,ISTL)
N1=N2+1
1 NW=NW+1
RETURN
END
SUBROUTINE WRDA(NN,KZ,NF,NS,FMGP,ISTR,ISTL)
DIMENSION FMGP(ISTR,ISTL)
NW=NN
N1=NF
DO 1 K=1,KZ
N2=N1+NS-1
WRITE(11,rec=NW)((FMGP(I,J),I=N1,N2),J=1,ISTL)
N1=N2+1
1 NW=NW+1
RETURN
END
SUBROUTINE MOVS(N1,N2,LENTA,FMGP,ISTR,ISTL)
DIMENSION FMGP(ISTR,ISTL)
K1=N1+1
K2=N1+N2
M=1
DO 2 I=K1,K2
DO 1 J=1,LENTA
FMGP(M,J)=FMGP(I,J)
1 CONTINUE
2 M=M+1
RETURN
END
SUBROUTINE SOL1(N1,N2,LENTA,FMGP,ISTR,ISTL,P,N)
DIMENSION FMGP(ISTR,ISTL),P(N)
II=N2-N1+1
DO 30 I=1,II
IF(FMGP(I,1))30,20,30
20 FMGP(I,1)=1.*1E7
30 CONTINUE
DO 5 N3=N1,N2
II=N3-N1+1
I=II
DO 4 L=2,LENTA
I=I+1
IF(FMGP(II,L))1,4,1
1 C=FMGP(II,L)/FMGP(II,1)
J=0
DO 3 K=L,LENTA
J=J+1
IF(FMGP(II,K))2,3,2

Страницы