Esope source files

   1 : C THSIG2    SOURCE    CHAT      05/01/13    03:38:04     5004           
   2 :       SUBROUTINE THSIG2(SHP,DJAC,XSTRS,NBNN,LRE,REL,N,R)           
   3 : C=====================================================================          
   4 : C                                                                               
   5 : C   CALCULE SOMME H,I * SIG I J *H ,J 
   6 : C   MODE DE FOURIER                                          
   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    N             = Mode de Fourier
  16 : C    R             = rayon                      
  17 : C                                                                               
  18 : C  SORTIES                                                                      
  19 : C    REL(LRE,LRE)=MATRICE DE RIGIDITE                                           
  20 : C    JEANVOINE AOUT 95                                                            
  21 : C=======================================================================        
  22 : 
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8(A-H,O-Z)                                                        
  25 :       DIMENSION SHP(6,*),XSTRS(*),REL(LRE,*)  
  26 :                                              
  27 :       DO 100 IA=1,NBNN                                                          
  28 :         DO 200 IB=1,NBNN
  29 : C
  30 :           CC1
  31 :      1     = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  32 :      1     + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  33 :      1     + ((1+(N*N))/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  34 : C
  35 :           CC2
  36 :      1     = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  37 :      1     + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  38 :      1     + ((N*N)/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  39 : C
  40 :           CC3
  41 :      1     = SHP(2,IA)*XSTRS(1)*SHP(2,IB) + SHP(3,IA)*XSTRS(2)*SHP(3,IB)
  42 :      1     + SHP(2,IA)*XSTRS(4)*SHP(3,IB) + SHP(3,IA)*XSTRS(4)*SHP(2,IB)
  43 :      1     + ((1+(N*N))/(R**2))*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  44 : C
  45 :           CC4 = 2*N/(R*R)*SHP(1,IA)*XSTRS(3)*SHP(1,IB)
  46 : C
  47 :            REL(1+3*(IA - 1),1+3*(IB - 1)) = DJAC*CC1 + 
  48 :      1     REL(1+3*(IA - 1),1+3*(IB - 1))
  49 :            REL(2+3*(IA - 1),2+3*(IB - 1)) = DJAC*CC2 +
  50 :      1     REL(2+3*(IA - 1),2+3*(IB - 1))
  51 :            REL(3+3*(IA - 1),3+3*(IB - 1)) = DJAC*CC3 +
  52 :      1     REL(3+3*(IA - 1),3+3*(IB - 1))
  53 :            REL(3+3*(IA - 1),1+3*(IB - 1)) = DJAC*CC4 +
  54 :      1     REL(3+3*(IA - 1),1+3*(IB - 1))
  55 :            REL(1+3*(IA - 1),3+3*(IB - 1)) = DJAC*CC4 +
  56 :      1     REL(1+3*(IA - 1),3+3*(IB - 1))
  57 : C
  58 :   200   CONTINUE                   
  59 :   100 CONTINUE
  60 :                                                                 
  61 :       RETURN                                                                    
  62 :       END                                                                       
  63 :  
  64 :