Télécharger la plaquette Cast3M


fr en

Télécharger comcri.eso

Retour à la liste

Numérotation des lignes : 

comcri

  1. C COMCRI    SOURCE    CB215821  24/04/12    21:15:20     11897          
  2.  
  3.       SUBROUTINE COMCRI(iqmod,ipcon,IPMINT,indeso,ipil,
  4.      &                  insupp,iwrk53,iretou)
  5. *
  6. *  cree les melval produits par le modele
  7. *  cree les deche associes
  8. *  range dans les pilnec qui conviennent
  9. *----------------------------------------------------------------
  10.       IMPLICIT INTEGER(I-N)
  11.       IMPLICIT REAL*8(A-H,O-Z)
  12.  
  13. -INC PPARAM
  14. -INC CCOPTIO
  15. -INC SMMODEL
  16. -INC SMCHAML
  17. -INC DECHE
  18.  
  19.       imodel = iqmod
  20.       liluc  = ipil
  21.  
  22.       if (formod(1).EQ.'METALLURGIE     ') then
  23. * variables internes -> 20 (Proportions des phases)
  24.          jluc = 20
  25.          mocomp = liluc(jluc,1)
  26.          lscont = liluc(jluc,2)
  27.          call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  28.      &               IPMINT,insupp,iwrk53,jluc,iretou)
  29.          return
  30.       endif
  31.  
  32.  
  33. * caracteristiques materiau -> 13
  34.          jluc = 13
  35.          mocomp = liluc(jluc,1)
  36.          lscont = liluc(jluc,2)
  37.          call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  38.      &               IPMINT,insupp,iwrk53,jluc,iretou)
  39. * caracteristiques geometriques -> 14
  40.          jluc = 14
  41.          mocomp = liluc(jluc,1)
  42.          lscont = liluc(jluc,2)
  43.          call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  44.      &               IPMINT,insupp,iwrk53,jluc,iretou)
  45.  
  46. * temperature -> 15
  47. *        jluc = 15
  48. *        mocomp = jnomid
  49. *        lscont = jilnec
  50. *        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  51. *     &              IPMINT,insupp,iwrk53,jluc,iretou)
  52. * -> 26 et + : Cas pouvant arriver ?
  53.          IF(LILUC(/1).GT.25) THEN
  54.            jluc = 26
  55.            mocomp = liluc(jluc,1)
  56.            lscont = liluc(jluc,2)
  57.            call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  58.      &                 IPMINT,insupp,iwrk53,jluc,iretou)
  59.          ENDIF
  60.  
  61. *-------------------------------------------------
  62.       if(FORMOD(1).EQ.'MECANIQUE       ' .OR. 
  63.      &   FORMOD(1).EQ.'POREUX          ') then
  64.  
  65. * contraintes -> 11
  66.        jluc = 11
  67.        mocomp = liluc(jluc,1)
  68.        lscont = liluc(jluc,2)
  69.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  70.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  71. * deformations -> 12
  72.        jluc = 12
  73.        mocomp = liluc(jluc,1)
  74.        lscont = liluc(jluc,2)
  75.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  76.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  77. * variables internes -> 20
  78.        jluc = 20
  79.        mocomp = liluc(jluc,1)
  80.        lscont = liluc(jluc,2)
  81.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  82.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  83. * deformations inelastiques -> 24
  84.        jluc = 24
  85.        mocomp = liluc(jluc,1)
  86.        lscont = liluc(jluc,2)
  87.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  88.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  89. *--------------------------------------------
  90.       ELSE if (formod(1).eq.'LIAISON         ') then
  91. * variables internes -> 20
  92.        jluc = 20
  93.        mocomp = liluc(jluc,1)
  94.        lscont = liluc(jluc,2)
  95.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  96.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  97. *--------------------------------------------
  98.       ELSE if (formod(1).eq.'MELANGE         ') then
  99. * phases -> 23
  100.          jluc = 23
  101.          mocomp = liluc(jluc,1)
  102.          lscont = liluc(jluc,2)
  103.          call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  104.      &               IPMINT,insupp,iwrk53,jluc,iretou)
  105. *--------------------------------------------
  106.       ELSE if (formod(1).eq.'DIFFUSION       ') then
  107. * courant
  108.        jluc = 11
  109.        mocomp = liluc(jluc,1)
  110.        lscont = liluc(jluc,2)
  111.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  112.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  113. * gradient
  114.        jluc = 12
  115.        mocomp = liluc(jluc,1)
  116.        lscont = liluc(jluc,2)
  117.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  118.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  119. * force
  120.        jluc = 7
  121.        mocomp = liluc(jluc,1)
  122.        lscont = liluc(jluc,2)
  123.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  124.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  125.  
  126. * variables internes -> 20
  127.        jluc = 20
  128.        mocomp = liluc(jluc,1)
  129.        lscont = liluc(jluc,2)
  130.        call COMCRE(iqmod,ipcon,indeso,mocomp,lscont,
  131.      &             IPMINT,insupp,iwrk53,jluc,iretou)
  132. *--------------------------------------------
  133.       ENDIF
  134. *--------------------------------------------
  135.  
  136.       return
  137.       END
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  

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