Esope source files

   1 : C COHILL    SOURCE    STRU      08/12/22    21:15:05     6246           
   2 : c     Sous-Programme du module MISTRAL0
   3 : C     --------------------------------------------------------------------------
   4 :       SUBROUTINE COHILL (TT,PSI, ID, PCOHI, NP, H,DHDTT,DHDPSI)
   5 : C     --------------------------------------------------------------------------
   6 : C     Calcul des coefficients de Hill H(I),
   7 : C     et de leur derivees partielles DHDTT(I) et DHDPSI(I) si ID = 1,
   8 : C     en fonction de :
   9 : C     - la temperature TT,
  10 : C     - la variable de durcissement du a l'irradiation PSI,
  11 : C     - des parametres transmis dans PCOHI,
  12 : C     NP (incremente ici) etant en entree le nombre
  13 : C     de ceux correspondant aux deformations de natures IDPP < IDP.
  14 : C     --------------------------------------------------------------------------
  15 :       IMPLICIT INTEGER (I-N)
  16 :       IMPLICIT REAL*8 (A-H, O-Z)
  17 :       DIMENSION PCOHI(1:*)
  18 : c     (NPCHM = 6*10, NPCOHM = NPCHM*NDP)
  19 :       DIMENSION H(1:*),DHDTT(1:*),DHDPSI(1:*)
  20 :       DO I = 1,6
  21 :        NP = NP+1
  22 :        KH = NINT(PCOHI(NP))
  23 :        IF (KH.EQ.1) THEN
  24 :         NP = NP+1
  25 :         H(I) = PCOHI(NP)
  26 :         IF (ID.EQ.1) THEN
  27 :          DHDTT(I) = 0.D0
  28 :          DHDPSI(I) = 0.D0
  29 :         END IF
  30 :        ELSE IF (KH.EQ.2) THEN
  31 :         H(I) = PCOHI(NP+1)+(PCOHI(NP+2)+PCOHI(NP+3)*TT)*TT
  32 :         IF (ID.EQ.1) THEN
  33 :          DHDTT(I) = PCOHI(NP+2)+2.*PCOHI(NP+3)*TT
  34 :          DHDPSI(I) = 0.D0
  35 :         END IF
  36 :         NP = NP+3
  37 :        ELSE IF (KH.EQ.3) THEN
  38 :         HI0 = PCOHI(NP+1)
  39 :         HI1 = PCOHI(NP+2)
  40 :         HI2 = PCOHI(NP+3)
  41 :         HI3 = PCOHI(NP+4)
  42 :         HI4 = PCOHI(NP+5)
  43 :         HI5 = PCOHI(NP+6)
  44 :         HI6 = PCOHI(NP+7)
  45 :         HI7 = PCOHI(NP+8)
  46 :         HI8 = PCOHI(NP+9)
  47 :         EMHI7F = EXP(-HI7*PSI)
  48 :         EMHI8F = EXP(-HI8*PSI)
  49 :         HIUN = HI1+(HI2+HI3*TT)*TT
  50 :         HIDE = HI4+(HI5+HI6*TT)*TT
  51 :           H(I) = HI0+HIUN*EMHI7F+HIDE*(1.-EMHI8F)
  52 :         IF (ID.EQ.1) THEN
  53 :          DHDTT(I) = (HI2+2.*HI3*TT)*EMHI7F+(HI5+2.*HI6*TT)*(1.-EMHI8F)
  54 :          DHDPSI(I) = -HI7*HIUN*EMHI7F+HI8*HIDE*EMHI8F
  55 :         END IF
  56 :         NP = NP+9
  57 :        ELSE IF (KH.EQ.4) THEN
  58 :         HI0 = PCOHI(NP+1)
  59 :         HI1 = PCOHI(NP+2)
  60 :         HI2 = PCOHI(NP+3)
  61 :         HI3 = PCOHI(NP+4)
  62 :         HI4 = PCOHI(NP+5)
  63 :         HI5 = PCOHI(NP+6)
  64 :         HI6 = PCOHI(NP+7)
  65 :         HI7 = PCOHI(NP+8)
  66 :         HI8 = PCOHI(NP+9)
  67 :         EMHPSI = EXP(-HI7*PSI**HI8)
  68 :         HIUN = HI1+(HI2+HI3*TT)*TT
  69 :         HIDE = HI4+(HI5+HI6*TT)*TT
  70 :           H(I) = HI0+HIUN*EMHPSI+HIDE*(1.-EMHPSI)
  71 :         IF (ID.EQ.1) THEN
  72 :          DHDTT(I) = (HI2+2.*HI3*TT)*EMHPSI+(HI5+2.*HI6*TT)*(1.-EMHPSI)
  73 :          DHDPSI(I) = -(HIUN-HIDE)*HI7*HI8*PSI**(HI8-1.)*EMHPSI
  74 :         END IF
  75 :         NP = NP+9
  76 :        END IF
  77 :       END DO
  78 :         RETURN
  79 :         END
  80 :