bmatpb.for

板弯曲分析有限元程序

C$DEBUG
C$LARGE
      SUBROUTINE BMATPS !(PP139)
 
      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
	BMATX9ISTRE,IEVAB)=0.0
   10 CONTINUE
      JGASH=0
	DO 20 INODE=1NNODE
	IGASH=JGASH+1
	BMATX(4,IGASH)=CARTD(1,INODE)
	BMATX(5,IGASH)=CARTD(2,INODE)
	IGASH=IGASH+1
	JGASH=JGASH+1
	BMATX(1,IGASH)=-CARTD(1,INODE)
	BMATX(3,IGASH)=-CARTD(2,INODE)
	BMATX(4,IGASH)=-SHAPE(INODE)
	BMATX(2,JGASH)=-CARTD(2,INODE)
	BMATX(3,JGASH)=-CARTD(1,INODE)
	BMATX(5,JGASH)=-SHAPE(INODE)
   20 CONTINUE
      RETURN
	END





      NGASH=0
      DO 10 INODE=1,NNODE
      MGASH=NGASH+1
      NGASH=MGASH+1
      BMATX(1,MGASH)=CARTD(1,INODE)
      BMATX(1,NGASH)=0.0
      BMATX(2,MGASH)=0.0
      BMATX(2,NGASH)=CARTD(2,INODE)
      BMATX(3,MGASH)=CARTD(2,INODE)
      BMATX(3,NGASH)=CARTD(1,INODE)
   10 CONTINUE
      RETURN
      END