Télécharger initt.eso

Retour à la liste

Numérotation des lignes :

  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