Télécharger la plaquette Cast3M


fr en

Télécharger barkce.eso

Retour à la liste

Numérotation des lignes : 

barkce

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

© Cast3M 2003 - Tous droits réservés.
Mentions légales