1 : C MATPO2 SOURCE CHAT 05/01/13 01:35:20 5004 2 : SUBROUTINE MATPO2 ( TAU, DT, TEMP, DTEMP, PREDEF, DPRED, WMAT, 3 : & IERUT ) 4 : IMPLICIT INTEGER(I-N) 5 : IMPLICIT REAL*8(A-H,O-Z) 6 : C----------------------------------------------------------------------- 7 : C 8 : C DESCRIPTION FONCTIONNELLE : 9 : C ------------------------- 10 : C 11 : C Dependance UMAPO2 12 : C => Dependance en deuxieme niveau du module Utilisateur UMAT 13 : C Appelee dans le cas suivant : 14 : C ILOI = 22 : Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21 15 : C RESTRICTION A LA FORMULATION MASSIVE 16 : C Instanciation des composantes de materiau au fur et a 17 : C mesure des iterations internes. 18 : C Le modele comprend cette fois des parametres externes. 19 : C Les composantes de materiau sont evaluees par les 20 : C memes fonctions externes que celles utilisees par 21 : C l'operateur VARI. 22 : C MODELE 'MAQUETTE' 23 : C 24 : C Calcul des proprietes de materiau au point TAU 25 : C 26 : C ENTREES : TAU, DT, TEMP, DTEMP, PREDEF, DPRED 27 : C SORTIES : WMAT, IERUT 28 : C 29 : C----------------------------------------------------------------------- 30 : C 31 : C Composantes de materiau 32 : C ----------------------- 33 : C WMAT ( 1) : 'YOUN' WMAT ( 5) : 'SMAX' 34 : C WMAT ( 2) : 'NU ' WMAT ( 6) : 'AF0 ' 35 : C WMAT ( 7) : 'AF1 ' 36 : C WMAT ( 4) : 'ALPH' WMAT ( 8) : 'AF2 ' 37 : C WMAT ( 9) : 'AF3 ' 38 : C WMAT (10) : 'AF4 ' 39 : C WMAT (11) : 'AF5 ' 40 : C WMAT (12) : 'AF6 ' 41 : C 42 : C N.B. WMAT(3) est laissee libre en coherence avec la declaration 43 : C des composantes de materiau du modele, correspond a 'RHO ' 44 : C qui n'est pas utilisee par l'application TOUTATIS 45 : C 46 : C Parametres externes 47 : C ------------------- 48 : C TEMP : 'T ' 49 : C PREDEF( 1) : 'PORO' 50 : C PREDEF( 2) : 'YOGC' 51 : C PREDEF( 3) : 'ALPC' 52 : C PREDEF( 4) : 'TFIS' 53 : C PREDEF( 5) : 'TU02' 54 : C PREDEF( 6) : 'FACF' 55 : C PREDEF( 7) : 'DSIU' 56 : C PREDEF( 8) : 'DGRA' 57 : C 58 : C DTEMP : increment de 'T ' sur le pas de temps 59 : C DPRED ( 1) : increment de 'PORO' sur le pas de temps 60 : C DPRED ( 2) : increment de 'YOGC' sur le pas de temps - INUTILE 61 : C DPRED ( 3) : increment de 'ALPC' sur le pas de temps - INUTILE 62 : C DPRED ( 4) : increment de 'TFIS' sur le pas de temps 63 : C DPRED ( 5) : increment de 'TUO2' sur le pas de temps 64 : C DPRED ( 6) : increment de 'FACF' sur le pas de temps 65 : C DPRED ( 7) : increment de 'DSIU' sur le pas de temps 66 : C DPRED ( 8) : increment de 'DGRA' sur le pas de temps 67 : C 68 : C----------------------------------------------------------------------- 69 : C Arguments de l'interface 70 : C 71 : REAL*8 TAU, DT, TEMP, DTEMP, PREDEF(*), DPRED(*), WMAT(*) 72 : INTEGER IERUT 73 : C 74 : C Variables locales 75 : C 76 : INTEGER NBPAR 77 : REAL*8 VALPAR(5) 78 : REAL*8 DDT, TC, POROC, YOGCC, ALPCC, TFISC, TUO2C, FACFC, 79 : & DSIUC, DGRAC 80 : C 81 : REAL*8 FVYOUN, FVNU, FVALPH, FVSMAX, FVAF0, FVAF1, FVAF2, 82 : & FVAF3, FVAF4, FVAF5, FVAF6 83 : EXTERNAL FVYOUN, FVNU, FVALPH, FVSMAX, FVAF0, FVAF1, FVAF2, 84 : & FVAF3, FVAF4, FVAF5, FVAF6 85 : C 86 : C------------------- Debut du code executable -------------------------- 87 : C 88 : IERUT = 0 89 : C 90 : IF (INT(DPRED(2)).NE.0) THEN 91 : IERUT = -22002 92 : RETURN 93 : ENDIF 94 : IF (INT(DPRED(3)).NE.0) THEN 95 : IERUT = -22003 96 : RETURN 97 : ENDIF 98 : C 99 : C 1 - Evaluation des parametres au point TAU 100 : C 101 : DDT = TAU/DT 102 : C 103 : TC = TEMP + ( DDT * DTEMP ) 104 : POROC = PREDEF(1) + ( DDT * DPRED(1) ) 105 : YOGCC = PREDEF(2) 106 : ALPCC = PREDEF(3) 107 : TFISC = PREDEF(4) + ( DDT * DPRED(4) ) 108 : TUO2C = PREDEF(5) + ( DDT * DPRED(5) ) 109 : FACFC = PREDEF(6) + ( DDT * DPRED(6) ) 110 : DSIUC = PREDEF(7) + ( DDT * DPRED(7) ) 111 : DGRAC = PREDEF(8) + ( DDT * DPRED(8) ) 112 : C 113 : C 2 - Evaluation des proprietes de materiau au point TAU 114 : C 115 : C.....Module d'Young 'YOUN' 116 : C 117 : NBPAR = 3 118 : VALPAR(1) = TC 119 : VALPAR(2) = POROC 120 : VALPAR(3) = YOGCC 121 : WMAT(1) = FVYOUN(VALPAR,NBPAR,IERUT) 122 : IF (IERUT.GT.0) THEN 123 : IERUT = -22000-IERUT 124 : RETURN 125 : ENDIF 126 : C 127 : C.....Coefficient de Poisson 'NU ' 128 : C 129 : WMAT(2) = FVNU(VALPAR,NBPAR,IERUT) 130 : IF (IERUT.GT.0) THEN 131 : IERUT = -22000-IERUT 132 : RETURN 133 : ENDIF 134 : C 135 : C.....Modele de fluage : contrainte de reference 'SMAX' 136 : C 137 : WMAT(5) = FVSMAX(VALPAR,NBPAR,IERUT) 138 : IF (IERUT.GT.0) THEN 139 : IERUT = -22000-IERUT 140 : RETURN 141 : ENDIF 142 : C 143 : C.....Coefficient de dilatation thermique 'ALPH' 144 : C 145 : NBPAR = 2 146 : VALPAR(2) = ALPCC 147 : WMAT(4) = FVALPH(VALPAR,NBPAR,IERUT) 148 : IF (IERUT.GT.0) THEN 149 : IERUT = -22000-IERUT 150 : RETURN 151 : ENDIF 152 : C 153 : C.....Modele de fluage : coefficient du polynome 'AF0 ' 154 : C 155 : NBPAR = 3 156 : VALPAR(1) = TFISC 157 : VALPAR(2) = TUO2C 158 : VALPAR(3) = FACFC 159 : WMAT(6) = FVAF0(VALPAR,NBPAR,IERUT) 160 : IF (IERUT.GT.0) THEN 161 : IERUT = -22000-IERUT 162 : RETURN 163 : ENDIF 164 : C 165 : C.....Modele de fluage : coefficient du polynome 'AF2 ' 166 : C 167 : WMAT(8) = FVAF2(VALPAR,NBPAR,IERUT) 168 : IF (IERUT.GT.0) THEN 169 : IERUT = -22000-IERUT 170 : RETURN 171 : ENDIF 172 : C 173 : C.....Modele de fluage : coefficient du polynome 'AF4 ' 174 : C 175 : WMAT(10) = FVAF4(VALPAR,NBPAR,IERUT) 176 : IF (IERUT.GT.0) THEN 177 : IERUT = -22000-IERUT 178 : RETURN 179 : ENDIF 180 : C 181 : C.....Modele de fluage : coefficient du polynome 'AF5 ' 182 : C 183 : WMAT(11) = FVAF5(VALPAR,NBPAR,IERUT) 184 : IF (IERUT.GT.0) THEN 185 : IERUT = -22000-IERUT 186 : RETURN 187 : ENDIF 188 : C 189 : C.....Modele de fluage : coefficient du polynome 'AF6 ' 190 : C 191 : WMAT(12) = FVAF6(VALPAR,NBPAR,IERUT) 192 : IF (IERUT.GT.0) THEN 193 : IERUT = -22000-IERUT 194 : RETURN 195 : ENDIF 196 : C 197 : C.....Modele de fluage : coefficient du polynome 'AF3 ' 198 : C 199 : NBPAR = 4 200 : VALPAR(4) = DSIUC 201 : WMAT(9) = FVAF3(VALPAR,NBPAR,IERUT) 202 : IF (IERUT.GT.0) THEN 203 : IERUT = -22000-IERUT 204 : RETURN 205 : ENDIF 206 : C 207 : C.....Modele de fluage : coefficient du polynome 'AF1 ' 208 : C 209 : NBPAR = 5 210 : VALPAR(5) = DGRAC 211 : WMAT(7) = FVAF1(VALPAR,NBPAR,IERUT) 212 : IF (IERUT.GT.0) THEN 213 : IERUT = -22000-IERUT 214 : RETURN 215 : ENDIF 216 : C 217 : RETURN 218 : END 219 : 220 :
© Cast3M 2003 - All rights reserved.
Disclaimer