C BSIG      SOURCE    CHAT      05/01/12    21:42:58     5004

C=======================================================================
C=                            B S I G                                  =
C=                            -------                                  =
C=                                                                     =
C=  Fonction :                                                         =
C=  ----------                                                         =
C=  Calcul de la contribution d'un point de Gauss aux forces internes  =
C=  B*Sigma d'un element fini                                          =
C=                                                                     =
C=  Parametres :  (E)=Entree  (S)=Sortie                               =
C=  ------------                                                       =
C=   BGENE    (E)   Matrice de gradient (B) calculee au point de Gauss =
C=   XSTRS    (E)   Valeurs des contraintes au point de Gauss          =
C=   NSTRS    (E)   Nombre de composantes de contraintes/deformations  =
C=   LRE      (E)   Nombre de DDL de l'element fini                    =
C=   DJAC     (S)   Jacobien*poids au point de Gauss etudie            =
C=   FORCE   (E/S)  Forces internes aux noeuds de l'element            =
C=======================================================================

      SUBROUTINE BSIG (BGENE,XSTRS,NSTRS,LRE,DJAC,XFORC)

      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)

-INC CCREEL

      DIMENSION BGENE(NSTRS,*),XSTRS(*),XFORC(*)

      DO i=1,LRE
        cc=XZero
        DO j=1,NSTRS
          cc=cc+BGENE(j,i)*XSTRS(j)
        ENDDO
        XFORC(i)=cc*DJAC+XFORC(i)
      ENDDO

      RETURN
      END