Esope source files

   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 : -INC CCOPTIO
  27 : *     
  28 :       DIMENSION CARAC(*)
  29 :       KERRE=0
  30 : *     
  31 :       IF(CARAC(3).GT.1.D-9.AND.ABS(CARAC(4)).LT.1.D-9.AND.
  32 :      1     ABS(CARAC(5)).LT.1.D-9.AND.ABS(CARAC(6)).LT.1.D-9) THEN
  33 :          MOTERR(1:4)='VECT'
  34 :          MOTERR(5:8)='CARA'
  35 :          KERRE=77
  36 :          RETURN
  37 :       ENDIF
  38 :       PRES=CARAC(7)
  39 :       CISA=CARAC(8)
  40 : C     
  41 : C  VECTEUR LOCAL
  42 : C     
  43 :       CARAC(7)=CARAC(4)
  44 :       CARAC(8)=CARAC(5)
  45 :       CARAC(9)=CARAC(6)
  46 :       IF(ICAS.EQ.0) RETURN
  47 : C     
  48 : C  CALCUL DES SECTION ET INERTIES
  49 : C  AM 29/4/94  ON CALCULE LES SECTIONS REDUITES
  50 : C              ET ON LES MULTIPLIE PAR CISA
  51 : C     
  52 :       EPAIS=CARAC(1)
  53 :       REXT=CARAC(2)
  54 :       RINT=REXT-EPAIS
  55 :       RACO=CARAC(3)
  56 :       XIN=XPI*(REXT**4-RINT**4)*0.25D00
  57 :       CARAC(4)=XPI*(REXT**2-RINT**2)
  58 :       CARAC(1)=2.D00*XIN
  59 :       CARAC(5)=CARAC(4)*0.5D0*CISA
  60 :       CARAC(6)=CARAC(5)
  61 : *     
  62 :       IF(ICAS.EQ.2) GO TO 20
  63 : *     
  64 :       CARAC(2)=XIN
  65 :       CARAC(3)=XIN
  66 :       RETURN
  67 :  20   CONTINUE
  68 :       XK=1.D00
  69 :       IF(RACO.EQ.0.D00) GO TO 10
  70 :       XLAM=EPAIS*RACO*4.D00/(REXT+RINT)**2
  71 :       IF(XLAM.EQ.0.D00) GO TO 999
  72 :       XK=1.65D00/XLAM
  73 :       IF(PRES.EQ.0.D00) GO TO 30
  74 : C     
  75 : C     EFFET DE LA PRESSION SUR LE COEF DE FLEXIBILITE
  76 : C     
  77 :       RMOY=REXT-EPAIS/2.D0
  78 :       XX=6.D0*(RMOY/EPAIS)**(4.D0/3.D0)*(RACO/RMOY)**(1.D0/3.D0)
  79 :       XCO=PRES*RMOY*XX/EPAIS/YOUNG
  80 :       XK=XK/(1.D0+XCO)
  81 :  30   CONTINUE
  82 :       IF(XK.LT.1.D00) XK=1.D00
  83 :  10   CONTINUE
  84 :       CARAC(2)=XIN/XK
  85 :       CARAC(3)=CARAC(2)
  86 :       RETURN
  87 :  999  KERRE=1
  88 :       RETURN
  89 :       END
  90 : 
  91 : 
  92 :  
  93 :  
  94 :