o o o 470 12 2 20 30 C- C C C C C C C C- 40 WRITE(NOT,4003) E(2),NU(2),THICK(2) WRITE(*,4004) E(3),NU(3) WRITE(NOT,4004) E(3),NU(3) ELSE DO 12 1=1,M WRITE(*,5003) I,E(I),NU(I),THICK(I) WRITE(NOT,5003) I,E(I),NU(I),THICK(I) WRITE(*,5004) NLAYS,E(NLAYS),NU(NLAYS) WRITE(NOT,5004) NLAYS,E(NLAYS),NU(NLAYS) ENDIF READ NUMBER OF LOADS AND THEIR PARAMETERS BI NLOAD = 2 NZEP = .FALSE. BI NZEQ = .FALSE. BI IDENT = IREF1 DO 30 1=1,NLOAD BI LDSTRS(I) = 500.0 RADIUS(I) =4.5 HOSTR(I) =0.0 PSI(I) = 0.0 XCI) = 0.0 PSI(I)=.0174533*PS1(1) BI IF(LDSTRS(I).GT.ACCUR(1)) NZEP= .TRUE. BI IF(HOSTRd) .GT.ACCUR(l)) NZEQ = .TRUE. BI IF(IDENT.EQ.IREF1) GO TO 20 BI IF(IDENT.NE.IREF2) WRITE(NOUT,9040) LDSTRS(I),HOSTR(I) BI GO TO 30 BI LDSTRS(I) = LDSTRS(I)/(3.14159*RADIUS(I)*RADIUS(I)) BI HOSTR(I) = HOSTR(I)/(3.14159*RADIUS(I)*RADIUS(I)) BI CONTINUE BI Y(1) = -10.0 Y(2) = 10.0 B1 TEST ON OBVIOUS MISTAKES IN SYSTEMS DATA-BI CARDS. BI WHEN IRED > 0 THE REDUCED SPRINGCOMPLIAN- BI CE WAS READ. BI A NON-VANISHING SLIPRESISTANCE IS SUBSTI- BI TUTED TO PREVENT RIGID-BODY MOTION OF THE BI TOPLAYERS BI ----BI DO 50 J = 1,NLAYS BI IF((1.0-NU(J)).LT.ACCUR(1)) GO TO 410 BI IF(E(J).LT.ACCUR(1)) GO TO 420 BI IF(J.EQ.NLAYS) GO TO 50 BI IF(IRED.EQ.0) GO TO 40 BI ALK(J) = AK(J) BI IF(ALK(J).LT.1000.0.OR..NOT.NZEQ) GO TO 50 BI ALK(J) = 1000.0 BI AK(J) = 1000.0 BI GO TO 50 BI ALK(J) = AK(J)*E(J)/(1.0+NU(J)) BI IF(ALK(J).LT.1000.0.OR..NOT.NZEQ) GO TO 50 BI ALK(J) = 1000.0 BI