Liren Large Software Subsidy 9

home *** CD-ROM | disk | FTP | other *** search

/ Liren Large Software Subsidy 9 / 09.iso / e / e003 / 10.ddi / PASST.FOR < prev next >

Wrap

Text File | 1987-06-23 | 66.4 KB | 830 lines

SUBROUTINE PASS 00161290 IMPLICIT REAL*8(A-H,O-Z) 00161300 REAL*8 NPAR 00161310 COMMON /MISC/ NBLOCK,NEQB,LL,NF,LB 00161320 COMMON /JUNK/ DUK(200),KKK(4),NDYN,NRJUNK(49) R0161330 COMMON /BAND/ KOPT,NRBAND(7) R0161340 DIMENSION T(3) 00161350 COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00161360 & ,RRELPA(24) R0161361 COMMON /OUT/ IDUMM(4),IOSIG,IODISP,NROUT(4) R0161370 COMMON /HEADIN/ TITLE1(8),TITLE2(2),TITLE3(4),RRHEAD(21) R0161380 COMMON A(1) 00161390 COMMON /TRASH/ TIM(400),RRTRAS(90) R0161400 COMMON /PREP/ XD(2),KSKIP,RRPREP(8) R0161410 LL=NF 00161420 NEQB=NEQ 00161430 LB=1 00161440 100 FORMAT(8F10.0) 00161450 WRITE(6, 110) LL 00161460 110 FORMAT (1X ,20X,38HTHERE WILL BE OUTPUT AT THE FOLLOWING ,I5,7H TI00161470 1MES.//) 00161480 IF (LL.LE.400) GO TO 130 00161490 WRITE(6, 120) 00161500 120 FORMAT (1H0,20X,34HTOO MANY TIMES HAVE BEEN SELECTED.//) 00161510 KSKIP=1 00161520 130 READ(5, 100) (TIM(I),I=1,LL) 00161530 TL=-1.0 00161540 DO 140 I=1,LL 00161550 IF (TIM(I).LT.0.0) GO TO 150 00161560 IF (TIM(I).LE.TL) GO TO 150 00161570 140 TL=TIM(I) 00161580 GO TO 170 00161590 150 WRITE(6, 160) 00161600 160 FORMAT (/20X,49HERROR - TIMES MUST BE POS. AND IN ASCENDING ORDER/00161610 1) 00161620 KSKIP=1 00161630 170 WRITE(6, 180) (TIM(I),I=1,LL) 00161640 180 FORMAT (20X,6G15.6) 00161650 IF(KSKIP.EQ.1) RETURN 00161660 NT14=14 00161670 ND=1 00161680 NV=N1+NEQ 00161690 REWIND NT14 00161700 NREC=2*((NEQ-1)/5+1) 00161710 NT=2 00161720 REWIND NT 00161730 LY=0 00161740 TH=-1.0 00161750 TL=-10.0 00161760 TIME=-10.0 00161770 LFL=0 00161780 DO 260 I=1,LL 00161790 IF (I.EQ.1) GO TO 190 00161800 IF (TIM(I).LE.TIM(LY)) GO TO 260 00161810 IF((TIM(I)-TIM(LY)).LT.(TH-TL)) TIM(I)=TH*(2.0+1.0E-8)-TL 00161820 190 CONTINUE 00161830 IF(LFL.EQ.1) GO TO 200 00161840 TL=TIME 00161850 CALL RDDIS (TIME,1,NT14,KFL,0,I) 00161860 IF (KFL.EQ.1) GO TO 270 00161870 TH=TIME 00161880 200 CALL RDDIS (A(N1),NEQ,NT14,KFL,0,I) 00161890 IF (KFL.EQ.1) GO TO 320 00161900 CALL RDDIS (A(NV),NEQ,NT14,KFL,0,NREC) 00161910 CALL RDDIS (A(NV),NEQ,NT14,KFL,0,NREC) 00161920 IF (KFL.EQ.1) GO TO 320 00161930 TL=TIME 00161940 CALL RDDIS (TIME,1,NT14,KFL,0,I) 00161950 IF (KFL.EQ.1) GO TO 270 00161960 TH=TIME 00161970 IF (I.GT.1) GO TO 210 00161980 IF (TIM(I).LT.TL) GO TO 230 00161990 210 DIFF1=TIM(I)-TL 00162000 DIFF2=TIME-TIM(I) 00162010 IF (TIM(I).GE.TL.AND.TIM(I).LE.TH) GO TO 220 00162020 GO TO 200 00162030 220 IF (DIFF1.GT.DIFF2) GO TO 240 00162040 230 IF(LY.EQ.0) GO TO 231 00162050 IF(TIM(LY).EQ.TL) GO TO 260 00162060 231 CONTINUE 00162070 LY=LY+1 00162080 TIM(LY)=TL 00162090 LFL=1 00162100 GO TO 250 00162110 240 JFL=0 00162120 IF(LY.EQ.0) GO TO 241 00162130 IF(TIM(LY).EQ.TH) JFL=1 00162140 IF(JFL.EQ.1) GO TO 242 00162150 241 LY=LY+1 00162160 TIM(LY)=TH 00162170 242 CONTINUE 00162180 LFL=0 00162190 CALL RDDIS (A(N1),NEQ,NT14,KFL,0,I) 00162200 IF (KFL.EQ.1) GO TO 320 00162210 CALL RDDIS (A(NV),NEQ,NT14,KFL,0,NREC) 00162220 CALL RDDIS (A(NV),NEQ,NT14,KFL,0,NREC) 00162230 IF (KFL.EQ.1) GO TO 320 00162240 IF(JFL.EQ.1) GO TO 260 00162250 250 CALL RDWRT (NT,A(N1),NEQ,13,I) 00162260 260 CONTINUE 00162270 GO TO 280 00162280 270 CONTINUE 00162290 CALL RDWRT (NT,A(N1),NEQ,13,I) 00162300 LY=LY+1 00162310 TIM(LY)=TL 00162320 280 REWIND NT 00162330 IF(LY.EQ.1) GO TO 281 00162340 IF(TIM(LY).EQ.TIM(LY-1)) LY=LY-1 00162350 281 CONTINUE 00162360 LL=LY 00162370 NBLOCK=LL 00162380 CALL SECOND (T(1)) 00162390 N2=N1+NUMNP*3 00162400 N3=N2+6 00162410 IF (IODISP.EQ.1) TITHOL=TITLE3(3) 00162420 CALL DPRNTD (A(N1),A(N2),A(N3),NUMNP,LL,NEQ,NT) 00162430 IF (IODISP.EQ.1) REWIND 6 00162440 IF (IODISP.EQ.1) TITLE3(3)=TITHOL 00162450 CALL SECOND (T(2)) 00162460 NADD=1 00162470 IF (KOPT.GT.0) NADD=NUMNP 00162480 N2=N1 00162490 N2A=N2+NEQ 00162500 N3=N2A+NADD 00162510 NDYN=0 00162520 KL=1 00162530 IF (IOSIG.EQ.1) TITHOL=TITLE3(3) 00162540 CALL DSTRES (A(N3),LL,NEQ,A(KL),A(N2A),NADD) 00162550 IF (IOSIG.EQ.1) REWIND 6 00162560 IF (IODISP.EQ.1) WRITE (6,300) 00162570 IF (IOSIG.EQ.1) WRITE (6,310) 00162580 IF (IOSIG.EQ.1) TITLE3(3)=TITHOL 00162590 CALL SECOND (T(3)) 00162600 T(1)=T(2)-T(1) 00162610 T(2)=T(3)-T(2) 00162620 WRITE (6,290) T(1),T(2) 00162630 RETURN 00162640 290 FORMAT (27H1....TIME LOG (CPU SECONDS)///33H PRINT DISPLACEMENTS.00162650 1.......... ,F8.2//33H COMPUTE STRESSES.............. ,F8.2//) 00162660 300 FORMAT (///20X,34(1H*)/20X,34HDISPLACEMENTS WILL NOT BE PRINTED./200162670 10X,34(1H*)//) 00162680 310 FORMAT (///20X,29(1H*)/20X,29HSTRESSES WILL NOT BE PRINTED./20X,2900162690 1(1H*)//) 00162700 320 WRITE(6, 330) NT14 00162710 330 FORMAT (///20X,35HAN ERROR WAS FOUND ON RESTART TAPE ,I4///) 00162720 KSKIP=1 00162730 RETURN 00162740 END 00162750 SUBROUTINE RDDIS (A,N,NT,K,J,NREC) 00196520 IMPLICIT REAL*8(A-H,O-Z) 00196530 DIMENSION A(N) 00196540 K=0 00196550 IF (J.GT.0) GO TO 40 00196560 READ (NT,10,END=20) A 00196570 GO TO 30 00196580 10 FORMAT (5E16.9) 00196590 20 K=1 00196600 30 RETURN 00196610 40 DO 50 I=1,NREC 00196620 READ (NT,10,END=20)X 00196630 50 CONTINUE 00196640 RETURN 00196650 END 00196660 SUBROUTINE DPRNTD (ID,D,B,NUMNP,LL,NEQ,NT) 00062880 IMPLICIT REAL*8(A-H,O-Z) 00062890 REAL*8 ID 00062900 COMMON /PREP/ XMX,XAD,KSKIP,NDYN,I1,RRPREP(7) R0062910 INTEGER DUMMY 00062920 COMMON /OUT/ NRES,NSTR,NDIS,DUMMY,IOSIG,IODISP,NROUT(4) R0062930 COMMON/ELPAR/XPAR(14),NDUM(8),MTOT,NDUM1(6),NUMEL,NUMEL2,NELPA(41)R0062940 COMMON A(1) 00062950 DIMENSION ID(NUMNP,3), B(NEQ), D(6) 00062960 IF (NDIS.LT.0) RETURN 00062970 REWIND NT 00062980 REWIND 8 00062990 READ (8) ID 00063000 REWIND 17 00063010 REWIND 18 00063020 NDPBLK=(MTOT-6-NUMNP)/6 00063030 NBLK=(NUMNP-1)/NDPBLK+1 00063040 KK=1 00063050 IF (NDIS.GT.0) WRITE (NDIS,20) KK,LL,NDYN,NDIS,NSTR,NUMNP,NUMEL,NU00063060 1MEL2 00063070 DO 30 LD=1,LL 00063080 READ (NT) B 00063090 DO 30 N=1,NUMNP 00063100 DO 10 I=1,6 00063110 D(I)=0.0 00063120 CALL UNPKID (ID,NUMNP,W,WX,2,N,I) 00063130 NNN=W 00063140 IF (NNN.LT.1) GO TO 10 00063150 D(I)=B(NNN) 00063160 10 CONTINUE 00063170 20 FORMAT (8I5) 00063180 WRITE (18) D 00063190 30 CONTINUE 00063200 K=1 00063210 N2=K+NUMNP 00063220 N3=N2+6 00063230 N4=N3+6*NDPBLK 00063240 CALL WRDISD (A(K),A(N2),A(N3),NUMNP,LL,NDPBLK,NDIS,NBLK) 00063250 RETURN 00063260 END 00063270 SUBROUTINE WRDISD (NORD,A,B,NUMNP,LL,NDPBLK,NDIS,NBLK) 00327500 REAL*8 A,B 00327510 DIMENSION NORD(NUMNP), A(6), B(NDPBLK,6) 00327520 COMMON /TRASH/ TIM(400),RRTRAS(580) R0327530 COMMON /OUT/ KDUMMY(9),KROT 00327540 COMMON /BAND/ KOPT,NRBAND(7) R0327550 NT1=17 00327560 NT2=18 00327570 REWIND NT2 00327580 IF (KOPT.GT.0) READ (NT1) NORD 00327590 DO 70 M=1,LL 00327600 WRITE(6, 10) TIM(M) 00327610 10 FORMAT (1X ,/20X,24HDISPLACEMENTS AT TIME = ,G15.6//) 00327620 WRITE (6,110) 00327630 KSHF2=0 00327640 KSHF=1-NDPBLK 00327650 KNT=0 00327660 DO 60 I=1,NBLK 00327670 KOUNT=0 00327680 KK=NUMNP+1 00327690 KSHF=KSHF+NDPBLK 00327700 KSHF2=KSHF2+NDPBLK 00327710 DO 30 KK=1,NUMNP 00327720 READ (NT2) A 00327730 KCH=KK 00327740 IF (KOPT.GT.0) KCH=NORD(KK) 00327750 IF (KCH.LT.KSHF.OR.KCH.GT.KSHF2) GO TO 30 00327760 KNT=KNT+1 00327770 KOUNT=KOUNT+1 00327780 NSHFT=KCH-KSHF+1 00327790 DO 20 K=1,6 00327800 20 B(NSHFT,K)=A(K) 00327810 IF (KOUNT.EQ.NDPBLK.OR.KNT.EQ.NUMNP) GO TO 40 00327820 30 CONTINUE 00327830 40 KNT1=(I-1)*NDPBLK 00327840 DO 50 J=1,NDPBLK 00327850 KNT1=KNT1+1 00327860 IF (KNT1.GT.NUMNP) GO TO 70 00327870 WRITE (6,80) KNT1,(B(J,K),K=1,6) 00327880 IF (NDIS.GT.0) WRITE (NDIS,90) KNT1,(B(J,K),K=1,3) 00327890 IF(NDIS.GT.0.AND.KROT.EQ.2) WRITE(NDIS,90) KNT1,(B(J,K),K=4,6) 00327900 WRITE(33,200)KNT1 00327910 WRITE(33,201)(B(J,K),K=1,3) 00327920 200 FORMAT(I5) 00327930 201 FORMAT(3(F20.10)) 00327940 50 CONTINUE 00327950 60 CONTINUE 00327960 KBL=0 00327970 WRITE(33,200)KBL 00327980 70 IF (KOPT.GT.0) WRITE (6,100) 00327990 80 FORMAT (1H0,I4,1P3E12.3,3E11.2) 00328000 90 FORMAT (I10,7E10.4/(8E10.4)) 00328010 100 FORMAT (//3X,46H*** NOTE *** NODE NUMBERS ARE ORIGINAL NUMBERS//) 00328020 110 FORMAT (40H0.......NODE DISPLACEMENTS AND ROTATIONS//5H NODE,7X, 00328030 11HX,11X,1HY,11X,1HZ,9X,2HXX,9X,2HYY,9X,2HZZ) 00328040 RETURN 00328050 END 00328060 SUBROUTINE DSTRES (D,LL,NEQ,SA,NORD,NADD) 00064960 IMPLICIT REAL*8(A-H,O-Z) 00064970 REAL*8 NPAR 00064980 DIMENSION D(NEQ) 00064990 DIMENSION SA(1) 00065000 COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,MEQ00065010 1,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN 00065020 INTEGER DUMMY 00065030 COMMON /TRASH/ TIM(400),RRTRAS(90) R0065040 COMMON /OUT/ NRES,NSTR,NDIS,DUMMY,IOSIG,IODISP,NROUT(4) R0065050 DIMENSION NORD(NADD) 00065060 COMMON /JUNK/ SIG(200),MM,L,KDU,NTAG,NDYN,NRJUNK(49) R0065070 COMMON /GPS/ NEQ4(10),NRGPS(10) R0065080 IF (NSTR.LT.0) RETURN 00065090 NT1=1 00065100 NT3=3 00065110 NT9=17 00065120 IF (NADD.GT.1) REWIND NT9 00065130 IF (NADD.GT.1) READ (NT9) NORD 00065140 NT=LL 00065150 LB=1 00065160 LH=0 00065170 CALL RDWRT (NT3,SA,1,6,J) 00065180 NELAST=0 00065190 N=2 00065200 IF (NSTR.GT.0) WRITE (NSTR,10) N,LL,NDYN,NDIS,NSTR,NUMNP,NUMEL,NUM00065210 1EL2 00065220 10 FORMAT (8I5) 00065230 NTU=2 00065240 REWIND NTU 00065250 DO 260 II=1,NT 00065260 WRITE(6, 20) TIM(II) 00065270 20 FORMAT (1X ,/20X,19HSTRESSES AT TIME = ,G15.6//) 00065280 LT=LH+1 00065290 LLT=1-LT 00065300 LH=LT+LB-1 00065310 IF (LH.GT.LL) LH=LL 00065320 IF (NSTR.GT.0) WRITE (NSTR,10) LT,LH 00065330 READ (NTU) D 00065340 LK=LH-LT+1 00065350 CALL RDWRT (NT1,SA,1,6,J) 00065360 NPAR(1)=0 00065370 NUME=NUMEL+NUMEL2 00065380 DO 260 MM=1,NUME 00065390 CALL RDWRT (NT1,SA,NEMN,0,KOUNT) 00065400 Z=SA(KOUNT) 00065410 IF (Z.NE.7.OR.MM.GT.NUMEL) GO TO 60 00065420 NELAST=7 00065430 WRITE (6,30) MM 00065440 30 FORMAT (/20X,7HELEMENT,I5,19H IS A BLANK ELEMENT/) 00065450 IF (NDYN.EQ.3) CALL RDWRT (NT3,Z,1,1,J) 00065460 DO 40 L=LT,LH 00065470 40 IF (NSTR.GT.0) WRITE (NSTR,50) L 00065480 50 FORMAT (3X,1H1,I2,2X,2H 7,6G10.4) 00065490 GO TO 260 00065500 60 CONTINUE 00065510 NS1=SA(KOUNT-1) 00065520 ND1=SA(KOUNT-2) 00065530 K=0 00065540 KSTR=ND1*NS1+ND1 00065550 DO 250 L=LT,LH 00065560 K=K+1 00065570 LMIN1=(L-1)*NS1+KSTR 00065580 DO 70 N=1,NS1 00065590 NPN=N+ND1 00065600 SIG(N)=0.0 00065610 DO 70 J=1,ND1 00065620 NELM=NPN+(J-1)*NS1 00065630 JJ=SA(J) 00065640 IF (JJ.LE.0) GO TO 70 00065650 SIG(N)=SIG(N)+SA(NELM)*D(JJ) 00065660 70 CONTINUE 00065670 MTYPE=Z 00065680 NTAG=1 00065690 IF (MTYPE.NE.NELAST) NTAG=0 00065700 NELAST=MTYPE 00065710 GO TO (80,90,100,110,120,130,140,150,160,210,170,170,170,23000065720 1 ), MTYPE 00065730 80 CALL DTRUSS 00065740 GO TO 240 00065750 90 CONTINUE 00065760 SIG(200)=NS1 00065770 CALL DBEAM 00065780 GO TO 240 00065790 100 CALL DPLANE 00065800 GO TO 240 00065810 110 CALL DAXIS 00065820 GO TO 240 00065830 120 NELM=KSTR 00065840 SIG(13)=SA(NELM+4) 00065850 SIG(14)=SA(NELM+6) 00065860 SIG(15)=NS1 00065870 CALL DTHRED 00065880 GO TO 240 00065890 130 CALL DSHELL 00065900 GO TO 240 00065910 140 SIG(150)=DABS(SA(KSTR+1)) 00065920 SIG(151)=DABS(SA(KSTR+2)) 00065930 CALL DBOUND(NORD,NADD) 00065940 GO TO 240 00065950 150 CALL DPLANE 00065960 GO TO 240 00065970 160 CONTINUE 00065980 SIG(200)=NS1 00065990 CALL DELBST 00066000 GO TO 240 00066010 170 NELM=KSTR-1 00066020 JJ=(NS1-16)/4 00066030 IF (JJ.LE.0) JJ=1 00066040 KK=NS1 00066050 DO 180 N=1,JJ 00066060 NELM=NELM+4 00066070 KK=KK+1 00066080 180 SIG(KK)=SA(NELM) 00066090 SIG(150)=NS1 00066100 IF (NS1.LE.4) GO TO 200 00066110 JJ=JJ+1 00066120 JF=NS1/4 00066130 DO 190 N=JJ,JF 00066140 NELM=NELM+4 00066150 KK=KK+1 00066160 SIG(KK)=SA(NELM) 00066170 KK=KK+1 00066180 190 SIG(KK)=SA(NELM+1) 00066190 200 SIG(151)=Z 00066200 CALL ST2DD (NORD,NADD) 00066210 GO TO 240 00066220 210 NELM=KSTR-2 00066230 JJ=NS1/6 00066240 KK=NS1 00066250 DO 220 N=1,JJ 00066260 NELM=NELM+6 00066270 KK=KK+1 00066280 220 SIG(KK)=SA(NELM) 00066290 SIG(150)=NS1 00066300 CALL ST3DD (NORD,NADD) 00066310 GO TO 240 00066320 230 CALL SIXSTD 00066330 240 CONTINUE 00066340 250 CONTINUE 00066350 260 CONTINUE 00066360 RETURN 00066370 END 00066380 SUBROUTINE DTRUSS 00066920 IMPLICIT REAL*8 (A-H,O-Z) 00066930 REAL*8 NPAR 00066940 COMMON A(1) 00066950 COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00066960 1,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN 00066970 COMMON /JUNK/ SIG(12),EXTRA(188),MM,L,K,NTAG,NDYN,NRJUNK(49) R0066980 COMMON /OUT/ NRES,NSTR,NDIS,NROUT(7) R0066990 10 IF (NTAG.EQ.0) WRITE (6,30) 00067000 WRITE (6,40) MM,L,SIG(1),SIG(2) 00067010 NTAG=1 00067020 IF (NSTR.GT.0) WRITE (NSTR,20) L,SIG(1),SIG(2) 00067030 20 FORMAT (3X,1H2,I2,2X,2H 1,6G10.4) 00067040 RETURN 00067050 30 FORMAT (23H1 TRUSS MEMBER ACTIONS//46H0 MEMBER LOAD ST00067060 1RESS FORCE/) 00067070 40 FORMAT (2I8,F15.0,F15.0) 00067080 END 00067090 SUBROUTINE DBEAM 00052470 IMPLICIT REAL*8 (A-H,O-Z) 00052480 COMMON /JUNK/ SIG(26),EXTRA(174),MM,L,K,NTAG,NDYN,NRJUNK(49) R0052490 COMMON /OUT/ NRES,NSTR,NDIS,NBMSTR,NROUT(6) R0052500 10 IF (NTAG.EQ.0) WRITE (6,30) 00052510 IF (NTAG.EQ.0.AND.NBMSTR.EQ.1) WRITE (6,40) 00052520 IF (NTAG.EQ.0) WRITE (6,50) 00052530 NS=EXTRA(174) 00052540 IF (NS.GT.12) WRITE (6,60) MM,L,(SIG(I1),I1=1,6),(SIG(I2),I2=13,1900052550 1),(SIG(I3),I3=7,12),(SIG(I4),I4=20,26) 00052560 IF (NS.LE.12) WRITE (6,70) MM,L,(SIG(I1),I1=1,12) 00052570 NTAG=1 00052580 IF (NSTR.GT.0) WRITE (NSTR,20) NS,L,(SIG(I),I=1,NS) 00052590 20 FORMAT (I4,I2,2X,2H 2,7G10.4/(8G10.4)) 00052600 RETURN 00052610 30 FORMAT (/59H1....BEAM FORCES(LBS), MOMENTS(IN-LBS), AND STRESSES(P00052620 1SI) //1X,7HELEMENT,2X,4HLOAD,2X,7HSTATION,15X,5HAXIAL,2(6X,5HSHEA00052630 2R),4X,9HTORSION ,2(5X,13HB E N D I N G,4X)/2X,6HNUMBER,2X,4HCASE,00052640 311X,5HFORCE,11X,2HR1,9X,2HR2,9X,2HR3,9X,2HM1,20X,2HM2,20X,2HM3) 00052650 40 FORMAT (25X,6HSTRESS,7X,5HR1/A1,3X,8HAT Q3/B3,4X,8HAT Q2/B2,5X,4H-00052660 1 - ,1X,2X,9H AT +C3,2X,9H AT -C3,2X,9H AT +C2,2X,9H AT -C200052670 2) 00052680 50 FORMAT (1X) 00052690 60 FORMAT (4X,I4,2X,I4,4X,5HEND-I,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X00052700 1,F11.0/23X,2X,6HSTRESS,1X,3F11.0,11X,4F11.0/18X,5HEND-J,2X,6HFORCE00052710 2 ,1X,4F11.0,11X,F11.0,11X,F11.0/23X,2X,6HSTRESS,1X,3F11.0,11X,4F1100052720 3.0) 00052730 70 FORMAT (4X,I4,2X,I4,4X,5HEND-I,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X00052740 1,F11.0/18X,5HEND-J,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X,F11.0) 00052750 END 00052760 SUBROUTINE DPLANE 00062090 IMPLICIT REAL*8 (A-H,O-Z) 00062100 COMMON /OUT/ NRES,NSTR,NDIS,NROUT(7) R0062110 COMMON /JUNK/ SIG(12),EXTRA(188),MM,L,K,NTAG,NDYN,NRJUNK(49) R0062120 10 IF (NTAG.EQ.0) WRITE (6,40) 00062130 CC=(SIG(1)+SIG(2))/2.0 00062140 BB=(SIG(1)-SIG(2))/2. 00062150 CR=DSQRT(BB**2+SIG(4)**2) 00062160 SIG(5)=CC+CR 00062170 SIG(6)=CC-CR 00062180 SIG(7)=0.0 00062190 EF=(SIG(3)-SIG(5))**2+(SIG(5)-SIG(6))**2+(SIG(6)-SIG(3))**2 00062200 EF=DSQRT(EF/2.) 00062210 IF ((BB.EQ.0.0).AND.(SIG(4).EQ.0.0)) GO TO 20 00062220 SIG(7)=28.648D0*DATAN2(SIG(4),BB) 00062230 20 WRITE (6,50) MM,L,(SIG(I),I=1,7),EF 00062240 NTAG=1 00062250 IF (NSTR.GT.0) WRITE (NSTR,30) L,SIG(5),SIG(6),EF,(SIG(I),I=1,4) 00062260 30 FORMAT (3X,1H7,I2,4H12 3,7G10.4) 00062270 RETURN 00062280 40 FORMAT (1X ,37H PLANE STRESS/STRAIN ELEMENT STRESSES// 00062290 112H0EL.NO. LOAD,7X,8HR-STRESS,7X, 00062300 18HZ-STRESS,7X,8HT-STRESS,6X,9HRZ-STRESS,5X,10HMAX-STRESS,5X, 00062310 210HMIN-STRESS,3X,5HANGLE,5X,6HSIG-EF/) 00062320 50 FORMAT (2I6,6F15.0,F8.3,F10.0) 00062330 END 00062340 SUBROUTINE DAXIS 00052180 IMPLICIT REAL*8 (A-H,O-Z) 00052190 REAL*8 NPAR 00052200 COMMON /JUNK/ SIG(200),MM,L,K,NTAG,NDYN,NRJUNK(49) R0052210 COMMON /ELPAR/ NPAR(14),NUMNP,MBAND,NELTYP,N1,N2,N3,N4,N5,MTOT,NEQ00052220 1,N2P,N3P,NMRI,NTRI,N1P,NUMEL,NUMEL2,KZ(20,2),NEMN 00052230 COMMON /OUT/ NRES,NSTR,NDIS,NROUT(7) R0052240 10 IF (NTAG.EQ.0) WRITE (6,40) 00052250 CC=(SIG(1)+SIG(2))/2.0 00052260 BB=(SIG(1)-SIG(2))/2. 00052270 CR=DSQRT(BB**2+SIG(4)**2) 00052280 SIG(5)=CC+CR 00052290 SIG(6)=CC-CR 00052300 SIG(7)=0.0 00052310 EF=(SIG(3)-SIG(5))**2+(SIG(5)-SIG(6))**2+(SIG(6)-SIG(3))**2 00052320 EF=DSQRT(EF/2.) 00052330 IF ((BB.EQ.0.0).AND.(SIG(4).EQ.0.0)) GO TO 20 00052340 SIG(7)=28.648D0*DATAN2(SIG(4),BB) 00052350 20 WRITE (6,50) MM,L,(SIG(I),I=1,7),EF 00052360 IF (NSTR.GT.0) WRITE (NSTR,30) L,SIG(5),SIG(6),EF,(SIG(I),I=1,4) 00052370 30 FORMAT (3X,1H7,I2,4H12 4,7G10.4) 00052380 NTAG=1 00052390 RETURN 00052400 40 FORMAT (1X ,30H AXISYMMETRIC ELEMENT STRESSES// 00052410 112H0EL.NO. LOAD,7X,8HR-STRESS, 00052420 17X,8HZ-STRESS,7X,8HT-STRESS,6X,9HRZ-STRESS,5X,10HMAX-STRESS,5X, 00052430 210HMIN-STRESS,3X,5HANGLE,5X,6HSIG-EF/) 00052440 50 FORMAT(2I6,6F15.0,F8.3,F10.0) 00052450 END 00052460 SUBROUTINE DTHRED 00066390 IMPLICIT REAL*8(A-H,O-Z) 00066400 REAL*8 NS 00066410 REAL*8 IS1,IS2,IS3,IS4 00066420 DIMENSION SPR(6) 00066430 DIMENSION EFS(2) 00066440 COMMON /JUNK/ SIG(15),EXRA(185),MM,L,K,NTAG,NDYN,NRJUNK(49) R0066450 COMMON /OUT/ NRES,NSTR,NDIS,NROUT(7) R0066460 EQUIVALENCE (IS3,SIG(13)), (IS4,SIG(14)), (NS,SIG(15)) 00066470 10 CONTINUE 00066480 IF (NTAG.EQ.0) WRITE (6,50) 00066490 NG=NS 00066500 IS1=IS3 00066510 IS2=IS4 00066520 IF (IS1.LT.1) IS1=IS1*100.+0.001 00066530 IF (IS2.LT.1) IS2=IS2*100.+0.001 00066540 NF1=IS1 00066550 NF2=IS2 00066560 IS1=NF1 00066570 IS2=NF2 00066580 IF (NF1.EQ.0) NF1=7 00066590 IF (NF2.EQ.0) NF2=7 00066600 NF=10*NF1+NF2 00066610 IF (NF.EQ.77) NF=71 00066620 EFS(2)=0.0 00066630 SPR(4)=0.0 00066640 SPR(5)=0.0 00066650 CALL DPRIST (IS1,IS2,SIG,SPR,NS) 00066660 II=NG/6 00066670 DO 30 J=1,II 00066680 M=6*J-6 00066690 EFS(J)=(SIG(M+1)-SIG(M+2))**2+(SIG(M+2)-SIG(M+3))**2+(SIG(M+3)-00066700 1 SIG(M+1))**2 00066710 DO 20 I=4,6 00066720 MI=M+I 00066730 20 EFS(J)=EFS(J)+6.*SIG(MI)**2 00066740 30 EFS(J)=DSQRT(EFS(J)/2.) 00066750 WRITE (6,60) MM,L,IS1,(SIG(I),I=1,6),(SPR(I),I=1,3),EFS(1) 00066760 IF (NG.EQ.12) WRITE (6,70) IS2,(SIG(I),I=7,12),(SPR(I),I=4,6),EFS(00066770 12) 00066780 II=6+NG 00066790 IF (NSTR.GT.0) WRITE (NSTR,40) II,L,NF,SPR(1),SPR(2),EFS(1),SPR(4)00066800 1,SPR(5),EFS(2),(SIG(I),I=1,NG) 00066810 40 FORMAT (I4,I2,I2,2H 5,7G10.4/(8G10.4)) 00066820 RETURN 00066830 50 FORMAT (1X ,35H .....8-NODE SOLID ELEMENT STRESSES// 00066840 113H ELEM. LOAD F, 00066850 162HACE SIG-XX SIG-YY SIG-ZZ SIG-XY SIG-YZ SIG-ZX,00066860 238H SIG-MAX SIG-MIN S2 OR SIG-EF/16H NO. NO. NO.,84X,00066870 35HANGLE/) 00066880 60 FORMAT (I5,I4,F5.0,1X,10F10.0) 00066890 70 FORMAT (9X,F5.0,1X,10F10.0) 00066900 END 00066910 SUBROUTINE DPRIST (IS1,IS2,SIG,SPR,NS) 00062350 IMPLICIT REAL*8 (A-H,O-Z) 00062360 REAL*8 IS1,IS2,NS 00062370 DIMENSION SIG(12), SPR(6), IS(2), SG(6) 00062380 IS(1)=IS1 00062390 IS(2)=IS2 00062400 NNS=1 00062410 IF (NS.EQ.12.) NNS=2 00062420 DO 60 N=1,NNS 00062430 IN=3*N-3 00062440 II=IN*2 00062450 IF (IS(N).EQ.0) GO TO 10 00062460 CC=(SIG(II+1)+SIG(II+2))/2. 00062470 BB=(SIG(II+1)-SIG(II+2))/2. 00062480 CR=DSQRT(BB**2+SIG(II+4)**2) 00062490 SPR(IN+1)=CC+CR 00062500 SPR(IN+2)=CC-CR 00062510 SPR(IN+3)=0. 00062520 IF (BB.NE.0.) SPR(IN+3)=28.648*DATAN2(SIG(II+4),BB) 00062530 GO TO 60 00062540 10 CC=(SIG(II+1)+SIG(II+2)+SIG(II+3))/3. 00062550 DO 20 I=1,3 00062560 SG(I)=SIG(II+I)-CC 00062570 20 SG(I+3)=SIG(II+I+3) 00062580 C2=(SG(1)**2+SG(2)**2+SG(3)**2)*.5+SG(4)**2+SG(5)**2+SG(6)**2 00062590 C3=SG(1)*(SG(2)*SG(3)-SG(5)*SG(5))+SG(4)*(SG(5)*SG(6)-SG(4)*SG(00062600 1 3))+SG(6)*(SG(4)*SG(5)-SG(2)*SG(6)) 00062610 IF (C2.EQ.0.0) C2=1.0E-08 00062620 T=DSQRT(C2/1.5) 00062630 A=C3*1.414214/T**3 00062640 AXQ=A+1.0 00062650 IF (AXQ.GT.0.0.AND.AXQ.LT.2.0) A=DACOS(A)/3.0 00062660 IF (AXQ.GE.2.0) A=0.0 00062670 IF (AXQ.LE.0.0) A=1.0471976 00062680 T=T*1.414214 00062690 SPR(IN+1)=T*DCOS(A) 00062700 SPR(IN+2)=T*DCOS(A+2.0944) 00062710 SPR(IN+3)=T*DCOS(A-2.0944) 00062720 DO 30 I=2,3 00062730 IF (SPR(IN+1).GT.SPR(IN+I)) GO TO 30 00062740 C3=SPR(IN+1) 00062750 SPR(IN+1)=SPR(IN+I) 00062760 SPR(IN+I)=C3 00062770 30 CONTINUE 00062780 IF (SPR(IN+2).LE.SPR(IN+3)) GO TO 40 00062790 C3=SPR(IN+2) 00062800 SPR(IN+2)=SPR(IN+3) 00062810 SPR(IN+3)=C3 00062820 40 DO 50 I=1,3 00062830 50 SPR(IN+I)=SPR(IN+I)+CC 00062840 60 CONTINUE 00062850 RETURN 00062860 END 00062870 SUBROUTINE DBOUND(NORD,NADD) 00052770 IMPLICIT REAL*8 (A-H,O-Z) 00052780 COMMON /JUNK/ SIG(12),EXRA(188),MM,L,K,NTAG,NDYN,NRJUNK(49) R0052790 DIMENSION NORD(NADD) 00052800 10 IF (NTAG.EQ.0) WRITE (6,20) 00052810 NN=EXRA(138)*10000.+.001 00052820 IF(EXRA(138).GT.1.0) NN=EXRA(139)*10000.+.001 00052830 IF(NADD.GT.1.AND.NN.LE.NADD.AND.NN.GT.0) NN=NORD(NN) 00052840 IF(NN.LE.0) NN=0 00052850 WRITE(6,30) MM,NN,L,SIG(1),SIG(2) 00052860 RETURN 00052870 20 FORMAT(1X ,39H CONSTRAINT FORCE - BOUNDARY ELEMENTS // 00052880 $77H EL. NUMBER NODE N LOAD CASE FORCE(LBS.) M00052890 $OMENT(IN-LBS.) /) 00052900 30 FORMAT(1X,I10,I12,I11,4X,2F15.0) 00052910 END 00052920 SUBROUTINE DELBST 00055930 IMPLICIT REAL*8(A-H,O-Z) 00055940 COMMON /JUNK/ SIG(39),EXR(161),MM,L,K,NTAG,RRJUNK(25) R0055950 COMMON /OUT/ N,NSTR,NDIS,NBMSTR,NROUT(6) R0055960 IF (NTAG.EQ.0) WRITE (6,20) 00055970 IF (NTAG.EQ.0.AND.NBMSTR.EQ.1) WRITE (6,30) 00055980 IF (NTAG.EQ.0) WRITE (6,10) 00055990 10 FORMAT (1X) 00056000 20 FORMAT (85H1E L B O W - F O R C E S(LBS), M O M E N T S(IN-LBS), A00056010 1 N D S T R E S S E S(PSI) //1X,7HELEMENT,2X,4HLOAD,2X,7HSTATIO00056020 2N,15X,5HAXIAL,2(6X,5HSHEAR),4X,9HTORSION ,2(5X,13HB E N D I N G,400056030 3X)/2X,6HNUMBER,2X,4HCASE,11X,5HFORCE,11X,2HRX,9X,2HRY,9X,2HRZ,9X,200056040 4HMX,20X,2HMY,20X,2HMZ) 00056050 30 FORMAT (25X,6HSTRESS,7X,5HRX/A1,3X,8HAT Q3/B3,4X,8HAT Q2/B2,5X,4H-00056060 1 - ,1X,2X,9H AT +C3,2X,9H AT -C3,2X,9H AT +C2,2X,9H AT -C200056070 2) 00056080 NS=EXR(161) 00056090 IF (NS.GT.18) WRITE (6,50) MM,L,(SIG(I1),I1=1,6),(SIG(I2),I2=19,2500056100 1),(SIG(I3),I3=7,12),(SIG(I4),I4=26,32),(SIG(I5),I5=13,18),(SIG(I6)00056110 2,I6=33,39) 00056120 IF (NS.LE.18) WRITE (6,60) MM,L,(SIG(I1),I1=1,18) 00056130 NTAG=1 00056140 IF (NSTR.GT.0) WRITE (NSTR,40) NS,L,(SIG(I),I=1,NS) 00056150 40 FORMAT (I4,I2,2X,2H 9,7G10.4/(8G10.4)) 00056160 50 FORMAT (4X,I4,2X,I4,4X,5HEND-I,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X00056170 1,F11.0/23X,2X,6HSTRESS,1X,3F11.0,11X,4F11.0/17X,6HCENTER,2X,6HFORC00056180 2E ,1X,4F11.0,11X,F11.0,11X,F11.0/23X,2X,6HSTRESS,1X,3F11.0,11X,4F100056190 31.0/18X,5HEND-J,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X,F11.0/23X,2X,600056200 4HSTRESS,1X,3F11.0,11X,4F11.0) 00056210 60 FORMAT (4X,I4,2X,I4,4X,5HEND-I,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X00056220 1,F11.0/17X,6HCENTER,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X,F11.0/18X,00056230 25HEND-J,2X,6HFORCE ,1X,4F11.0,11X,F11.0,11X,F11.0) 00056240 RETURN 00056250 END 00056260 SUBROUTINE ST2DD (NORD,NADD) 00262790 IMPLICIT REAL*8(A-H,O-Z) 00262800 COMMON /JUNK/ SIG(200),MM,L,K,NTAG,RRJUNK(25) R0262810 COMMON /BAND/ KOPT,NRBAND(7) R0262820 COMMON /OUT/ NRES,NSTR,NROUT(8) R0262830 COMMON /PREP/ XMX,XAD,NRPREP(17) R0262840 DIMENSION NORD(NADD) 00262850 COMMON /QTSARG/ SIGEX(100),X(2),RRQTSA(898) R0262860 DIMENSION C(4), D(4), Q(3) 00262870 DATA Q/4HNODE,4HC.G.,4HI.P./ 00262880 JJ=SIG(150)/4 00262890 MTYP=SIG(151) 00262900 KL=0 00262910 IF (NTAG.GT.0) GO TO 70 00262920 IF (MTYP.EQ.11.AND.NTAG.EQ.0) WRITE (6,10) 00262930 IF (MTYP.EQ.12.AND.NTAG.EQ.0) WRITE (6,20) 00262940 IF (MTYP.EQ.13.AND.NTAG.EQ.0) WRITE (6,30) 00262950 10 FORMAT (1H0,32H .... AXISYMMETRIC ELEMENTS ..../) 00262960 20 FORMAT (1H0,32H .... PLANE STRESS ELEMENTS ..../) 00262970 30 FORMAT (1H0,32H .... PLANE STRAIN ELEMENTS ..../) 00262980 IF (NTAG.EQ.0.AND.MTYP.NE.11) WRITE (6,40) 00262990 IF (NTAG.EQ.0.AND.MTYP.EQ.11) WRITE (6,50) 00263000 40 FORMAT (12H0EL.NO. LOAD,12X, 00263010 1 90HY-STRESS Z-STRESS YZ-STRESS T-ST00263020 1RESS MAX-STRESS MIN-STRESS ANGLE SIG-EF LOCATION/) 00263030 50 FORMAT (12H0EL.NO. LOAD,12X, 00263040 1 90HR-STRESS Z-STRESS RZ-STRESS T-ST00263050 1RESS MAX-STRESS MIN-STRESS ANGLE SIG-EF LOCATION/) 00263060 IF (NTAG.EQ.0.AND.KOPT.GT.1.AND.NADD.EQ.1) WRITE (6,60) 00263070 60 FORMAT (5X,39HNOTE - NODES ARE NEW RENUMBERED NODES /) 00263080 70 KK=SIG(150) 00263090 JL=1000 00263100 IF (JJ.GT.1) JL=(KK-16)/4 00263110 DO 140 I=1,JJ 00263120 KK=KK+1 00263130 FACE=SIG(KK)*10000.+.001 00263140 NF=FACE 00263150 LL=NF 00263160 A=Q(2) 00263170 IF (I.EQ.1) GO TO 80 00263180 IF (NADD.GT.1.AND.I.LE.JL) NF=NORD(LL) 00263190 A=Q(1) 00263200 IF (I.EQ.1) A=Q(2) 00263210 IF (I.LE.JL) GO TO 80 00263220 A=Q(3) 00263230 NF=I-JL 00263240 LL=SIG(KK) 00263250 IF (SIG(KK).LT.0.0) LL=LL-1 00263260 X(1)=(SIG(KK)-LL-XAD)*XMX 00263270 KK=KK+1 00263280 LL=SIG(KK) 00263290 IF (SIG(KK).LT.0.0) LL=LL-1 00263300 X(2)=(SIG(KK)-LL-XAD)*XMX 00263310 80 CONTINUE 00263320 LL=4*(I-1) 00263330 DO 90 J=1,4 00263340 LL=LL+1 00263350 90 D(J)=SIG(LL) 00263360 CALL MXMND (D,C(1),C(2),C(3),C(4)) 00263370 DO 100 K=1,4 00263380 KL=KL+1 00263390 100 SIGEX(KL)=C(K) 00263400 IF (I.EQ.1) WRITE (6,120) MM,L,D,C,A 00263410 IF (I.GT.1.AND.I.LE.JL) WRITE (6,130) D,C,A,NF 00263420 IF (I.GT.JL) SIGEX(KL-1)=X(1) 00263430 IF (I.GT.JL) SIGEX(KL)=X(2) 00263440 IF (I.GT.JL) GO TO 140 00263450 IF (I.GT.JL.AND.L.EQ.1) WRITE (6,110) X(1),X(2),D,C,A,NF 00263460 IF (I.GT.JL.AND.L.GT.1) WRITE (6,130) D,C,A,NF 00263470 110 FORMAT (3H Y=,F8.2,3H Z=,F8.2,F10.0,5F11.0,F7.2,F9.0,2X,A4,I5) 00263480 120 FORMAT (2I6,9X,6F11.0,F7.2,F9.0,2X,A4,I5) 00263490 130 FORMAT (21X,6F11.0,F7.2,F9.0,2X,A4,I5) 00263500 140 CONTINUE 00263510 KK=SIG(150) 00263520 LL=JJ*4 00263530 NS=KK+LL 00263540 IF (NSTR.GT.0) WRITE (NSTR,150) NS,L,MTYP,(SIGEX(I),I=1,LL),(SIG(I00263550 1),I=1,KK) 00263560 150 FORMAT (I4,I2,2X,I2,7G10.4/(8G10.4)) 00263570 RETURN 00263580 END 00263590 SUBROUTINE MXMND (STRESS,P1,P2,AG,EF) 00152110 IMPLICIT REAL*8 (A-H,O-Z) 00152120 DIMENSION STRESS(1) 00152130 CC=(STRESS(1)+STRESS(2))*0.5E0 00152140 BB=(STRESS(1)-STRESS(2))*0.5E0 00152150 CR=DSQRT(BB**2+STRESS(3)**2) 00152160 P1=CC+CR 00152170 P2=CC-CR 00152180 AG=45.0E0 00152190 IF (DABS(BB).LT.1.0D-8) GO TO 10 00152200 AG=28.648D0*DATAN2(STRESS(3),BB) 00152210 10 EF=(P1-P2)**2+(P2-STRESS(4))**2+(STRESS(4)-P1)**2 00152220 EF=DSQRT(EF/2.0) 00152230 RETURN 00152240 END 00152250 SUBROUTINE ST3DD (NORD,NADD) 00264960 IMPLICIT REAL*8(A-H,O-Z) 00264970 COMMON /JUNK/ SIG(200),MM,L,K,NTAG,RRJUNK(25) R0264980 COMMON /BAND/ KOPT,NRBAND(7) R0264990 COMMON /OUT/ NR,NSTR,NROUT(8) R0265000 COMMON /QTSARG/ SIGEX(90),RRQTSA(910) R0265010 DIMENSION C(6), D(12), Q(2) 00265020 DIMENSION NORD(NADD) 00265030 DATA Q/4HNODE,4HC.G./ 00265040 JJ=SIG(150)/6. 00265050 IF (NTAG.GT.0) GO TO 40 00265060 IF (NTAG.EQ.0) WRITE (6,10) 00265070 10 FORMAT (1H0,39H .... 3-D SOLID ELEMENTS 8-20NODE .... /) 00265080 IF (NTAG.EQ.0.AND.KOPT.GT.1.AND.NADD.EQ.1) WRITE (6,20) 00265090 20 FORMAT (5X,35H NODE NUMBERS ARE RENUMBERED NODES /) 00265100 IF (NTAG.EQ.0) WRITE (6,30) 00265110 30 FORMAT (118H ELEM. LOAD LOC. SIG-XX SIG-YY SIG-ZZ SIG-00265120 1XY SIG-YZ SIG-ZX SIG-MAX SIG-MIN SIG-MID SIG-EF/)00265130 40 KK=SIG(150) 00265140 A=0. 00265150 B=0. 00265160 E=6. 00265170 KL=0 00265180 DO 100 I=1,JJ 00265190 KK=KK+1 00265200 FACE=SIG(KK)*10000.+.001 00265210 NF=FACE 00265220 LL=NF 00265230 IF (NADD.GT.1.AND.I.GT.1) NF=NORD(LL) 00265240 LL=6*(I-1) 00265250 DO 50 J=1,6 00265260 LL=LL+1 00265270 50 D(J)=SIG(LL) 00265280 CALL DPRIST (A,B,D,C,E) 00265290 EFS=(C(1)-C(2))**2+(C(2)-C(3))**2+(C(3)-C(1))**2 00265300 EFS=DSQRT(EFS/2.) 00265310 DO 60 K=1,3 00265320 KL=KL+1 00265330 60 SIGEX(KL)=C(K) 00265340 KL=KL+1 00265350 SIGEX(KL)=EFS 00265360 IF (I.GT.1) GO TO 80 00265370 WRITE (6,70) MM,L,Q(2),(D(J),J=1,6),(C(J),J=1,3),EFS 00265380 70 FORMAT (I5,I4,1X,A4,1X,10F10.0) 00265390 GO TO 100 00265400 80 WRITE (6,90) Q(1),NF,(D(J),J=1,6),(C(J),J=1,3),EFS 00265410 90 FORMAT (5X,A4,I5,1X,10F10.0) 00265420 100 CONTINUE 00265430 KK=SIG(150) 00265440 LL=JJ*4 00265450 NS=KK+LL 00265460 IF (NSTR.GT.0) WRITE (NSTR,110) NS,L,(SIGEX(I),I=1,LL),(SIG(I),I=100265470 1,KK) 00265480 110 FORMAT (I4,I2,2X,2H10,7G10.4/(8G10.4)) 00265490 RETURN 00265500 END 00265510 SUBROUTINE SIXSTD 00238470 IMPLICIT REAL*8(A-H,O-Z) 00238480 COMMON /OUT/ NRES,NSTR,NROUT(8) R0238490 COMMON /JUNK/ SIG(6),EXTRA(194),MM,L,K,NTAG,RRJUNK(25) R0238500 IF (NTAG.EQ.0) WRITE (6,20) 00238510 WRITE (6,10) MM,L,SIG 00238520 NTAG=1 00238530 RETURN 00238540 10 FORMAT (2I6,6F15.1) 00238550 20 FORMAT (//10X,48H SIX BY SIX STIFFNESS ELEMENT FORCES AND MOMENTS/00238560 120X,18H(LOCAL DIRECTIONS)/ 00238570 112H0EL.NO. LOAD,5X,10HFORCE - XX,5X,10HFORCE - YY, 00238580 25X,10HFORCE - ZZ,4X,11HMOMENT - XX,4X,11HMOMENT - YY,4X, 00238590 311HMOMENT - ZZ) 00238600 END 00238610