1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
|
npoutre = nelem - nconnec - nbbar
lele=0
je=1
IF ( npoutre.GT.0 ) THEN
1100 CONTINUE
READ(15,*) je,(node(i,je),i=1,npele),
& itype(je),loibet(je),etype(je)
IF (itype(je).GT.ntype) THEN
PRINT*,'ERREUR dans la section de lelement ',je
STOP
ENDIF
IF (etype(je).EQ.1 .AND. loibet(je).GT.nbloibet) THEN
PRINT*,'ERREUR dans la loi de lelement ',je
STOP
ENDIF
incele=0
IF (incele .LE. 0) THEN
incele=1
ENDIF
1200 CONTINUE
lele= lele+1
IF (je-lele) 1500,1500,1300
1300 CONTINUE
itype(lele)=itype(lele-1)
loibet(lele)=loibet(lele-1)
etype(lele)=etype(lele-1)
DO 10 i=1,npele
node(i,lele)=node(i,lele-1)+incele
10 CONTINUE
1500 CONTINUE
IF (je-lele) 1600,1600,1200
1600 CONTINUE
IF (npoutre-je) 1700,2000,1100
1700 CONTINUE
STOP
2000 CONTINUE
nodmax=0
DO 20 ii=1,npoutre
DO 20 jj=1,npele
DO 20 ll=1,npele
noddiff=IABS(node(jj,ii)-node(ll,ii))
IF (noddiff .GT. nodmax) nodmax=noddiff
20 CONTINUE
ENDIF |
Partager