flowpl.for

来自「有限元计算程序」· FOR 代码 · 共 33 行

C$DEBUG
C$LARGE
	SUBROUTINE FLOWPL(AVECT,ABETA,DVECT,NTYPE,PROPS,LPROP,
     .                  NSTR1,MMATS)                    !(Page 243)
C******************************************************************
C
C****THIS SUBROUTINE EVALUATES THE PLASTIC D VECTOR
C
C******************************************************************
	DIMENSION AVECT(4),DVECT(4),PROPS(MMATS,7)
	YOUNG=PROPS(LPROP,1)
	POISS=PROPS(LPROP,2)
	HARDS=PROPS(LPROP,6)
	FMUL1=YOUNG/(1.0+POISS)
	IF(NTYPE.EQ.1) GO TO 60
	FMUL2=YOUNG*POISS*(AVECT(1)+AVECT(2)+AVECT(4))/((1.0+POISS)*
     .      (1.0-2.0*POISS))
	DVECT(1)=FMUL1*AVECT(1)+FMUL2
	DVECT(2)=FMUL1*AVECT(2)+FMUL2
	DVECT(3)=0.5*AVECT(3)*YOUNG/(1.0+POISS)
	DVECT(4)=FMUL1*AVECT(4)+FMUL2
	GO TO 70
  60  FMUL3=YOUNG*POISS*(AVECT(1)+AVECT(2))/(1.0-POISS*POISS)
	DVECT(1)=FMUL1*AVECT(1)+FMUL3
	DVECT(2)=FMUL1*AVECT(2)+FMUL3
	DVECT(3)=0.5*AVECT(3)*YOUNG/(1.0+POISS)
	DVECT(4)=FMUL1*AVECT(4)+FMUL3
  70  DENOM=HARDS
	DO 80 ISTR1=1,NSTR1
  80  DENOM=DENOM+AVECT(ISTR1)*DVECT(ISTR1)
	ABETA=1.0/DENOM
	RETURN
	END