Download bsig.eso

Back to the list

   1 : C BSIG      SOURCE    CHAT      05/01/12    21:42:58     5004           
   2 : 
   3 : C=======================================================================
   4 : C=                            B S I G                                  =
   5 : C=                            -------                                  =
   6 : C=                                                                     =
   7 : C=  Fonction :                                                         =
   8 : C=  ----------                                                         =
   9 : C=  Calcul de la contribution d'un point de Gauss aux forces internes  =
  10 : C=  B*Sigma d'un element fini                                          =
  11 : C=                                                                     =
  12 : C=  Parametres :  (E)=Entree  (S)=Sortie                               =
  13 : C=  ------------                                                       =
  14 : C=   BGENE    (E)   Matrice de gradient (B) calculee au point de Gauss =
  15 : C=   XSTRS    (E)   Valeurs des contraintes au point de Gauss          =
  16 : C=   NSTRS    (E)   Nombre de composantes de contraintes/deformations  =
  17 : C=   LRE      (E)   Nombre de DDL de l'element fini                    =
  18 : C=   DJAC     (S)   Jacobien*poids au point de Gauss etudie            =
  19 : C=   FORCE   (E/S)  Forces internes aux noeuds de l'element            =
  20 : C=======================================================================
  21 : 
  22 :       SUBROUTINE BSIG (BGENE,XSTRS,NSTRS,LRE,DJAC,XFORC)
  23 : 
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8 (A-H,O-Z)
  26 : 
  27 : -INC CCREEL
  28 : 
  29 :       DIMENSION BGENE(NSTRS,*),XSTRS(*),XFORC(*)
  30 : 
  31 :       DO i=1,LRE
  32 :         cc=XZero
  33 :         DO j=1,NSTRS
  34 :           cc=cc+BGENE(j,i)*XSTRS(j)
  35 :         ENDDO
  36 :         XFORC(i)=cc*DJAC+XFORC(i)
  37 :       ENDDO
  38 : 
  39 :       RETURN
  40 :       END
  41 : 
  42 :  
  43 :  

© Cast3M 2003 - All rights reserved.
Disclaimer