Télécharger la plaquette Cast3M


fr en

Télécharger model6.eso

Retour à la liste

Numérotation des lignes : 

model6

  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