Sources Esope | Cast3M

model2
C MODEL2    SOURCE    MB234859  25/08/26    21:15:11     12343          
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      IF (NBPRO.GT.N1MAX) GOTO 99
 
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
      NBPRO=NBPRO + 2
      IF (NBPRO.GT.N1MAX) GOTO 99
C     Restauration du comportement si FUSION
      LESPRO(NBPRO-1) = 'FUSION'
C     CRItere Pression
      LESPRO(NBPRO)   = 'CRIP'
 
C  14 - Noms de types d'ELEMENTS FINIS
C -------------------------------------
C     -----------
C     Elements 1D
C     -----------
CCC   IF (IFOMOD.GE.3.AND.IFOMOD.LE.5) THEN
      IF (IDIM.EQ.1) THEN
        NBTEF=2
        IF (NBTEF.GT.N2MAX) GOTO 99
        LESTEF(1)='M1D2'
        LESTEF(2)='M1D3'
      ELSE
C       -----------------
C       Elements 2D et 3D
C       -----------------
        NBTEF=18
        IF (NBTEF.GT.N2MAX) GOTO 99
        LESTEF( 1)='SEG2'
        LESTEF( 2)='SEG3'
        LESTEF( 3)='TRI3'
        LESTEF( 4)='TRI6'
        LESTEF( 5)='QUA4'
        LESTEF( 6)='QUA8'
        LESTEF( 7)='BARR'
        LESTEF( 8)='TRIH'
        LESTEF( 9)='JOI2'
        LESTEF(10)='JOI3'
        LESTEF(11)='BAR3'
        LESTEF(12)='QUAH'
        LESTEF(13)='Q4RI'
        LESTEF(14)='Q8RI'
        LESTEF(15)='XQ4R'
        LESTEF(16)='ZCO2'
        LESTEF(17)='COS2'
        LESTEF(18)='COA2'
        IF (IFOMOD.LT.2) THEN
          NBTEF=NBTEF+2
          IF (NBTEF.GT.N2MAX) GOTO 99
          LESTEF(NBTEF-1)='COQ2'
          LESTEF(NBTEF  )='POLY'
          IF (IFOMOD.EQ.0.OR.IFOMOD.EQ.1) THEN
            NBTEF=NBTEF+1
            IF (NBTEF.GT.N2MAX) GOTO 99
            LESTEF(NBTEF)='CERC'
          ELSE IF (IFOMOD.EQ.-1) THEN
            NBTEF=NBTEF+14
            IF (NBTEF.GT.N2MAX) GOTO 99
            LESTEF(NBTEF-13)='POI1'
            LESTEF(NBTEF-12)='TIMO'
            LESTEF(NBTEF-11)='POUT'
            LESTEF(NBTEF-10)='JOI1'
            LESTEF(NBTEF- 9)='TRIS'
            LESTEF(NBTEF- 8)='QUAS'
            LESTEF(NBTEF- 7)='POIS'
            LESTEF(NBTEF- 6)='SEGS'
            LESTEF(NBTEF- 5)='POJS'
            LESTEF(NBTEF- 4)='TRH6'
            LESTEF(NBTEF- 3)='JCT3'
            LESTEF(NBTEF- 2)='JCI4'
            LESTEF(NBTEF- 1)='JGI2'
            LESTEF(NBTEF   )='CIFL'
          ENDIF
CCC     ELSE  IF (IFOMOD.EQ.2) THEN
        ELSE
          NBTEF=NBTEF+43
          IF (NBTEF.GT.N2MAX) GOTO 99
          LESTEF(NBTEF-42)='CUB8'
          LESTEF(NBTEF-41)='CU20'
          LESTEF(NBTEF-40)='PRI6'
          LESTEF(NBTEF-39)='PR15'
          LESTEF(NBTEF-38)='TET4'
          LESTEF(NBTEF-37)='TE10'
          LESTEF(NBTEF-36)='PYR5'
          LESTEF(NBTEF-35)='PY13'
          LESTEF(NBTEF-34)='COQ3'
          LESTEF(NBTEF-33)='COQ4'
          LESTEF(NBTEF-32)='COQ6'
          LESTEF(NBTEF-31)='COQ8'
          LESTEF(NBTEF-30)='DKT '
          LESTEF(NBTEF-29)='POUT'
          LESTEF(NBTEF-28)='TUYA'
          LESTEF(NBTEF-27)='LISP'
          LESTEF(NBTEF-26)='TUFI'
          LESTEF(NBTEF-25)='DST '
          LESTEF(NBTEF-24)='JOI4'
          LESTEF(NBTEF-23)='TUYO'
          LESTEF(NBTEF-22)='TIMO'
          LESTEF(NBTEF-21)='TRIS'
          LESTEF(NBTEF-20)='QUAS'
          LESTEF(NBTEF-19)='POIS'
          LESTEF(NBTEF-18)='LISM'
          LESTEF(NBTEF-17)='JOT3'
          LESTEF(NBTEF-16)='BAEX'
          LESTEF(NBTEF-15)='LIA2'
          LESTEF(NBTEF-14)='CUBH'
          LESTEF(NBTEF-13)='POJS'
          LESTEF(NBTEF-12)='Q4RI'
          LESTEF(NBTEF-11)='Q8RI'
          LESTEF(NBTEF-10)='JGT3'
          LESTEF(NBTEF- 9)='JGI4'
          LESTEF(NBTEF- 8)='SEGS'
          LESTEF(NBTEF- 7)='SHB8'
          LESTEF(NBTEF- 6)='POI1'
          LESTEF(NBTEF- 5)='XC8R'
          LESTEF(NBTEF- 4)='JOI1'
          LESTEF(NBTEF- 3)='ZCO3'
          LESTEF(NBTEF- 2)='ZCO4'
          LESTEF(NBTEF- 1)='C20R'
          LESTEF(NBTEF   )='P15R'
        ENDIF
C       -----------
C       Elements incompressibles : inutilisables en serie de FOURIER
C       -----------
        IF (IFOMOD.NE.1) THEN
          NBTEF=NBTEF+13
          IF (NBTEF.GT.N2MAX) GOTO 99
          LESTEF(NBTEF-12)='ICT3'
          LESTEF(NBTEF-11)='ICQ4'
          LESTEF(NBTEF-10)='ICT6'
          LESTEF(NBTEF- 9)='ICQ8'
          LESTEF(NBTEF- 8)='ICC8'
          LESTEF(NBTEF- 7)='ICT4'
          LESTEF(NBTEF- 6)='ICP6'
          LESTEF(NBTEF- 5)='IC20'
          LESTEF(NBTEF- 4)='IC10'
          LESTEF(NBTEF- 3)='IC15'
          LESTEF(NBTEF- 2)='ICY5'
          LESTEF(NBTEF- 1)='IC13'
          LESTEF(NBTEF   )='BBAR'
        ENDIF
      ENDIF
      RETURN
 
C  15 - Traitement des ERREURS
C -----------------------------
  99  CALL ERREUR(5)
 
      END