Esope source files

   1 : C BSIG8     SOURCE    CHAT      05/01/12    21:42:51     5004           
   2 :       SUBROUTINE BSIG8(BGENE,NBNN,XSTRS,DJAC,XFORC)                             
   3 : C=======================================================================        
   4 : C     CALCULE B*SIGMA    POUR LE COQ8 ET LE COQ6                                
   5 : C                                                                               
   6 : C  ENTREES                                                                      
   7 : C    BGENE(5,LRE)   : MATRICE B                                                 
   8 : C    NBNN           : NOMBRE DE NOEUDS                                          
   9 : C    XSTRS(5)=CONTRAINTES   EN UN POINT DE GAUSS                                
  10 : C    XFORC(LRE)     : FORCES NODALES                                            
  11 : C    DJAC           : JACOBIEN  AU NPOINT DE GAUSS                              
  12 : C  SORTIES                                                                      
  13 : C    XFORC(LRE)     : FORCES NODALES                                            
  14 : C            EBERSOLT   AOUT 86                                                 
  15 : C=======================================================================        
  16 :       IMPLICIT INTEGER(I-N)
  17 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  18 :       PARAMETER(X577=.577350269189626D0,UNSIX=.166666666666666D0)               
  19 :       PARAMETER(XZER=0.D0)                                                      
  20 :       DIMENSION BGENE(5,*),XSTRS(*),XFORC(*)                                    
  21 : C                                                                               
  22 : C                                                                               
  23 :       DO 50 IA=1,NBNN                                                           
  24 :       IDEB=(IA-1)*6                                                             
  25 :       DO 60 IAA=1,6                                                             
  26 :       CC=XZER                                                                   
  27 : C                                                                               
  28 :       DO 100 IB=1,5                                                             
  29 :       CC=CC+BGENE(IB,IDEB+IAA)*XSTRS(IB)                                        
  30 :   100 CONTINUE                                                                  
  31 :       XFORC(IDEB+IAA)=XFORC(IDEB+IAA)+CC*DJAC                                   
  32 :    60 CONTINUE                                                                  
  33 :    50 CONTINUE                                                                  
  34 : C                                                                               
  35 :       RETURN                                                                    
  36 :       END                                                                       
  37 :