Esope source files

   1 : C BARKCE    SOURCE    KK2000    05/05/25    21:15:06     5088           
   2 :       SUBROUTINE BARKCE(REL,LRE,CARAC,COORD,VROT,WORK,KERRE)                                      
   3 : C-----------------------------------------------------------------------        
   4 : C  ROUTINE DE CALCUL DE LA MATRICE DE RAIDEUR CENTRIFUGE ELEMENTAIRE 
   5 : C  DE L ELEMENT DE BARRE (REPERE TOURNANT)                                                                    
   6 : C                                                                               
   7 : C  ENTREE                                                                       
   8 : C          LRE     DIMENSION DU TABLEAU REL                                     
   9 : C          CARAC   CARACTERISTIQUE MATERIAU                                                  
  10 : C          COORD   COORDONNEES DES 2 NOEUDS    
  11 : C          VROT    VECTEUR VITESSE DE ROTATION                                 
  12 : C                                                                               
  13 : C  SORTIE                                                                       
  14 : C          REL     MATRICE DE COUPLAGE ELEMENTAIRE     
  15 : C
  16 : C      Didier COMBESCURE mars 2003
  17 : C                            
  18 : C-----------------------------------------------------------------------                
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  21 : -INC CCREEL                                                                   
  22 :       DIMENSION REL(LRE,*),CARAC(*),COORD(3,*),V(3),WORK(*),VROT(*) 
  23 :       DIMENSION ROME(3,3)                                     
  24 : C
  25 :       SEC0=CARAC(1)
  26 :       RHO=CARAC(2)
  27 : C 
  28 :       DO 1 I=1,36
  29 :   1   REL(I,1)=0.D00
  30 : C                                                                              
  31 :       V(3)=0.                                                                   
  32 :       DO 2 I=1,3                                                             
  33 :   2   V(I)=COORD(I,2)-COORD(I,1)                                                
  34 :       XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3)) 
  35 : C      WRITE(*,*) XL                                       
  36 :       XL=SQRT(XL)                                                               
  37 : C                                                                               
  38 :       FAC=SEC0*RHO*XL/6.                                                             
  39 :       FAC2=FAC*2. 
  40 : C 
  41 :       ROME(1,1) = (-1.)*((VROT(2)**2) + (VROT(3)**2))
  42 :       ROME(2,2) = (-1.)*((VROT(1)**2) + (VROT(3)**2))
  43 :       ROME(3,3) = (-1.)*((VROT(1)**2) + (VROT(2)**2))
  44 :       ROME(1,2) = VROT(1)*VROT(2)
  45 :       ROME(1,3) = VROT(1)*VROT(3)
  46 :       ROME(2,3) = VROT(2)*VROT(3)
  47 :       ROME(2,1) = ROME(1,2) 
  48 :       ROME(3,1) = ROME(1,3) 
  49 :       ROME(3,2) = ROME(2,3) 
  50 : C
  51 : C  REMPLISSAGE DE LA MATRICE
  52 : C  FX1
  53 :       REL(1,1)=FAC2*ROME(1,1)
  54 :       REL(1,2)=FAC2*ROME(1,2)
  55 :       REL(1,3)=FAC2*ROME(1,3)
  56 :       REL(1,4)=FAC*ROME(1,1)
  57 :       REL(1,5)=FAC*ROME(1,2)
  58 :       REL(1,6)=FAC*ROME(1,3)
  59 : C  FY1
  60 :       REL(2,1)=FAC2*ROME(2,1)
  61 :       REL(2,2)=FAC2*ROME(2,2)
  62 :       REL(2,3)=FAC2*ROME(2,3)
  63 :       REL(2,4)=FAC*ROME(2,1)
  64 :       REL(2,5)=FAC*ROME(2,2)
  65 :       REL(2,6)=FAC*ROME(2,3)
  66 : C  FZ1
  67 :       REL(3,1)=FAC2*ROME(3,1)
  68 :       REL(3,2)=FAC2*ROME(3,2)
  69 :       REL(3,3)=FAC2*ROME(3,3)
  70 :       REL(3,4)=FAC*ROME(3,1)
  71 :       REL(3,5)=FAC*ROME(3,2)
  72 :       REL(3,6)=FAC*ROME(3,3)
  73 : C  FX2
  74 :       REL(4,1)=FAC*ROME(1,1)
  75 :       REL(4,2)=FAC*ROME(1,2)
  76 :       REL(4,3)=FAC*ROME(1,3)
  77 :       REL(4,4)=FAC2*ROME(1,1)
  78 :       REL(4,5)=FAC2*ROME(1,2)
  79 :       REL(4,6)=FAC2*ROME(1,3)
  80 : C  FY2
  81 :       REL(5,1)=FAC*ROME(2,1)
  82 :       REL(5,2)=FAC*ROME(2,2)
  83 :       REL(5,3)=FAC*ROME(2,3)
  84 :       REL(5,4)=FAC2*ROME(2,1)
  85 :       REL(5,5)=FAC2*ROME(2,2)
  86 :       REL(5,6)=FAC2*ROME(2,3)
  87 : C  FZ2
  88 :       REL(6,1)=FAC2*ROME(3,1)
  89 :       REL(6,2)=FAC2*ROME(3,2)
  90 :       REL(6,3)=FAC2*ROME(3,3)
  91 :       REL(6,4)=FAC*ROME(3,1)
  92 :       REL(6,5)=FAC*ROME(3,2)
  93 :       REL(6,6)=FAC*ROME(3,3)
  94 : C
  95 :       KERRE = 0
  96 : C                                                      
  97 :       RETURN                                                                    
  98 :       END                                                                       
  99 :  
 100 :  
 101 :