zero.for

本程序是D.R.J.OWEN主编的FINITE ELEMENTS IN PLASTICITYG一书中第八章例题的有限元程序

C$DEBUG
C$LARGE
	SUBROUTINE ZERO(ELOAD,MELEM,MEVAB,MPOIN,MTOTG,MTOTV,NDOFN,NELEM,
     .          NEVAB,NGAUS,NSTR1,NTOTG,NTOTV,NVFIX,STRSG,TDISP,
     .          VIVEL,VISTN,TTIME,TLOAD,TREAC,MVFIX) !(Page 297)
C*********************************************************************
C
C*** THIS SUBROUTINE INITIALISES VARIOUS ARRAYS TO ZERO
C
C*********************************************************************
	DIMENSION ELOAD(MELEM,MEVAB),STRSG(4,MTOTG),TDISP(MTOTV),
     .          TLOAD(MELEM,MEVAB),TREAC(MVFIX,2),VIVEL(5,MTOTG),
     .          VISTN(4,MTOTG)
	TTIME=0.0
	TFACT=0.0
	DO 30 IELEM=1,NELEM
	DO 30 IEVAB=1,NEVAB
	ELOAD(IELEM,IEVAB)=0.0
  30  TLOAD(IELEM,IEVAB)=0.0
	DO 40 ITOTV=1,NTOTV
  40  TDISP(ITOTV)=0.0
	DO 50 IVFIX=1,NVFIX
	DO 50 IDOFN=1,NDOFN
  50  TREAC(IVFIX,IDOFN)=0.0
	DO 60 ITOTG=1,NTOTG
	VIVEL(5,ITOTG)=0.0
	DO 60 ISTR1=1,NSTR1
	VISTN(ISTR1,ITOTG)=0.0
	VIVEL(ISTR1,ITOTG)=0.0
  60  STRSG(ISTR1,ITOTG)=0.0
	RETURN
	END