dilth
C DILTH SOURCE STRU 08/12/22 21:15:09 6246 c Sous-Programme du module MISTRAL0 C -------------------------------------------------------------------------- C -------------------------------------------------------------------------- C Calcul des vitesses de dilatation thermique VEPSTH(IJ) C en fonction de : C - la temperature TT, C - sa derivee par rapport au temps TTP, C - et des parametres transmis dans PDILT, C pour un materiau orthotrope, dans la base principale d'orthotropie. C -------------------------------------------------------------------------- IMPLICIT INTEGER (I-N) IMPLICIT REAL*8 (A-H, O-Z) DIMENSION PDILT(1:*) DIMENSION VEPSTH(1:*) ISOTH = NINT(PDILT(1)) IP = 1 IF (ISOTH.EQ.1) THEN C ----------------- C Materiau isotrope C DO IJ = 1,3 END DO ELSE IF (ISOTH.EQ.0) THEN C ------------------- C Materiau orthotrope C DO IJ = 1,3 VEPSTH(IJ) = (PDILT(IP+1)+(PDILT(IP+2)+PDILT(IP+3)*TT)*TT)*TTP IP = IP+3 END DO ELSE IF (ISOTH.EQ.-1) THEN C --------------------------- C Pas de dilatation thermique C DO IJ = 1,3 VEPSTH(IJ) = 0. END DO END IF C ------ DO IJ = 4,6 VEPSTH(IJ) = 0. END DO RETURN END
© Cast3M 2003 - Tous droits réservés.
Mentions légales