| C*** EVALUATE SHEAR STRESS RESULTANTS60 IF(IFSHE.EQ.0) RETURN | STRP | 27 |
| STRP | 28 |
| ESHXX=DGRAD(1)-XZROT | STRP | 29 |
| ESHYY=DGRAD(4)-YZROT | STRP | 30 |
| STRES(4)=DSHER(1,1)*ESHXX | STRP | 31 |
| STRES(5)=DSHER(2,2)*ESHYY | STRP | 32 |
| RETURN | STRP | 33 |
| END | STRP | 34 |
# 9.5.15 Subroutine SUBMP
This subroutine evaluates $[B_{i}]^{T}D[B_{j}]detJ\times Gauss$ weights and is used in the evaluation of the element stiffness matrices.
```csv
SUBROUTINE SUBMP (BIMAT,BJMAT,DAREA,DMATX,ESTIF,INODE, SUBP 1
. JNODE,NCOLI,NROIJ,NCOLJ) SUBP 2
C**************************SUBP 3
C SUBP 4
C*** CARRY OUT MATRIX MULTIPLICATION SUBP 5
C SUBP 6
C**************************SUBP 7
DIMENSION BIMAT(NROIJ,NCOLI),BJMAT(NROIJ,NCOLJ), SUBP 8
. DMATX(NROIJ,NROIJ),DBMAT(3,3), SUBP 9
. ESTIF(27,27),SBSTF(3,3) SUBP 10
C*** EVALUATE D*BJ SUBP 11
DO 10 J=1,NCOLJ SUBP 12
DO 10 I=1,NROIJ SUBP 13
DBMAT(I,J)=0.0 SUBP 14
DO 10 K=1,NROIJ SUBP 15
10 DBMAT(I,J)=DBMAT(I,J)+DMATX(I,K)*BJMAT(K,J) SUBP 16
C*** EVALUATE BIT*(D*BJ) SUBP 17
DO 20 J=1,NCOLJ SUBP 18
DO 20 I=1,NCOLI SUBP 19
SBSTF(I,J)=0.0 SUBP 20
DO 20 K=1,NROIJ SUBP 21
20 SBSTF(I,J)=SBSTF(I,J)+BIMAT(K,I)*DBMAT(K,J) SUBP 22
C*** ASSEMBLE SBSTF INTO ELEMENT STIFFNESS MATRIX SUBP 23
IFROW=0 SUBP 24
JFCOL=0 SUBP 25
IFROW=(INODE-1)*3+IFROW SUBP 26
JFCOL=(JNODE-1)*3+JFCOL SUBP 27
DO 30 I=1,NCOLI SUBP 28
IRSUB=IFROW+I SUBP 29
DO 30 J=1,NCOLJ SUBP 30
JCSUB=JFCOL+J SUBP 31
30 ESTIF(IRSUB,JCSUB)=ESTIF(IRSUB,JCSUB)+SBSTF(I,J)*DAREA SUBP 32
RETURN SUBP 33
END SUBP 34
```
# 9.5.16 Subroutines VZERO and ZEROMP
These routines simply set to zero the components of various vectors and arrays.
```txt
SUBROUTINE VZERO (NCOMP, VECTO) ZERO 1
C************************** ZERO 2
C ZERO 3
C*** ZEROES VECTOR VECTO ZERO 4
C ZERO 5
C************************** ZERO 6
DIMENSION VECTO(NCOMP) ZERO 7
DO 10 ICOMP=1,NCOMP ZERO 8
10 VECTO(ICOMP)=0.0 ZERO 9
RETURN ZERO 10
END ZERO 11
```
```txt
SUBROUTINE ZEROMP (EFFST,ELOAD,EPSTN,MELEM,MEVAB,MTOTG, ZERP 1
. MTOTV,MVFIX,NDOFN,NELEM,NEVAB,NGAUS, ZERP 2
. NTOTG,NTOTV,NVFIX,STRSG,TDISP,TFACT, ZERP 3
. TLOAD,TREAC) ZERP 4
C**************************ZERP 5
C ZERP 6
C*** ZERO EFFST,ELOAD,EPSTN,STRSG,TDISP,TFACT,TLOAD,TREAC ZERP 7
C ZERP 8
C**************************ZERP 9
DIMENSION ELOAD(MELEM,MEVAB),STRSG(5,MTOTG),TDISP(MTOTV), ZERP 10
. TLOAD(MELEM,MEVAB),TREAC(MVFIX,3),EPSTN(MTOTG), ZERP 11
. EFFST(MTOTG) ZERP 12
TFACT=0.0 ZERP 13
DO 30 IELEM=1,NELEM ZERP 14
DO 30 IEVAB=1,NEVAB ZERP 15
ELOAD(IELEM,IEVAB)=0.0 ZERP 16
30 TLOAD(IELEM,IEVAB)=0.0 ZERP 17
DO 40 ITOTV=1,NTOTV ZERP 18
40 TDISP(ITOTV)=0.0 ZERP 19
DO 50 IVFIX=1,NVFIX ZERP 20
DO 50 IDOFN=1,NDOFN ZERP 21
50 TREAC(IVFIX,IDOFN)=0.0 ZERP 22
DO 60 ITOTG=1,NTOTG ZERP 23
EPSTN(ITOTG)=0.0 ZERP 24
EFFST(ITOTG)=0.0 ZERP 25
DO 60 ISTR1=1,5 ZERP 26
60 STRSG(ISTR1,ITOTG)=0.0 ZERP 27
RETURN ZERP 28
END ZERP 29
```
# 9.6 Software for the layered approach
# 9.6.1 Overall program structure.
The overall program structure for the elasto-plastic Mindlin plate bending analysis program using the layered approach is given in Fig. 9.5. This program is named MINDLAY.
The program can solve problems of the same size as those solved by program MINDLIN. A maximum of 26 layers is allowed.
All new routines are now documented and these include: FEAM, DEPMPA, LAYMPA, MDMPA, OUTMPA, RESMPA, STIMPA and STRMPA. The outer routines, which have been described earlier, include ALGOR, BMATPB, CHECK1, CHECK2, ECHO, FRONT, INCREM, INPUT, JACOB2 and NODEXY.
The files which are used in the program are 5 (cardreader), 6 (lineprinter) and 1, 2, 3, 4, 8 (scratch files).
# 9.6.2 Subroutine FEAM
This routine organises the calling of the main routines in sequence.
