increm.for

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

C$DEBUG
C$LARGE
	SUBROUTINE INCREM(ELOAD,FIXED,IINCS,MELEM,MEVAB,MITER,
     .                  MTOTV,MVFIX,NDOFN,NELEM,NEVAB,NOUTP,
     .                  NOFIX,NTOTV,NVFIX,PRESC,RLOAD,TFACT,
     .                  TLOAD,TOLER)             !(Page 212)
C********************************************************************
C
C*** THIS SUBROUTINE INCREMENTS THE APPLIED LOADING
C
C********************************************************************
      DIMENSION ELOAD(MELEM,MEVAB),FIXED(MTOTV),IFFIX(MTOTV),
     .          NOUTP(2),NOFIX(MVFIX),
     .          PRESC(MVFIX,NDOFN),RLOAD(MELEM,MEVAB),TLOAD(MELEM,MEVAB)
	WRITE(6,900) IINCS
 900  FORMAT(1H 5X,17HINCREMENT NUMBER ,I5)
	READ(5,*) FACTO,TOLER,MITER,NOUTP(1),NOUTP(2)
 950  FORMAT(2F10.5,3I5)
	TFACT=TFACT+FACTO
	WRITE(6,960)TFACT,TOLER,MITER,NOUTP(1),NOUTP(2)
 960  FORMAT(1H 5X,13HLOAD FACTOR =,F10.5,5X,
     .24H CONVERGENCE TOLERANCE=,F10.5,//,5X,24HMAX. NO. OF ITERATIONS =,
     . I5,//27H INITIAL OUTPUT PARAMETER =,I5,5X,24HFINAL OUTPUT PARAMET
     .ER =,I5)
	DO 80 IELEM=1,NELEM
	DO 80 IEVAB=1,NEVAB
	ELOAD(IELEM,IEVAB)=ELOAD(IELEM,IEVAB)+RLOAD(IELEM,IEVAB)*FACTO
  80  TLOAD(IELEM,IEVAB)=TLOAD(IELEM,IEVAB)+RLOAD(IELEM,IEVAB)*FACTO
C
C*** INTERPRET FIXITY DATA IN VECTOR FORM
C
	DO 100 ITOTV=1,NTOTV
 100  FIXED(ITOTV)=0.0
	DO 110 IVFIX=1,NVFIX
	NLOCA=(NOFIX(IVFIX)-1)*NDOFN
	DO 110 IDOFN=1,NDOFN
	NGASH=NLOCA+IDOFN
	FIXED(NGASH)=PRESC(IVFIX,IDOFN)*FACTO
 110  CONTINUE
	RETURN
	END