Esope source files

   1 : C MTOKCE    SOURCE    BP208322  12/03/08    21:15:04     7294           
   2 : c       SUBROUTINE MTOKCE(LRE,NDDL,REL,VROT,RELB)                                
   3 :       SUBROUTINE MTOKCE(LRE,NDDL,REL,ROME,RELB)                                
   4 : C=======================================================================        
   5 : C                                                                               
   6 : C  Calcul de la matrice de raideur centrifuge à partir dela matrice de masse   
   7 : C  pour les éléments massifs en 3D
   8 : C
   9 : C     Didier COMBESCURE Mars 2003 
  10 : C     BP on sort le remplissage de ROME (cst) de la boucle (cf KCENT3.eso)
  11 : C                                                                               
  12 : C=======================================================================        
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  15 : c       DIMENSION REL(LRE,1),RELB(LRE,1),VROT(*),ROME(3,3)
  16 :       DIMENSION REL(LRE,1),RELB(LRE,1),ROME(3,3)
  17 : C
  18 :       NNOEU = LRE/NDDL
  19 : C
  20 : c       ROME(1,1) = (-1.)*((VROT(2)**2) + (VROT(3)**2))
  21 : c       ROME(2,2) = (-1.)*((VROT(1)**2) + (VROT(3)**2))
  22 : c       ROME(3,3) = (-1.)*((VROT(1)**2) + (VROT(2)**2))
  23 : c       ROME(1,2) = VROT(1)*VROT(2)
  24 : c       ROME(1,3) = VROT(1)*VROT(3)
  25 : c       ROME(2,3) = VROT(2)*VROT(3)
  26 : c       ROME(2,1) = ROME(1,2) 
  27 : c       ROME(3,1) = ROME(1,3) 
  28 : c       ROME(3,2) = ROME(2,3) 
  29 : C
  30 : C                                 
  31 :       DO   100 IN=1,NNOEU
  32 : C
  33 :         DO   100 IM=1,NNOEU
  34 : C
  35 :           DO   100 IDD1=1,NDDL
  36 : C
  37 :             DO   100 IDD2=1,NDDL
  38 : C
  39 :             RELB((3*IN)-3+IDD1,(3*IM)-3+IDD2)=
  40 :      . REL((3*IN)-3+IDD1,(3*IM)-3+IDD1)*ROME(IDD1,IDD2) 
  41 : C 
  42 :   100      CONTINUE     
  43 : C                                      
  44 :       RETURN                                                                    
  45 :       END    
  46 :                                         
  47 :  
  48 :  
  49 :