Télécharger la plaquette Cast3M


fr en

Télécharger tuycar.eso

Retour à la liste

Numérotation des lignes : 

tuycar

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

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