Sources Esope | Cast3M

tufifp
C TUFIFP    SOURCE    CHAT      05/01/13    03:55:41     5004
      SUBROUTINE TUFIFP(TESPI,AXX,FP,FM,FMP,FOP,FOM)
C=======================================================================
C   CALCULE LES FONCTIONS FP FM FMP    COEFF COMPLIANCE TUYAU FISSURE
C                         FOM FOP      FACTEUR DE FORME
C    ENTREES
C     TESPI = THETA SUR PI
C     AXX   = COEFFICIENT ZAHOR
C    SORTIES
C     FP FM FMP   =   COEFF COMPLIANCE
C     FOP FOM     =   FCTR FORME
C
C    EBERSOLT AVRIL 86
C=====================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION AL(19)
C
C   AL()=  COEFFICIENTS DES PUISSANCES DE TETA DANS FP FM FMP
C                               COEFFICIENTS ZAHOR
C      FP
      DATA AL(1 ),AL(2 )/6.09D0,12.03D0/
      DATA AL(3 ),AL(4 ),AL(5 )/11.36D0,51.70D0,67.24D0/
C      FM
      DATA AL(6 ),AL(7 )/5.26D0,1.69D0/
      DATA AL(8 ),AL(9 ),AL(10)/8.46D0,6.28D0,1.33D0/
C      FMP
      DATA AL(11),AL(12)/5.67D0,6.86D0/
      DATA AL(13),AL(14),AL(15)/9.81D0,25.95D0,9.46D0/
C    FM FP    FACT FORME
      DATA  AL(16),AL(17)/5.33D0,18.77D0/
      DATA  AL(18),AL(19)/4.60D0,2.64D0/
C
      TES1=TESPI**1.5D0
      TES2=TESPI**4.24D0
      TES3=TESPI**3
      TES4=TESPI**5.74D0
      TES5=TESPI**8.48D0
C
      FP = 1.D0 + AXX*(AL(1 )*TES1+AL(2 )*TES2)
     1     +AXX*AXX*(AL(3 )*TES3+AL(4 )*TES4+AL(5 )*TES5)
      FM = 1.D0 + AXX*(AL(6 )*TES1+AL(7 )*TES2)
     1     +AXX*AXX*(AL(8 )*TES3+AL(9 )*TES4+AL(10)*TES5)
      FMP= 1.D0 + AXX*(AL(11)*TES1+AL(12)*TES2)
     1     +AXX*AXX*(AL(13)*TES3+AL(14)*TES4+AL(15)*TES5)
C
      FOP = 1.D0 + AXX*(AL(16)*TES1+AL(17)*TES2)
      FOM = 1.D0 + AXX*(AL(18)*TES1+AL(19)*TES2)
      RETURN
      END