Esope source files

   1 : C THSIG1    SOURCE    CHAT      05/01/13    03:38:00     5004           
   2 :       SUBROUTINE THSIG1(SHP,DJAC,XSTRS,NBNN,LRE,REL,R)           
   3 : C=====================================================================          
   4 : C                                                                               
   5 : C   CALCULE SOMME H,I * SIG I J *H ,J 
   6 : C   MODE DE FOURIER 0                                          
   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=3*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 :      1     + (1.D0/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  42 : C
  43 :            REL(1+3*(IA - 1),1+3*(IB - 1)) = DJAC*CC1 + 
  44 :      1     REL(1+3*(IA - 1),1+3*(IB - 1))
  45 :            REL(2+3*(IA - 1),2+3*(IB - 1)) = DJAC*CC2 +
  46 :      1     REL(2+3*(IA - 1),2+3*(IB - 1))
  47 :            REL(3+3*(IA - 1),3+3*(IB - 1)) = DJAC*CC3 +
  48 :      1     REL(3+3*(IA - 1),3+3*(IB - 1))
  49 : C
  50 :   200   CONTINUE                   
  51 :   100 CONTINUE
  52 :                                                                 
  53 :       RETURN                                                                    
  54 :       END                                                                       
  55 :  
  56 :  
  57 :