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 : -INC CCOPTIO 10 : * 11 : CHARACTER*16 MOMODL(100) 12 : * 13 : CHARACTER*16 LESPRO(*) 14 : CHARACTER*4 LESTEF(*) 15 : * 16 : * Noms des comportements lineaires 17 : * 17 : CALL MODLIN (MOMODL,NMOD) 18 : * 19 : DO 1 I=1,NMOD 20 : LESPRO(I)=MOMODL(I) 21 : 1 CONTINUE 22 : NBPRO=NMOD 23 : * 24 : * Noms des materiaux pour le comportement lineaire 25 : * 25 : CALL MODELA (MOMODL,NMOD) 26 : * 27 : DO 2 I=NBPRO+1,NBPRO+NMOD 28 : LESPRO(I)=MOMODL(I-NBPRO) 29 : 2 CONTINUE 30 : NBPRO=NBPRO+NMOD 31 : * 32 : * Noms des comportements non lineaires 33 : * 34 : * Remarque 35 : * La routine MODNLI va retourner les noms des sept familles de 36 : * comportements non lineaires, dans l'ordre suivant : 37 : * 'PLASTIQUE', 'FLUAGE', 'VISCOPLASTIQUE', 'ENDOMMAGEMENT', 38 : * 'PLASTIQUE_ENDOM', 'NON_LINEAIRE' et 'VISCO_EXTERNE'. 39 : * Les noms de toutes les familles de comportements non lineaires 40 : * feront donc partie du vocabulaire admissible dans la definition 41 : * du modele, sans restriction. 42 : * Toutefois, DANS LE CAS DE LA FORMULATION 'POREUX', seuls les 43 : * materiaux des 6 PREMIERES FAMILLES seront acceptables : cf. appels 44 : * ci-dessous aux routines MODPLA, MODFLU, MODVIS et MODEND. 45 : * Il y a donc une restriction en formulation 'POREUX' par rapport 46 : * a la formulation 'MECANIQUE', mais on ne la fait pas apparaitre 47 : * au stade de la definition des familles admissibles. 48 : 48 : CALL MODNLI (MOMODL,NMOD) 49 : * 50 : DO 3 I=NBPRO+1,NBPRO+NMOD 51 : LESPRO(I)=MOMODL(I-NBPRO) 52 : 3 CONTINUE 53 : NBPRO=NBPRO+NMOD 54 : 55 : * 56 : * Noms des materiaux plastiques 57 : * 57 : CALL MODPLA (MOMODL,NMOD) 58 : * 59 : DO 4 I=NBPRO+1,NBPRO+NMOD 60 : LESPRO(I)=MOMODL(I-NBPRO) 61 : 4 CONTINUE 62 : NBPRO=NBPRO+NMOD 63 : * 64 : * Noms des materiaux fluages 65 : * 65 : CALL MODFLU (MOMODL,NMOD) 66 : * 67 : DO 5 I=NBPRO+1,NBPRO+NMOD 68 : LESPRO(I)=MOMODL(I-NBPRO) 69 : 5 CONTINUE 70 : NBPRO=NBPRO+NMOD 71 : * 72 : * Noms des materiaux viscoplastiques 73 : * 73 : CALL MODVIS (MOMODL,NMOD) 74 : * 75 : DO 6 I=NBPRO+1,NBPRO+NMOD 76 : LESPRO(I)=MOMODL(I-NBPRO) 77 : 6 CONTINUE 78 : NBPRO=NBPRO+NMOD 79 : * 80 : * Noms des materiaux endommageables 81 : * 81 : CALL MODEND (MOMODL,NMOD) 82 : * 83 : DO 7 I=NBPRO+1,NBPRO+NMOD 84 : LESPRO(I)=MOMODL(I-NBPRO) 85 : 7 CONTINUE 86 : NBPRO=NBPRO+NMOD 87 : * 88 : * Noms des materiaux plastiques endommageables 89 : * 89 : CALL MODPLE (MOMODL,NMOD) 90 : * 91 : DO 8 I=NBPRO+1,NBPRO+NMOD 92 : LESPRO(I)=MOMODL(I-NBPRO) 93 : 8 CONTINUE 94 : NBPRO=NBPRO+NMOD 95 : * 96 : * Noms des materiaux elastiques non linéaires 97 : * 97 : CALL MODENL (MOMODL,NMOD) 98 : * 99 : DO 9 I=NBPRO+1,NBPRO+NMOD 100 : LESPRO(I)=MOMODL(I-NBPRO) 101 : 9 CONTINUE 102 : NBPRO=NBPRO+NMOD 103 : * 104 : IF(NBPRO.GT.N1MAX) THEN 104 : CALL ERREUR(550) 105 : GO TO 99 106 : END IF 107 : * 108 : * NOMS DE TYPES D'ELEMENTS FINIS 109 : * 110 : IF(IDIM.EQ.2) THEN 111 : NBTEF=9 112 : IF(NBTEF.GT.N2MAX) THEN 112 : CALL ERREUR(551) 113 : GO TO 99 114 : END IF 115 : LESTEF(1)='TRIP' 116 : LESTEF(2)='QUAP' 117 : LESTEF(3)='JOP3' 118 : LESTEF(4)='TRIQ' 119 : LESTEF(5)='QUAQ' 120 : LESTEF(6)='JOQ3' 121 : LESTEF(7)='TRIR' 122 : LESTEF(8)='QUAR' 123 : LESTEF(9)='JOR3' 124 : ELSE 125 : NBTEF=15 126 : IF(NBTEF.GT.N2MAX) THEN 126 : CALL ERREUR(551) 127 : GO TO 99 128 : END IF 129 : LESTEF(1)='CUBP' 130 : LESTEF(2)='TETP' 131 : LESTEF(3)='PRIP' 132 : LESTEF(4)='JOP6' 133 : LESTEF(5)='JOP8' 134 : LESTEF(6)='CUBQ' 135 : LESTEF(7)='TETQ' 136 : LESTEF(8)='PRIQ' 137 : LESTEF(9)='JOQ6' 138 : LESTEF(10)='JOQ8' 139 : LESTEF(11)='CUBR' 140 : LESTEF(12)='TETR' 141 : LESTEF(13)='PRIR' 142 : LESTEF(14)='JOR6' 143 : LESTEF(15)='JOR8' 144 : ENDIF 145 : RETURN 145 : 99 CALL ERREUR(5) 146 : RETURN 147 : END 148 : 149 : 150 : 151 : 152 :
© Cast3M 2003 - All rights reserved.
Disclaimer