Télécharger la plaquette Cast3M


fr en

Télécharger tuykar.eso

Retour à la liste

Numérotation des lignes : 

tuykar

  1. C TUYKAR    SOURCE    CHAT      05/01/13    03:57:16     5004
  2.       SUBROUTINE TUYKAR(CARAC,KERRE,ICAS,YOUNG)
  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          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    JUIN 85  MILLARD
  21. C    NOVEMBRE 89 EFFET DE PRESSION SUR COUDE M. PETIT
  22. C-----------------------------------------------------------------------
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25. -INC CCREEL
  26.  
  27. -INC PPARAM
  28. -INC CCOPTIO
  29. *
  30.       DIMENSION CARAC(*)
  31.       KERRE=0
  32. *
  33.       IF(CARAC(3).GT.1.D-9.AND.ABS(CARAC(4)).LT.1.D-9.AND.
  34.      1     ABS(CARAC(5)).LT.1.D-9.AND.ABS(CARAC(6)).LT.1.D-9) THEN
  35.          MOTERR(1:4)='VECT'
  36.          MOTERR(5:8)='CARA'
  37.          KERRE=77
  38.          RETURN
  39.       ENDIF
  40.       PRES=CARAC(7)
  41.       CISA=CARAC(8)
  42. C
  43. C  VECTEUR LOCAL
  44. C
  45.       CARAC(7)=CARAC(4)
  46.       CARAC(8)=CARAC(5)
  47.       CARAC(9)=CARAC(6)
  48.       IF(ICAS.EQ.0) RETURN
  49. C
  50. C  CALCUL DES SECTION ET INERTIES
  51. C  AM 29/4/94  ON CALCULE LES SECTIONS REDUITES
  52. C              ET ON LES MULTIPLIE PAR CISA
  53. C
  54.       EPAIS=CARAC(1)
  55.       REXT=CARAC(2)
  56.       RINT=REXT-EPAIS
  57.       RACO=CARAC(3)
  58.       XIN=XPI*(REXT**4-RINT**4)*0.25D00
  59.       CARAC(4)=XPI*(REXT**2-RINT**2)
  60.       CARAC(1)=2.D00*XIN
  61.       CARAC(5)=CARAC(4)*0.5D0*CISA
  62.       CARAC(6)=CARAC(5)
  63. *
  64.       IF(ICAS.EQ.2) GO TO 20
  65. *
  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(PRES.EQ.0.D00) GO TO 30
  76. C
  77. C     EFFET DE LA PRESSION SUR LE COEF DE FLEXIBILITE
  78. C
  79.       RMOY=REXT-EPAIS/2.D0
  80.       XX=6.D0*(RMOY/EPAIS)**(4.D0/3.D0)*(RACO/RMOY)**(1.D0/3.D0)
  81.       XCO=PRES*RMOY*XX/EPAIS/YOUNG
  82.       XK=XK/(1.D0+XCO)
  83.  30   CONTINUE
  84.       IF(XK.LT.1.D00) XK=1.D00
  85.  10   CONTINUE
  86.       CARAC(2)=XIN/XK
  87.       CARAC(3)=CARAC(2)
  88.       RETURN
  89.  999  KERRE=1
  90.       RETURN
  91.       END
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  

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