Esope source files

   1 : C MTOGYR    SOURCE    CHAT      05/01/13    01:53:58     5004           
   2 :       SUBROUTINE MTOGYR(LRE,NDDL,REL,VROT,RELB)                                
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C  Calcul de la matrice de raideur centrifuge à partir de la matrice de masse   
   6 : C  pour les éléments massifs en 3D
   7 : C
   8 : C     Didier COMBESCURE Mars 2003                                                
   9 : C                                                                               
  10 : C=======================================================================        
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  13 :       DIMENSION REL(LRE,1),RELB(LRE,1),VROT(*),ROME(3,3)
  14 : C
  15 :       NNOEU = LRE/NDDL
  16 : C
  17 : C
  18 :       ROME(1,1) = 0.
  19 :       ROME(2,2) = 0.
  20 :       ROME(3,3) = 0.
  21 :       ROME(1,2) = (-1.)*VROT(3)
  22 :       ROME(1,3) = VROT(2)
  23 :       ROME(2,3) = (-1.)*VROT(1)
  24 :       ROME(2,1) = VROT(3)
  25 :       ROME(3,1) = (-1.)*VROT(2)
  26 :       ROME(3,2) = VROT(1)
  27 : C
  28 : C                                 
  29 :       DO   100 IN=1,NNOEU
  30 : C
  31 :         DO   100 IM=1,NNOEU
  32 : C
  33 :           DO   100 IDD1=1,NDDL
  34 : C
  35 :             DO   100 IDD2=1,NDDL
  36 : C
  37 :             RELB((3*IN)-3+IDD1,(3*IM)-3+IDD2)=
  38 :      . REL((3*IN)-3+IDD2,(3*IM)-3+IDD2)*ROME(IDD1,IDD2) 
  39 : C 
  40 :   100      CONTINUE    
  41 : C                                   
  42 :       RETURN                                                                    
  43 :       END    
  44 :                                         
  45 :