Esope source files

   1 : C DILTH     SOURCE    STRU      08/12/22    21:15:09     6246           
   2 : c     Sous-Programme du module MISTRAL0
   3 : C     --------------------------------------------------------------------------   
   4 :       SUBROUTINE DILTH (TT,TTP, PDILT, VEPSTH)           
   5 : C     --------------------------------------------------------------------------
   6 : C     Calcul des vitesses de dilatation thermique VEPSTH(IJ)
   7 : C     en fonction de :
   8 : C     - la temperature TT,
   9 : C     - sa derivee par rapport au temps TTP,
  10 : C     - et des parametres transmis dans PDILT,
  11 : C     pour un materiau orthotrope, dans la base principale d'orthotropie.
  12 : C     --------------------------------------------------------------------------
  13 :       IMPLICIT INTEGER (I-N)
  14 :       IMPLICIT REAL*8 (A-H, O-Z)
  15 :       DIMENSION PDILT(1:*)
  16 :         DIMENSION VEPSTH(1:*)
  17 :       ISOTH = NINT(PDILT(1))
  18 :       IP = 1
  19 :         IF (ISOTH.EQ.1) THEN
  20 : C     -----------------
  21 : C     Materiau isotrope
  22 : C       
  23 :          ALPHA = PDILT(2)+(PDILT(3)+PDILT(4)*TT)*TT
  24 :          DO IJ = 1,3
  25 :           VEPSTH(IJ) = ALPHA*TTP
  26 :          END DO
  27 :         ELSE IF (ISOTH.EQ.0) THEN
  28 : C     -------------------
  29 : C     Materiau orthotrope
  30 : C       
  31 :          DO IJ = 1,3
  32 :           VEPSTH(IJ) = (PDILT(IP+1)+(PDILT(IP+2)+PDILT(IP+3)*TT)*TT)*TTP
  33 :         IP = IP+3
  34 :          END DO
  35 :         ELSE IF (ISOTH.EQ.-1) THEN
  36 : C     ---------------------------
  37 : C     Pas de dilatation thermique
  38 : C       
  39 :          DO IJ = 1,3
  40 :           VEPSTH(IJ) = 0.
  41 :          END DO
  42 :         END IF
  43 : C     ------
  44 :         DO IJ = 4,6
  45 :          VEPSTH(IJ) = 0.
  46 :         END DO
  47 :         RETURN
  48 :         END
  49 :