Télécharger la plaquette Cast3M


fr en

Télécharger initt.eso

Retour à la liste

Numérotation des lignes : 

initt

  1. C INITT     SOURCE    CHAT      05/01/13    00:38:23     5004
  2.       SUBROUTINE INITT(YOUNG,NYOUNG,XNU,NXNU,EN,NEN,EM,NEM,EK,NEK,ALPHA,
  3.      &NALPHA,BETA,NBETA,ER,NER,EA,NEA,ALPHAT,NALFAT,RHO,NRHO,SIGY,NSIGY,
  4.      &FKX,NKX,NPTEUR,NC,XMAT,NCOMAT,T0,TINF,TSUP,ITEST,TRUC,NCOURB)
  5. C
  6. C----------------------------------------------------------------
  7. C Objet: Initialiser les parametres scalaires caracteristiques du
  8. C        materiau a la temperature T0
  9. C----------------------------------------------------------------
  10. C
  11. C---------------------------------------------------------------------
  12. C Entree: YOUNG(NYOUNG) module d'Young en fonction de T
  13. C         XNU(NXNU) coefficient de Poisson en fonction de T
  14. C         SIGY(NSIGY) limite elastique en fonction de T
  15. C         RHO(NRHO) masse volumique en fonction de T
  16. C         ALPHAT(NALFAT) coeff de dilatation thermique en fonction de T
  17. C         EN(NEN)
  18. C         EM(NEM)
  19. C         EK(NEK)     courbes des autres coefficients caracteristiques
  20. C         ALPHA(NALPHA)  en fonction de T intervenant
  21. C         BETA(NBETA)   dans les lois d'evolution
  22. C         ER(NER)
  23. C         EA(NEA)
  24. C         FKX(NKX) fonction k(X) tabulee en fonction de la temperature
  25. C         NPTEUR(NC) tableau de pointeurs sur les courbes de k(X)
  26. C         T0 temperature d'interpolation
  27. C         ITEST = 1  listes de temperatures identiques pour coef. non lin.
  28. C               = 0  listes de temperatures differentes
  29. C----------------------------------------------------------------------
  30. C
  31. C--------------------------------------------------------------------
  32. C Sortie: XMAT(NCOMAT) tableau des parametres a la temperature T0
  33. C         TINF,TSUP bornes entre lesquelles est comprise T0 quand
  34. C         les listes de temperatures definissant la variation des
  35. C         coefficients non lineaires sont identiques
  36. C--------------------------------------------------------------------
  37.       IMPLICIT INTEGER(I-N)
  38.        IMPLICIT REAL*8(A-H,O-Z)
  39.        DIMENSION YOUNG(*),XNU(*),SIGY(*),RHO(*),ALPHAT(*),XMAT(*)
  40.        DIMENSION EN(*),EM(*),EK(*),ALPHA(*),BETA(*),ER(*),EA(*)
  41.        DIMENSION FKX(*),NPTEUR(*),TRUC(*)
  42. C
  43.        TINF=T0
  44.        TSUP=T0
  45.        CALL DERTRA(NYOUNG,YOUNG,T0,Y,YPRIM,TIN,TSU)
  46.        XMAT(1)=Y
  47.        CALL DERTRA(NXNU,XNU,T0,Y,YPRIM,TIN,TSU)
  48.        XMAT(2)=Y
  49.        CALL DERTRA(NEN,EN,T0,Y,YPRIM,TIN,TSU)
  50.        XMAT(3)=Y
  51.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  52.                   TINF=TIN
  53.                   TSUP=TSU
  54.                 ENDIF
  55.        CALL DERTRA(NEM,EM,T0,Y,YPRIM,TIN,TSU)
  56.        XMAT(4)=Y
  57.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  58.                   TINF=TIN
  59.                   TSUP=TSU
  60.                 ENDIF
  61.        CALL DERTRA(NEK,EK,T0,Y,YPRIM,TIN,TSU)
  62.        XMAT(5)=Y
  63.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  64.                   TINF=TIN
  65.                   TSUP=TSU
  66.                 ENDIF
  67.        CALL DERTRA(NALPHA,ALPHA,T0,Y,YPRIM,TIN,TSU)
  68.        XMAT(6)=Y
  69.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  70.                   TINF=TIN
  71.                   TSUP=TSU
  72.                 ENDIF
  73.        CALL DERTRA(NBETA,BETA,T0,Y,YPRIM,TIN,TSU)
  74.        XMAT(7)=Y
  75.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  76.                   TINF=TIN
  77.                   TSUP=TSU
  78.                 ENDIF
  79.        CALL DERTRA(NER,ER,T0,Y,YPRIM,TIN,TSU)
  80.        XMAT(8)=Y
  81.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  82.                   TINF=TIN
  83.                   TSUP=TSU
  84.                 ENDIF
  85.        CALL DERTRA(NEA,EA,T0,Y,YPRIM,TIN,TSU)
  86.        XMAT(9)=Y
  87.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  88.                   TINF=TIN
  89.                   TSUP=TSU
  90.                 ENDIF
  91.        CALL DERTRA(NALFAT,ALPHAT,T0,Y,YPRIM,TIN,TSU)
  92.        XMAT(10)=Y
  93.        CALL DERTRA(NRHO,RHO,T0,Y,YPRIM,TIN,TSU)
  94.        XMAT(12)=Y
  95.        CALL DERTRA(NSIGY,SIGY,T0,Y,YPRIM,TIN,TSU)
  96.        XMAT(13)=Y
  97.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  98.                   TINF=TIN
  99.                   TSUP=TSU
  100.                 ENDIF
  101.        X=0.
  102.        CALL VALPAR(FKX,NPTEUR,NKX,NC,X,Y,YPRIM,T0,TIN,TSU,TRUC,NCOURB)
  103.                 IF (ITEST.EQ.1.AND.TIN.NE.TSU) THEN
  104.                   TINF=TIN
  105.                   TSUP=TSU
  106.                 ENDIF
  107.        RETURN
  108.        END
  109.  
  110.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales