Esope source files

   1 : C BARCOR    SOURCE    KK2000    05/05/25    21:15:02     5088           
   2 :       SUBROUTINE BARCOR(REL,LRE,CARAC,COORD,VROT,WORK,KERRE)                                      
   3 : C-----------------------------------------------------------------------        
   4 : C  ROUTINE DE CALCUL DE LA MATRICE DE COUPLAGE GYROSCOPIQUE 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 : C
  29 :       DO 1 I=1,36
  30 :   1   REL(I,1)=0.D00
  31 : C                                                                              
  32 :       V(3)=0.                                                                   
  33 :       DO 2 I=1,3                                                             
  34 :   2   V(I)=COORD(I,2)-COORD(I,1)                                                
  35 :       XL=(V(1)*V(1)+V(2)*V(2)+V(3)*V(3)) 
  36 : C      WRITE(*,*) XL                                       
  37 :       XL=SQRT(XL)                                                               
  38 : C                                                                               
  39 :       FAC=SEC0*RHO*XL/6.                                                             
  40 :       FAC2=FAC*2. 
  41 : C 
  42 :       ROME(1,1) = 0.
  43 :       ROME(2,2) = 0.
  44 :       ROME(3,3) = 0.
  45 :       ROME(1,2) = (-1.)*VROT(3)
  46 :       ROME(1,3) = VROT(2)
  47 :       ROME(2,3) = (-1.)*VROT(1)
  48 :       ROME(2,1) = VROT(3)
  49 :       ROME(3,1) = (-1.)*VROT(2)
  50 :       ROME(3,2) = VROT(1)
  51 : C
  52 : C  REMPLISSAGE DE LA MATRICE
  53 : C  FX1
  54 :       REL(1,1)=FAC2*ROME(1,1)
  55 :       REL(1,2)=FAC2*ROME(1,2)
  56 :       REL(1,3)=FAC2*ROME(1,3)
  57 :       REL(1,4)=FAC*ROME(1,1)
  58 :       REL(1,5)=FAC*ROME(1,2)
  59 :       REL(1,6)=FAC*ROME(1,3)
  60 : C  FY1
  61 :       REL(2,1)=FAC2*ROME(2,1)
  62 :       REL(2,2)=FAC2*ROME(2,2)
  63 :       REL(2,3)=FAC2*ROME(2,3)
  64 :       REL(2,4)=FAC*ROME(2,1)
  65 :       REL(2,5)=FAC*ROME(2,2)
  66 :       REL(2,6)=FAC*ROME(2,3)
  67 : C  FZ1
  68 :       REL(3,1)=FAC2*ROME(3,1)
  69 :       REL(3,2)=FAC2*ROME(3,2)
  70 :       REL(3,3)=FAC2*ROME(3,3)
  71 :       REL(3,4)=FAC*ROME(3,1)
  72 :       REL(3,5)=FAC*ROME(3,2)
  73 :       REL(3,6)=FAC*ROME(3,3)
  74 : C  FX2
  75 :       REL(4,1)=FAC*ROME(1,1)
  76 :       REL(4,2)=FAC*ROME(1,2)
  77 :       REL(4,3)=FAC*ROME(1,3)
  78 :       REL(4,4)=FAC2*ROME(1,1)
  79 :       REL(4,5)=FAC2*ROME(1,2)
  80 :       REL(4,6)=FAC2*ROME(1,3)
  81 : C  FY2
  82 :       REL(5,1)=FAC*ROME(2,1)
  83 :       REL(5,2)=FAC*ROME(2,2)
  84 :       REL(5,3)=FAC*ROME(2,3)
  85 :       REL(5,4)=FAC2*ROME(2,1)
  86 :       REL(5,5)=FAC2*ROME(2,2)
  87 :       REL(5,6)=FAC2*ROME(2,3)
  88 : C  FZ2
  89 :       REL(6,1)=FAC2*ROME(3,1)
  90 :       REL(6,2)=FAC2*ROME(3,2)
  91 :       REL(6,3)=FAC2*ROME(3,3)
  92 :       REL(6,4)=FAC*ROME(3,1)
  93 :       REL(6,5)=FAC*ROME(3,2)
  94 :       REL(6,6)=FAC*ROME(3,3)
  95 : C
  96 : C      WRITE(*,*) 'barcor.eso'
  97 : C      WRITE(*,*) FAC
  98 :       KERRE = 0
  99 : C                                                      
 100 :       RETURN                                                                    
 101 :       END                                                                       
 102 :  
 103 :  
 104 :