Sources Esope | Cast3M

bsig8
C BSIG8     SOURCE    CHAT      05/01/12    21:42:51     5004
      SUBROUTINE BSIG8(BGENE,NBNN,XSTRS,DJAC,XFORC)
C=======================================================================
C     CALCULE B*SIGMA    POUR LE COQ8 ET LE COQ6
C
C  ENTREES
C    BGENE(5,LRE)   : MATRICE B
C    NBNN           : NOMBRE DE NOEUDS
C    XSTRS(5)=CONTRAINTES   EN UN POINT DE GAUSS
C    XFORC(LRE)     : FORCES NODALES
C    DJAC           : JACOBIEN  AU NPOINT DE GAUSS
C  SORTIES
C    XFORC(LRE)     : FORCES NODALES
C            EBERSOLT   AOUT 86
C=======================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      PARAMETER(X577=.577350269189626D0,UNSIX=.166666666666666D0)
      PARAMETER(XZER=0.D0)
      DIMENSION BGENE(5,*),XSTRS(*),XFORC(*)
C
C
      DO 50 IA=1,NBNN
      IDEB=(IA-1)*6
      DO 60 IAA=1,6
      CC=XZER
C
      DO 100 IB=1,5
      CC=CC+BGENE(IB,IDEB+IAA)*XSTRS(IB)
  100 CONTINUE
      XFORC(IDEB+IAA)=XFORC(IDEB+IAA)+CC*DJAC
   60 CONTINUE
   50 CONTINUE
C
      RETURN
      END