Numérotation des lignes :

C MODEL2    SOURCE    PASCAL    20/10/07    21:15:11     10746           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 LINEAIRESC --------------------------------------      CALL MODLIN(MOMODL,NMOD)      DO i=1,NMOD        LESPRO(i)=MOMODL(i)      ENDDO      NBPRO=NMOD C  2 - Noms des MATERIAUx pour le comportement LINEAIREC ------------------------------------------------------      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 LINEAIRESC ------------------------------------------      CALL MODNLI(MOMODL,NMOD)      DO i=NBPRO+1,NBPRO+NMOD        LESPRO(i)=MOMODL(i-NBPRO)      ENDDO      NBPRO=NBPRO+NMOD C  4 - Noms des modeles PLASTIQUESC ---------------------------------      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 FLUAGEC --------------------------------      CALL MODFLU(MOMODL,NMOD)      DO i=NBPRO+1,NBPRO+NMOD        LESPRO(i)=MOMODL(i-NBPRO)      ENDDO      NBPRO=NBPRO+NMOD C  6 - Noms des modeles VISCOPLASTIQUESC --------------------------------------      CALL MODVIS(MOMODL,NMOD)      DO i=NBPRO+1,NBPRO+NMOD        LESPRO(i)=MOMODL(i-NBPRO)      ENDDO      NBPRO=NBPRO+NMOD C  7 - Noms des modeles ENDOMMAGEABLESC -------------------------------------      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 ENDOMMAGEABLESC ------------------------------------------------      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 LINEAIRESC -----------------------------------------------      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 EXTERNESC ------------------------------------------------      CALL MODVIX(MOMODL,NMOD)      DO i=NBPRO+1,NBPRO+NMOD        LESPRO(i)=MOMODL(i-NBPRO)      ENDDO      NBPRO=NBPRO+NMOD C  12 - Noms des modeles LIAISONSC ------------------------------------------------      CALL MODLIA(MOMODL,NMOD)      DO i=NBPRO+1,NBPRO+NMOD        LESPRO(i)=MOMODL(i-NBPRO)      ENDDO      NBPRO=NBPRO+NMOD C  13 - Noms des modeles IMPEDANCEC ------------------------------------------------      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 FINISC -------------------------------------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) THENC  Les elements incompressibles ne sont pas utilisables pour le calculC  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'        ENDIFC  DIMENSION 3 :C ---------------      ELSE IF (IDIM.EQ.3) THEN        NBTEF=54        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'      ENDIF      RETURN C  12 - Traitement des ERREURSC -----------------------------  99  CALL ERREUR(5)      RETURN       END

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