dbe.for

板弯曲分析有限元程序

C$DEBUG
C$LARGE
      SUBROUTINE DBE !(PP96)
C
C***  CALCULATES D X B
C
      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)
      DO 10 ISTRE=1,NSTRE
      DO 10 IEVAB=1,NEVAB
      DBMAT(ISTRE,IEVAB)=0.0
      DO 10 JSTRE=1,NSTRE
      DBMAT(ISTRE,IEVAB)=DBMAT(ISTRE,IEVAB)+
     * DMATX(ISTRE,JSTRE)*BMATX(JSTRE,IEVAB)
   10 CONTINUE
      RETURN
      END