Esope source files

   1 : C TUFIFP    SOURCE    CHAT      05/01/13    03:55:41     5004           
   2 :       SUBROUTINE TUFIFP(TESPI,AXX,FP,FM,FMP,FOP,FOM)
   3 : C=======================================================================
   4 : C   CALCULE LES FONCTIONS FP FM FMP    COEFF COMPLIANCE TUYAU FISSURE
   5 : C                         FOM FOP      FACTEUR DE FORME
   6 : C    ENTREES
   7 : C     TESPI = THETA SUR PI
   8 : C     AXX   = COEFFICIENT ZAHOR
   9 : C    SORTIES
  10 : C     FP FM FMP   =   COEFF COMPLIANCE
  11 : C     FOP FOM     =   FCTR FORME
  12 : C     
  13 : C    EBERSOLT AVRIL 86
  14 : C=====================================================================
  15 :       IMPLICIT INTEGER(I-N)
  16 :       IMPLICIT REAL*8(A-H,O-Z)
  17 :       DIMENSION AL(19)
  18 : C     
  19 : C   AL()=  COEFFICIENTS DES PUISSANCES DE TETA DANS FP FM FMP
  20 : C                               COEFFICIENTS ZAHOR
  21 : C      FP
  22 :       DATA AL(1 ),AL(2 )/6.09D0,12.03D0/
  23 :       DATA AL(3 ),AL(4 ),AL(5 )/11.36D0,51.70D0,67.24D0/
  24 : C      FM
  25 :       DATA AL(6 ),AL(7 )/5.26D0,1.69D0/
  26 :       DATA AL(8 ),AL(9 ),AL(10)/8.46D0,6.28D0,1.33D0/
  27 : C      FMP
  28 :       DATA AL(11),AL(12)/5.67D0,6.86D0/
  29 :       DATA AL(13),AL(14),AL(15)/9.81D0,25.95D0,9.46D0/
  30 : C    FM FP    FACT FORME
  31 :       DATA  AL(16),AL(17)/5.33D0,18.77D0/
  32 :       DATA  AL(18),AL(19)/4.60D0,2.64D0/
  33 : C     
  34 :       TES1=TESPI**1.5D0
  35 :       TES2=TESPI**4.24D0
  36 :       TES3=TESPI**3
  37 :       TES4=TESPI**5.74D0
  38 :       TES5=TESPI**8.48D0
  39 : C     
  40 :       FP = 1.D0 + AXX*(AL(1 )*TES1+AL(2 )*TES2)
  41 :      1     +AXX*AXX*(AL(3 )*TES3+AL(4 )*TES4+AL(5 )*TES5)
  42 :       FM = 1.D0 + AXX*(AL(6 )*TES1+AL(7 )*TES2)
  43 :      1     +AXX*AXX*(AL(8 )*TES3+AL(9 )*TES4+AL(10)*TES5)
  44 :       FMP= 1.D0 + AXX*(AL(11)*TES1+AL(12)*TES2)
  45 :      1     +AXX*AXX*(AL(13)*TES3+AL(14)*TES4+AL(15)*TES5)
  46 : C     
  47 :       FOP = 1.D0 + AXX*(AL(16)*TES1+AL(17)*TES2)
  48 :       FOM = 1.D0 + AXX*(AL(18)*TES1+AL(19)*TES2)
  49 :       RETURN
  50 :       END
  51 : 
  52 : 
  53 :  
  54 :