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

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