Esope source files

   1 : C THSIG3    SOURCE    CHAT      05/01/13    03:38:08     5004           
   2 :       SUBROUTINE THSIG3(SHP,DJAC,XSTRS,NBNN,LRE,REL,R)           
   3 : C=====================================================================          
   4 : C                                                                               
   5 : C   CALCULE SOMME H,I * SIG I J *H ,J 
   6 : C   AXISYMETRIQUE                                          
   7 : C                                                                               
   8 : C  ENTREE                                                                       
   9 : C    SHP(6,NBNO)   = H   H,X   H,Y  H,Z                                         
  10 : C    DJAC          = JACOBIEN                                                    
  11 : C    XSTRS(NSTRS)  = CONTRAINTES                                                
  12 : C    NBNN          = NOMBRE DE NOEUDS                                                                                             
  13 : C    LRE           = TAILLE DE LA MATRICE DE RIGIDITE                            
  14 : C    REL(LRE,LRE)  = MATRICE DE RIGIDITE    (LRE=2*NBNN) 
  15 : C    R             = rayon                      
  16 : C                                                                               
  17 : C  SORTIES                                                                      
  18 : C    REL(LRE,LRE)=MATRICE DE RIGIDITE                                           
  19 : C    JEANVOINE AOUT 95                                                            
  20 : C=======================================================================        
  21 : 
  22 :       IMPLICIT INTEGER(I-N)
  23 :       IMPLICIT REAL*8(A-H,O-Z)                                                        
  24 :       DIMENSION SHP(6,*),XSTRS(*),REL(LRE,*)  
  25 :                                              
  26 :       DO 100 IA=1,NBNN                                                          
  27 :         DO 200 IB=1,NBNN
  28 : C
  29 :           CC1
  30 :      1     = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  31 :      1     + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  32 :      1     + (1.D0/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  33 : C
  34 :           CC2
  35 :      1     = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  36 :      1     + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  37 : C
  38 :           CC3
  39 :      1     = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  40 :      1     + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  41 : C
  42 :            REL(1+2*(IA - 1),1+2*(IB - 1)) = DJAC*CC1 + 
  43 :      1     REL(1+2*(IA - 1),1+2*(IB - 1))
  44 :            REL(2+2*(IA - 1),2+2*(IB - 1)) = DJAC*CC2 +
  45 :      1     REL(2+2*(IA - 1),2+2*(IB - 1))
  46 : C
  47 :   200   CONTINUE                   
  48 :   100 CONTINUE
  49 :                                                                 
  50 :       RETURN                                                                    
  51 :       END                                                                       
  52 :  
  53 :  
  54 :