Télécharger model6.eso

Retour à la liste

Numérotation des lignes :

  1. C MODEL6 SOURCE CHAT 05/01/13 01:48:31 5004
  2. SUBROUTINE MODEL6(LESPRO,NBPRO,LESTEF,NBTEF,N1MAX,N2MAX)
  3. *----------------------------------------------------------------------
  4. * OBJET MODELE - CAS DE LA FORMULATION 'POREUX'
  5. * DEFINITION DES NOMS DE MODELES DE MATERIAU
  6. * ET DES NOMS DE TYPES D'ELEMENTS FINIS
  7. *----------------------------------------------------------------------
  8. IMPLICIT INTEGER(I-N)
  9. -INC CCOPTIO
  10. *
  11. CHARACTER*16 MOMODL(100)
  12. *
  13. CHARACTER*16 LESPRO(*)
  14. CHARACTER*4 LESTEF(*)
  15. *
  16. * Noms des comportements lineaires
  17. *
  18. CALL MODLIN (MOMODL,NMOD)
  19. *
  20. DO 1 I=1,NMOD
  21. LESPRO(I)=MOMODL(I)
  22. 1 CONTINUE
  23. NBPRO=NMOD
  24. *
  25. * Noms des materiaux pour le comportement lineaire
  26. *
  27. CALL MODELA (MOMODL,NMOD)
  28. *
  29. DO 2 I=NBPRO+1,NBPRO+NMOD
  30. LESPRO(I)=MOMODL(I-NBPRO)
  31. 2 CONTINUE
  32. NBPRO=NBPRO+NMOD
  33. *
  34. * Noms des comportements non lineaires
  35. *
  36. * Remarque
  37. * La routine MODNLI va retourner les noms des sept familles de
  38. * comportements non lineaires, dans l'ordre suivant :
  39. * 'PLASTIQUE', 'FLUAGE', 'VISCOPLASTIQUE', 'ENDOMMAGEMENT',
  40. * 'PLASTIQUE_ENDOM', 'NON_LINEAIRE' et 'VISCO_EXTERNE'.
  41. * Les noms de toutes les familles de comportements non lineaires
  42. * feront donc partie du vocabulaire admissible dans la definition
  43. * du modele, sans restriction.
  44. * Toutefois, DANS LE CAS DE LA FORMULATION 'POREUX', seuls les
  45. * materiaux des 6 PREMIERES FAMILLES seront acceptables : cf. appels
  46. * ci-dessous aux routines MODPLA, MODFLU, MODVIS et MODEND.
  47. * Il y a donc une restriction en formulation 'POREUX' par rapport
  48. * a la formulation 'MECANIQUE', mais on ne la fait pas apparaitre
  49. * au stade de la definition des familles admissibles.
  50.  
  51. CALL MODNLI (MOMODL,NMOD)
  52. *
  53. DO 3 I=NBPRO+1,NBPRO+NMOD
  54. LESPRO(I)=MOMODL(I-NBPRO)
  55. 3 CONTINUE
  56. NBPRO=NBPRO+NMOD
  57.  
  58. *
  59. * Noms des materiaux plastiques
  60. *
  61. CALL MODPLA (MOMODL,NMOD)
  62. *
  63. DO 4 I=NBPRO+1,NBPRO+NMOD
  64. LESPRO(I)=MOMODL(I-NBPRO)
  65. 4 CONTINUE
  66. NBPRO=NBPRO+NMOD
  67. *
  68. * Noms des materiaux fluages
  69. *
  70. CALL MODFLU (MOMODL,NMOD)
  71. *
  72. DO 5 I=NBPRO+1,NBPRO+NMOD
  73. LESPRO(I)=MOMODL(I-NBPRO)
  74. 5 CONTINUE
  75. NBPRO=NBPRO+NMOD
  76. *
  77. * Noms des materiaux viscoplastiques
  78. *
  79. CALL MODVIS (MOMODL,NMOD)
  80. *
  81. DO 6 I=NBPRO+1,NBPRO+NMOD
  82. LESPRO(I)=MOMODL(I-NBPRO)
  83. 6 CONTINUE
  84. NBPRO=NBPRO+NMOD
  85. *
  86. * Noms des materiaux endommageables
  87. *
  88. CALL MODEND (MOMODL,NMOD)
  89. *
  90. DO 7 I=NBPRO+1,NBPRO+NMOD
  91. LESPRO(I)=MOMODL(I-NBPRO)
  92. 7 CONTINUE
  93. NBPRO=NBPRO+NMOD
  94. *
  95. * Noms des materiaux plastiques endommageables
  96. *
  97. CALL MODPLE (MOMODL,NMOD)
  98. *
  99. DO 8 I=NBPRO+1,NBPRO+NMOD
  100. LESPRO(I)=MOMODL(I-NBPRO)
  101. 8 CONTINUE
  102. NBPRO=NBPRO+NMOD
  103. *
  104. * Noms des materiaux elastiques non linéaires
  105. *
  106. CALL MODENL (MOMODL,NMOD)
  107. *
  108. DO 9 I=NBPRO+1,NBPRO+NMOD
  109. LESPRO(I)=MOMODL(I-NBPRO)
  110. 9 CONTINUE
  111. NBPRO=NBPRO+NMOD
  112. *
  113. IF(NBPRO.GT.N1MAX) THEN
  114. CALL ERREUR(550)
  115. GO TO 99
  116. END IF
  117. *
  118. * NOMS DE TYPES D'ELEMENTS FINIS
  119. *
  120. IF(IDIM.EQ.2) THEN
  121. NBTEF=9
  122. IF(NBTEF.GT.N2MAX) THEN
  123. CALL ERREUR(551)
  124. GO TO 99
  125. END IF
  126. LESTEF(1)='TRIP'
  127. LESTEF(2)='QUAP'
  128. LESTEF(3)='JOP3'
  129. LESTEF(4)='TRIQ'
  130. LESTEF(5)='QUAQ'
  131. LESTEF(6)='JOQ3'
  132. LESTEF(7)='TRIR'
  133. LESTEF(8)='QUAR'
  134. LESTEF(9)='JOR3'
  135. ELSE
  136. NBTEF=15
  137. IF(NBTEF.GT.N2MAX) THEN
  138. CALL ERREUR(551)
  139. GO TO 99
  140. END IF
  141. LESTEF(1)='CUBP'
  142. LESTEF(2)='TETP'
  143. LESTEF(3)='PRIP'
  144. LESTEF(4)='JOP6'
  145. LESTEF(5)='JOP8'
  146. LESTEF(6)='CUBQ'
  147. LESTEF(7)='TETQ'
  148. LESTEF(8)='PRIQ'
  149. LESTEF(9)='JOQ6'
  150. LESTEF(10)='JOQ8'
  151. LESTEF(11)='CUBR'
  152. LESTEF(12)='TETR'
  153. LESTEF(13)='PRIR'
  154. LESTEF(14)='JOR6'
  155. LESTEF(15)='JOR8'
  156. ENDIF
  157. RETURN
  158. 99 CALL ERREUR(5)
  159. RETURN
  160. END
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  

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