Esope source files

   1 : C HILL      SOURCE    KK2000    05/05/25    21:15:17     5088           
   2 : c     Sous-Programme du module MISTRAL0
   3 : C     --------------------------------------------------------------------------
   4 :       SUBROUTINE HILL (TT,PSI,S, ID, PCOHI, NP, SEQ,DSESDS,
   5 :      &                 DSEDTT,DSDPSI)
   6 : C     --------------------------------------------------------------------------
   7 : C     Calcul de :
   8 : C     - la contrainte equivalente SEQ,
   9 : C     - ses derivees partielles DSESDS(IJ) par rapport aux S(IJ),
  10 : C     puis, si ID = 1, DSEDTT et DSDPSI par rapport a TT et PSI,
  11 : C     en fonction :
  12 : C     - de la temperature TT,
  13 : C     - de la variable de durcissement du a l'irradiation PSI,
  14 : C     - des contraintes S(IJ),
  15 : C     - des parametres transmis dans PCOHI,
  16 : C     NP (incremente dans COHILL) etant en entree le nombre
  17 : C     de ceux correspondant aux deformations de natures IDPP < IDP.
  18 : C     --------------------------------------------------------------------------
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8 (A-H, O-Z)
  21 :       DIMENSION S(1:*)
  22 :       DIMENSION PCOHI(1:*)
  23 :       DIMENSION H(1:6),DHDTT(1:6),DHDPSI(1:6)
  24 :         DIMENSION ISU(1:3)
  25 :       DIMENSION DSESDS(1:*)
  26 : C     ------------------------------------------------------------------
  26 :       CALL COHILL (TT,PSI, ID, PCOHI, NP, H,DHDTT,DHDPSI)
  27 :       SEQ = SQRT(
  28 :      &      H(1)*(S(2)-S(3))**2+H(2)*(S(3)-S(1))**2+H(3)*(S(1)-S(2))**2
  29 :      &     +2.*H(4)*S(4)**2+2.*H(5)*S(5)**2+2.*H(6)*S(6)**2 )
  30 :         IF (SEQ.NE.0.) THEN
  30 :          CALL PROMAH (H,S, DSESDS)
  31 :        DO I = 1,6
  32 :           DSESDS(I) = DSESDS(I)/SEQ
  33 :          END DO
  34 :         END IF
  35 : C     ------------------------------------------------------------------
  36 :       IF (ID.EQ.1) THEN
  37 :        ISU(1) = 2 
  38 :        ISU(2) = 3 
  39 :        ISU(3) = 1
  40 :        DSEDTT = 0.
  41 :        DSDPSI = 0.
  42 :        DO I = 1,6
  43 :         IF ((DHDTT(I).NE.0.D0).OR.(DHDPSI(I).NE.0.D0)) THEN
  44 :          IF (I.LE.3) THEN
  45 :           J = ISU(I)
  46 :             K = ISU(J)
  47 :             SS2 = (S(J)-S(K))**2
  48 :            ELSE
  49 :             SS2 = 2.*S(I)**2
  50 :          END IF
  51 :            DSEDTT = DSEDTT+SS2*DHDTT(I)
  52 :            DSDPSI = DSDPSI+SS2*DHDPSI(I)
  53 :         END IF
  54 :          END DO
  55 :        DSEDTT = DSEDTT/2./SEQ
  56 :        DSDPSI = DSDPSI/2./SEQ
  57 :       END IF
  58 : C     ------------------------------------------------------------------
  59 :         RETURN
  60 :         END
  61 :  
  62 :  
  63 :