modpb.for

板弯曲分析有限元程序

C$DEBUG
C$LARGE
      SUBROUTINE MODPS(LPROP) !(PP141)

      COMMON/CONTRO/NPOIN,NELEM,NNODE,NDOFN,NDIME,
     *NSTRE,NTYPE,NGAUS,NPROP,NMATS,NVFIX,NEVAB,
     *ICASE,NCASE,ITEMP,IPROB,NPROB
	COMMON/LGDATA/COORD(MPOIN,MDIME),PROPS(MMATS,
     *MPROP),PRESC(MVFIX,MDOFN),ASDIS(MTOTV),ELOAD
     *(MELEM,MEVAB),NOFIX(MVFIX),IFPRE(MVFIX,MDOFN),
     *LNODS(MELEM,MNODE),MATNO(MELEM)
	COMMON/WORK/ELCOD(MDIME,MNODE),SHAPE(MNODE),
     *DERIV(MDIME,MNODE),DMATX(MSTRE,MSTRE),CARTD
     *(MDIME,MNODE),DBMAT(MSTRE,MEVAB),BMATX(MSTRE,
     *MEVAB),SMATX(MSTRE,MEVAB,MGASP),POSGP(MGAUS),
     *WEIGP(MGAUS),GPCOD(MDIME,MGASP),NEROR(24)
C
      DO 10 ISTRE=1,NSTRE
      DO 10 JSTRE=1,NSTRE
	DMATX(ISTRE,JSTRE)=0.0
   10 CONTINUE
      YOUNG=PROPS(LPROP,1)
      POISS=PROPS(LPROP,2)
      THICK=PROPS(LPROP,3)
	DMATX(1,1)=YOUNG*THICK*THICK*THICK
     */(12.0*(1.0-POISS*POISS))
	DMATX(1,2)=POISS*DMATX(1,1)
	DMATX(2,2)=DMATX(1,1)
	DMATX(2,1)=DMATX(1,2)
      DMATX(3,3)=(1.0-POISS)*DMATX(1,1)/2.0
      DMATX(4,4)=YOUNG*THICK/(2.4*(1.0+POISS))
	DMATX(5,5)=DMATX(4,4)
	RETURN
	END