Télécharger matpo2.eso

Retour à la liste

Numérotation des lignes :

  1. C MATPO2 SOURCE CHAT 05/01/13 01:35:20 5004
  2. SUBROUTINE MATPO2 ( TAU, DT, TEMP, DTEMP, PREDEF, DPRED, WMAT,
  3. & IERUT )
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. C-----------------------------------------------------------------------
  7. C
  8. C DESCRIPTION FONCTIONNELLE :
  9. C -------------------------
  10. C
  11. C Dependance UMAPO2
  12. C => Dependance en deuxieme niveau du module Utilisateur UMAT
  13. C Appelee dans le cas suivant :
  14. C ILOI = 22 : Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21
  15. C RESTRICTION A LA FORMULATION MASSIVE
  16. C Instanciation des composantes de materiau au fur et a
  17. C mesure des iterations internes.
  18. C Le modele comprend cette fois des parametres externes.
  19. C Les composantes de materiau sont evaluees par les
  20. C memes fonctions externes que celles utilisees par
  21. C l'operateur VARI.
  22. C MODELE 'MAQUETTE'
  23. C
  24. C Calcul des proprietes de materiau au point TAU
  25. C
  26. C ENTREES : TAU, DT, TEMP, DTEMP, PREDEF, DPRED
  27. C SORTIES : WMAT, IERUT
  28. C
  29. C-----------------------------------------------------------------------
  30. C
  31. C Composantes de materiau
  32. C -----------------------
  33. C WMAT ( 1) : 'YOUN' WMAT ( 5) : 'SMAX'
  34. C WMAT ( 2) : 'NU ' WMAT ( 6) : 'AF0 '
  35. C WMAT ( 7) : 'AF1 '
  36. C WMAT ( 4) : 'ALPH' WMAT ( 8) : 'AF2 '
  37. C WMAT ( 9) : 'AF3 '
  38. C WMAT (10) : 'AF4 '
  39. C WMAT (11) : 'AF5 '
  40. C WMAT (12) : 'AF6 '
  41. C
  42. C N.B. WMAT(3) est laissee libre en coherence avec la declaration
  43. C des composantes de materiau du modele, correspond a 'RHO '
  44. C qui n'est pas utilisee par l'application TOUTATIS
  45. C
  46. C Parametres externes
  47. C -------------------
  48. C TEMP : 'T '
  49. C PREDEF( 1) : 'PORO'
  50. C PREDEF( 2) : 'YOGC'
  51. C PREDEF( 3) : 'ALPC'
  52. C PREDEF( 4) : 'TFIS'
  53. C PREDEF( 5) : 'TU02'
  54. C PREDEF( 6) : 'FACF'
  55. C PREDEF( 7) : 'DSIU'
  56. C PREDEF( 8) : 'DGRA'
  57. C
  58. C DTEMP : increment de 'T ' sur le pas de temps
  59. C DPRED ( 1) : increment de 'PORO' sur le pas de temps
  60. C DPRED ( 2) : increment de 'YOGC' sur le pas de temps - INUTILE
  61. C DPRED ( 3) : increment de 'ALPC' sur le pas de temps - INUTILE
  62. C DPRED ( 4) : increment de 'TFIS' sur le pas de temps
  63. C DPRED ( 5) : increment de 'TUO2' sur le pas de temps
  64. C DPRED ( 6) : increment de 'FACF' sur le pas de temps
  65. C DPRED ( 7) : increment de 'DSIU' sur le pas de temps
  66. C DPRED ( 8) : increment de 'DGRA' sur le pas de temps
  67. C
  68. C-----------------------------------------------------------------------
  69. C Arguments de l'interface
  70. C
  71. REAL*8 TAU, DT, TEMP, DTEMP, PREDEF(*), DPRED(*), WMAT(*)
  72. INTEGER IERUT
  73. C
  74. C Variables locales
  75. C
  76. INTEGER NBPAR
  77. REAL*8 VALPAR(5)
  78. REAL*8 DDT, TC, POROC, YOGCC, ALPCC, TFISC, TUO2C, FACFC,
  79. & DSIUC, DGRAC
  80. C
  81. C
  82. C------------------- Debut du code executable --------------------------
  83. C
  84. IERUT = 0
  85. C
  86. IF (INT(DPRED(2)).NE.0) THEN
  87. IERUT = -22002
  88. RETURN
  89. ENDIF
  90. IF (INT(DPRED(3)).NE.0) THEN
  91. IERUT = -22003
  92. RETURN
  93. ENDIF
  94. C
  95. C 1 - Evaluation des parametres au point TAU
  96. C
  97. DDT = TAU/DT
  98. C
  99. TC = TEMP + ( DDT * DTEMP )
  100. POROC = PREDEF(1) + ( DDT * DPRED(1) )
  101. YOGCC = PREDEF(2)
  102. ALPCC = PREDEF(3)
  103. TFISC = PREDEF(4) + ( DDT * DPRED(4) )
  104. TUO2C = PREDEF(5) + ( DDT * DPRED(5) )
  105. FACFC = PREDEF(6) + ( DDT * DPRED(6) )
  106. DSIUC = PREDEF(7) + ( DDT * DPRED(7) )
  107. DGRAC = PREDEF(8) + ( DDT * DPRED(8) )
  108. C
  109. C 2 - Evaluation des proprietes de materiau au point TAU
  110. C
  111. C.....Module d'Young 'YOUN'
  112. C
  113. NBPAR = 3
  114. VALPAR(1) = TC
  115. VALPAR(2) = POROC
  116. VALPAR(3) = YOGCC
  117. WMAT(1) = FVYOUN(VALPAR,NBPAR,IERUT)
  118. IF (IERUT.GT.0) THEN
  119. IERUT = -22000-IERUT
  120. RETURN
  121. ENDIF
  122. C
  123. C.....Coefficient de Poisson 'NU '
  124. C
  125. WMAT(2) = FVNU(VALPAR,NBPAR,IERUT)
  126. IF (IERUT.GT.0) THEN
  127. IERUT = -22000-IERUT
  128. RETURN
  129. ENDIF
  130. C
  131. C.....Modele de fluage : contrainte de reference 'SMAX'
  132. C
  133. WMAT(5) = FVSMAX(VALPAR,NBPAR,IERUT)
  134. IF (IERUT.GT.0) THEN
  135. IERUT = -22000-IERUT
  136. RETURN
  137. ENDIF
  138. C
  139. C.....Coefficient de dilatation thermique 'ALPH'
  140. C
  141. NBPAR = 2
  142. VALPAR(2) = ALPCC
  143. WMAT(4) = FVALPH(VALPAR,NBPAR,IERUT)
  144. IF (IERUT.GT.0) THEN
  145. IERUT = -22000-IERUT
  146. RETURN
  147. ENDIF
  148. C
  149. C.....Modele de fluage : coefficient du polynome 'AF0 '
  150. C
  151. NBPAR = 3
  152. VALPAR(1) = TFISC
  153. VALPAR(2) = TUO2C
  154. VALPAR(3) = FACFC
  155. WMAT(6) = FVAF0(VALPAR,NBPAR,IERUT)
  156. IF (IERUT.GT.0) THEN
  157. IERUT = -22000-IERUT
  158. RETURN
  159. ENDIF
  160. C
  161. C.....Modele de fluage : coefficient du polynome 'AF2 '
  162. C
  163. WMAT(8) = FVAF2(VALPAR,NBPAR,IERUT)
  164. IF (IERUT.GT.0) THEN
  165. IERUT = -22000-IERUT
  166. RETURN
  167. ENDIF
  168. C
  169. C.....Modele de fluage : coefficient du polynome 'AF4 '
  170. C
  171. WMAT(10) = FVAF4(VALPAR,NBPAR,IERUT)
  172. IF (IERUT.GT.0) THEN
  173. IERUT = -22000-IERUT
  174. RETURN
  175. ENDIF
  176. C
  177. C.....Modele de fluage : coefficient du polynome 'AF5 '
  178. C
  179. WMAT(11) = FVAF5(VALPAR,NBPAR,IERUT)
  180. IF (IERUT.GT.0) THEN
  181. IERUT = -22000-IERUT
  182. RETURN
  183. ENDIF
  184. C
  185. C.....Modele de fluage : coefficient du polynome 'AF6 '
  186. C
  187. WMAT(12) = FVAF6(VALPAR,NBPAR,IERUT)
  188. IF (IERUT.GT.0) THEN
  189. IERUT = -22000-IERUT
  190. RETURN
  191. ENDIF
  192. C
  193. C.....Modele de fluage : coefficient du polynome 'AF3 '
  194. C
  195. NBPAR = 4
  196. VALPAR(4) = DSIUC
  197. WMAT(9) = FVAF3(VALPAR,NBPAR,IERUT)
  198. IF (IERUT.GT.0) THEN
  199. IERUT = -22000-IERUT
  200. RETURN
  201. ENDIF
  202. C
  203. C.....Modele de fluage : coefficient du polynome 'AF1 '
  204. C
  205. NBPAR = 5
  206. VALPAR(5) = DGRAC
  207. WMAT(7) = FVAF1(VALPAR,NBPAR,IERUT)
  208. IF (IERUT.GT.0) THEN
  209. IERUT = -22000-IERUT
  210. RETURN
  211. ENDIF
  212. C
  213. RETURN
  214. END
  215.  
  216.  
  217.  

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