Télécharger la plaquette Cast3M


fr en

Télécharger tuycag.eso

Retour à la liste

Numérotation des lignes : 

tuycag

  1. C TUYCAG    SOURCE    CHAT      05/01/13    03:57:06     5004
  2.       SUBROUTINE TUYCAG(CARAC,KERRE,ICAS)
  3. C-----------------------------------------------------------------------
  4. C  ROUTINE DE CALCUL DES CARACTERISTIQUES GEOMETRIQUES DE L ELEMENT DE
  5. C  TUYAU DROIT OU COUDE - D'APRES TUYCAR (APPELE PAR GYRO2)
  6. C
  7. C  ENTREE
  8. C          CARAC   TABLEAU DE CARACTERISTIQUES GEOMETRIQUES
  9. C                  DONNEES PAR L UTILISATEUR
  10. C          ICAS    ON CALCULE PLUS OU MOINS DE CARACTERISTIQUES
  11. C                  SELON LES CAS :
  12. C                  0 = VECTEUR LOCAL
  13. C                  1 = 0 + SECTION + INERTIES NON MODIFIEES
  14. C                  2 = 0 + SECTION + INERTIES MODIFIEES
  15. C
  16. C  SORTIE
  17. C          CARAC   CARACTERISTIQUES DE LA POUTRE EQUIVALENTE AU TUYAU
  18. C          KERRE   CODE D ERREUR  O SI PAS DE PB
  19. C
  20. C    COMBESCURE Didier Mars 2003
  21. C-----------------------------------------------------------------------
  22.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8(A-H,O-Z)
  24. -INC CCREEL
  25.  
  26. -INC PPARAM
  27. -INC CCOPTIO
  28. *
  29.       DIMENSION CARAC(*)
  30.       KERRE=0
  31. C
  32. C  SI 'RACO' EST PRESENT, IL FAUT UN VECTEUR
  33. C
  34.       IF(CARAC(3).GT.1.D-9.AND.ABS(CARAC(4)).LT.1.D-9.AND.
  35.      1     ABS(CARAC(5)).LT.1.D-9.AND.ABS(CARAC(6)).LT.1.D-9) THEN
  36.          MOTERR(1:4)='VECT'
  37.          MOTERR(5:8)='CARA'
  38.          KERRE=77
  39.          RETURN
  40.       ENDIF
  41. C
  42. C  VECTEUR LOCAL
  43. C
  44.       CARAC(8)=CARAC(5)
  45.       CARAC(9)=CARAC(6)
  46.       CARAC(10)=CARAC(7)
  47. C
  48. C  VITESSE DE ROTATION
  49. C
  50.       CARAC(7)=CARAC(4)
  51.       IF(ICAS.EQ.0) RETURN
  52. C
  53. C  CALCUL DES SECTION ET INERTIES
  54. C  ON MET LES SECTIONS REDUITES A 0.D0
  55. C
  56.       CARAC(5)=0.D00
  57.       CARAC(6)=0.D00
  58.       EPAIS=CARAC(1)
  59.       REXT=CARAC(2)
  60.       RINT=REXT-EPAIS
  61.       RACO=CARAC(3)
  62.       XIN=XPI*(REXT**4-RINT**4)*0.25D00
  63.       CARAC(4)=XPI*(REXT**2-RINT**2)
  64.       CARAC(1)=2.D00*XIN
  65.       IF(ICAS.EQ.2) GO TO 20
  66.       CARAC(2)=XIN
  67.       CARAC(3)=XIN
  68.       RETURN
  69.  20   CONTINUE
  70.       XK=1.D00
  71.       IF(RACO.EQ.0.D00) GO TO 10
  72.       XLAM=EPAIS*RACO*4.D00/(REXT+RINT)**2
  73.       IF(XLAM.EQ.0.D00) GO TO 999
  74.       XK=1.65D00/XLAM
  75.       IF(XK.LT.1.D00) XK=1.D00
  76.  10   CONTINUE
  77.       CARAC(2)=XIN/XK
  78.       CARAC(3)=CARAC(2)
  79.       RETURN
  80.  999  KERRE=1
  81.       RETURN
  82.       END
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  

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