Sources Esope | Cast3M

model2
C MODEL2    SOURCE    JB251061  22/12/19    21:15:05     11522          
 
C=======================================================================
C=  Objet MODELE  -  Cas de la formulation 'MECANIQUE'                 =
C=                                                                     =
C=   Definition des noms de modeles de MATERIAU                        =
C=           et des noms de types d'ELEMENTS FINIS                     =
C=======================================================================
 
      SUBROUTINE MODEL2 (LESPRO,NBPRO,LESTEF,NBTEF,N1MAX,N2MAX)
 
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
 
 
-INC PPARAM
-INC CCOPTIO
 
      CHARACTER*16 MOMODL(200)
      CHARACTER*16 LESPRO(*)
      CHARACTER*4 LESTEF(*)
 
C  1 - Noms des comportements LINEAIRES
C --------------------------------------
      CALL MODLIN(MOMODL,NMOD)
      DO i=1,NMOD
        LESPRO(i)=MOMODL(i)
      ENDDO
      NBPRO=NMOD
 
C  2 - Noms des MATERIAUx pour le comportement LINEAIRE
C ------------------------------------------------------
      CALL MODELA(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  3 - Noms des comportements NON LINEAIRES
C ------------------------------------------
      CALL MODNLI(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  4 - Noms des modeles PLASTIQUES
C ---------------------------------
      CALL MODPLA(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  5 - Noms des modeles de FLUAGE
C --------------------------------
      CALL MODFLU(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  6 - Noms des modeles VISCOPLASTIQUES
C --------------------------------------
      CALL MODVIS(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  7 - Noms des modeles ENDOMMAGEABLES
C -------------------------------------
      CALL MODEND(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  8 - Noms des modeles PLASTIQUES ENDOMMAGEABLES
C ------------------------------------------------
      CALL MODPLE(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  9 - Noms des modeles ELASTIQUES NON LINEAIRES
C -----------------------------------------------
      CALL MODENL(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  10 - Noms des modeles VISCOPLASTIQUES EXTERNES
C ------------------------------------------------
      CALL MODVIX(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  12 - Noms des modeles LIAISONS
C ------------------------------------------------
      CALL MODLIA(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  13 - Noms des modeles IMPEDANCE
C ------------------------------------------------
      CALL MODIMP(MOMODL,NMOD)
      DO i=NBPRO+1,NBPRO+NMOD
        LESPRO(i)=MOMODL(i-NBPRO)
      ENDDO
      NBPRO=NBPRO+NMOD
 
C  restauration du comportement si FUSION
      LESPRO(NBPRO + 1) = 'FUSION'
      NBPRO = NBPRO + 1
 
C  CRItere Pression
      LESPRO(NBPRO + 1) = 'CRIP'
      NBPRO = NBPRO + 1
 
      IF (NBPRO.GT.N1MAX) GOTO 99
 
C  11 - Noms de types d'ELEMENTS FINIS
C -------------------------------------
C  DIMENSION 1 :
C ---------------
      IF (IDIM.EQ.1) THEN
        NBTEF=2
        IF (NBTEF.GT.N2MAX) GOTO 99
        LESTEF(1)='M1D2'
        LESTEF(2)='M1D3'
C  DIMENSION 2 :
C ---------------
C  Les EF sont definis suivant la valeur de IFOUR
      ELSE IF (IDIM.EQ.2) THEN
        IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
          NBTEF=19
          IF (NBTEF.GT.N2MAX) GOTO 99
          LESTEF(1)='TRI3'
          LESTEF(2)='TRI6'
          LESTEF(3)='QUA4'
          LESTEF(4)='QUA8'
          LESTEF(5)='COQ2'
          LESTEF(6)='BARR'
          LESTEF(7)='TRIH'
          LESTEF(8)='CERC'
          LESTEF(9)='JOI2'
          LESTEF(10)='JOI3'
          LESTEF(11)='POLY'
          LESTEF(12)='BAR3'
          LESTEF(13)='QUAH'
          LESTEF(14)='Q4RI'
          LESTEF(15)='Q8RI'
          LESTEF(16)='XQ4R'
          LESTEF(17)='ZCO2'
          LESTEF(18)='COS2'
          LESTEF(19)='COA2'
          IF (IFOUR.EQ.0) THEN
C  Les elements incompressibles ne sont pas utilisables pour le calcul
C  en serie de FOURIER
            NBTEF=22
            IF (NBTEF.GT.N2MAX) GOTO 99
            LESTEF(18)='ICT3'
            LESTEF(19)='ICQ4'
            LESTEF(20)='ICT6'
            LESTEF(21)='ICQ8'
            LESTEF(22)='BBAR'
          ENDIF
        ELSE
          NBTEF=37
          IF (NBTEF.GT.N2MAX) GOTO 99
          LESTEF(1)='TRI3'
          LESTEF(2)='TRI6'
          LESTEF(3)='QUA4'
          LESTEF(4)='QUA8'
          LESTEF(5)='COQ2'
          LESTEF(6)='BARR'
          LESTEF(7)='TRIH'
          LESTEF(8)='JOI2'
          LESTEF(9)='POI1'
          LESTEF(10)='JOI3'
          LESTEF(11)='ICT3'
          LESTEF(12)='ICQ4'
          LESTEF(13)='ICT6'
          LESTEF(14)='ICQ8'
          LESTEF(15)='POLY'
          LESTEF(16)='BAR3'
          LESTEF(17)='QUAH'
          LESTEF(18)='TIMO'
          LESTEF(19)='POUT'
          LESTEF(20)='TRIS'
          LESTEF(21)='QUAS'
          LESTEF(22)='POIS'
          LESTEF(23)='TRH6'
          LESTEF(24)='SEGS'
          LESTEF(25)='POJS'
          LESTEF(26)='Q4RI'
          LESTEF(27)='Q8RI'
          LESTEF(28)='JCT3'
          LESTEF(29)='JCI4'
          LESTEF(30)='JGI2'
          LESTEF(31)='CIFL'
          LESTEF(32)='XQ4R'
          LESTEF(33)='JOI1'
          LESTEF(34)='ZCO2'
          LESTEF(35)='COS2'
          LESTEF(36)='COA2'
          LESTEF(37)='BBAR'
        ENDIF
C  DIMENSION 3 :
C ---------------
      ELSE IF (IDIM.EQ.3) THEN
        NBTEF=56
        IF (NBTEF.GT.N2MAX) GOTO 99
        LESTEF( 1)='CUB8'
        LESTEF( 2)='CU20'
        LESTEF( 3)='PRI6'
        LESTEF( 4)='PR15'
        LESTEF( 5)='TET4'
        LESTEF( 6)='TE10'
        LESTEF( 7)='PYR5'
        LESTEF( 8)='PY13'
        LESTEF( 9)='COQ3'
        LESTEF(10)='COQ4'
        LESTEF(11)='COQ6'
        LESTEF(12)='COQ8'
        LESTEF(13)='DKT '
        LESTEF(14)='BARR'
        LESTEF(15)='POUT'
        LESTEF(16)='TUYA'
        LESTEF(17)='LISP'
        LESTEF(18)='TUFI'
        LESTEF(19)='DST '
        LESTEF(20)='JOI4'
        LESTEF(21)='TUYO'
        LESTEF(22)='TIMO'
        LESTEF(23)='TRIS'
        LESTEF(24)='QUAS'
        LESTEF(25)='POIS'
        LESTEF(26)='LISM'
        LESTEF(27)='JOT3'
        LESTEF(28)='BAR3'
        LESTEF(29)='BAEX'
        LESTEF(30)='LIA2'
        LESTEF(31)='CUBH'
        LESTEF(32)='POJS'
        LESTEF(33)='Q4RI'
        LESTEF(34)='Q8RI'
        LESTEF(35)='JGT3'
        LESTEF(36)='JGI4'
        LESTEF(37)='SEGS'
        LESTEF(38)='SHB8'
        LESTEF(39)='POI1'
        LESTEF(40)='XC8R'
        LESTEF(41)='JOI1'
        LESTEF(42)='ZCO3'
        LESTEF(43)='ZCO4'
        LESTEF(44)='COS2'
        LESTEF(45)='COA2'
        LESTEF(46)='ICC8'
        LESTEF(47)='ICT4'
        LESTEF(48)='ICP6'
        LESTEF(49)='IC20'
        LESTEF(50)='IC10'
        LESTEF(51)='IC15'
        LESTEF(52)='BBAR'
        LESTEF(53)='ICY5'
        LESTEF(54)='IC13'
        LESTEF(55)='C20R'
        LESTEF(56)='P15R'
      ENDIF
      RETURN
 
C  12 - Traitement des ERREURS
C -----------------------------
  99  CALL ERREUR(5)
      RETURN
 
      END