gaussq.for

板弯曲分析有限元程序

C$DEBUG
C$LARGE
      SUBROUTINE GAUSSQ !(PP67)

      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
      IF(NGAUS.GT.2) GO TO 10
      POSGP(1)=-0.577350269189626
      WEIGP(1)=1.0
      GO TO 20
   10 POSGP(1)=-0.774596669241483
      POSGP(2)=0.0
      WEIGP(1)=0.5555555555555556
      WEIGP(2)=0.8888888888888889
   20 KGAUS=NGAUS/2
      DO 30 IGASH=1,KGAUS
      JGASH=NGAUS+1-IGASH
      POSGP(JGASH)=-POSGP(IGASH)
      WEIGP(JGASH)=WEIGP(IGASH)
   30 CONTINUE
      RETURN
      END