1 : C FLUTRA SOURCE AM 09/12/08 21:15:20 6582 2 : SUBROUTINE FLUTRA(SIG0,NSTRS,DEPST,VAR0,NVARI,XMAT,NMATT, 3 : * IFOUR,TDT,IB,IGAU,TETA1,TETA2,ITHER,NMAT, 4 : * SIGF,VARF,CMAT,DEPSI,SIGI,DEP,EPSILO,EPSEQU, 5 : * EPSFLU,DSIGT,SIGEF,AN,BN,CN,KERRE) 6 : C 7 : C___________________________________________________________________ 8 : C 9 : C MODELE N. REVIRON 10 : C 11 : C DEFINITION DES DIFFERENTES VARIABLES ENTREES, SORTIES, 12 : C PARAMETRES 13 : C 14 : C___________________________________________________________________ 15 : C 16 : C NSTRS = TAILLE DU VECTEUR DES CONTRAINTES 17 : C 18 : C ENTREES 3D 19 : C 20 : C SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES 21 : C 1 : XX 22 : C 2 : YY 23 : C 3 : ZZ 24 : C 4 : XY 25 : C 5 : ZX 26 : C 6 : YZ 27 : C 28 : C SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES 29 : C (SE11,SE22,SE12) 30 : C idem 31 : C 32 : C DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES 33 : C idem 34 : C 35 : C ENTREES 2D CONTRAINTES PLANES 36 : C 37 : C SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES 38 : C 1 : XX 39 : C 2 : YY 40 : C (3 : ZZ) 41 : C 4 : XY 42 : C 43 : C SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES 44 : C (SE11,SE22,SE12) 45 : C idem 46 : C 47 : C DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES 48 : C 1 : XX 49 : C 2 : YY 50 : C (3 : ZZ) 51 : C 4 : XY 52 : C 53 : C DEPSI (UTILSER DANS LES SUBROUTINES) 54 : C 1 : XX 55 : C 2 : YY 56 : C 3 : XY 57 : C (4 : ZZ) 58 : C 59 : C ENTREES 2D DEFORMATIONS PLANES 60 : C 61 : C SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES 62 : C 1 : XX 63 : C 2 : YY 64 : C 3 : ZZ 65 : C 4 : XY 66 : C 67 : C SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES 68 : C (SE11,SE22,SE12) 69 : C idem 70 : C 71 : C DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES 72 : C 1 : XX 73 : C 2 : YY 74 : C (3 : ZZ = 0) 75 : C 4 : XY 76 : C 77 : C ENTREES AXISYMETRIE 78 : C 79 : C SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES 80 : C 1 : RR 81 : C 2 : ZZ 82 : C 3 : TT 83 : C 4 : RZ 84 : C 85 : C SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES 86 : C idem 87 : C 88 : C DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES 89 : C idem 90 : C___________________________________________________________ 91 : C 92 : C VAR0(NVARI) = VECTEUR DES VARIABLES INTERNES INITIALES 93 : C___________________________________________________________ 94 : C 95 : C VARIABLES INTERNES MODELE MECANIQUE ISOTROPE 96 : C 97 : C VAR0(1) = EPSTIM 98 : C VAR0(2) = DINI 99 : C VAR0(3) = EF11 100 : C VAR0(4) = EF22 101 : C VAR0(5) = EF33 102 : C VAR0(6) = EF12 103 : C VAR0(7) = EF31 104 : C VAR0(8) = EF23 105 : C VAR0(9) = SE11 106 : C VAR0(10) = SE22 107 : C VAR0(11) = SE12 108 : C__________________________________________________________ 109 : C 110 : C NVARI = TAILLE DU VECTEUR DE VARIABLES INTERNES 111 : C__________________________________________________________ 112 : C 113 : C XMAT(NVARI) = VECTEUR DES PARAMETRES DU MATERIAU 114 : C 115 : C PARAMETRES STANDARD 116 : C XMAT(1) = YOUN : MODULE D YOUNG 117 : C XMAT(2) = XNU : COEFFICIENT DE POISSON ELASTIQUE 118 : C XMAT(3) = RHO : MASSE VOLUMIQUE 119 : C XMAT(4) = ALPH : COEFFICIENT DE DILATATION THERMIQUE 120 : C 121 : C PARAMETRES MODELE MECANIQUE 122 : C XMAT(5) = KTR0 : KTR0 Mazars 123 : C XMAT(6) = ATRA : ATRA Mazars 124 : C XMAT(7) = BTRA : BTRA Mazars 125 : C XMAT(11) = A : COEF De Schutter Module d YOUNG 126 : C XMAT(12) = B : COEF De Schutter fc : non utilise 127 : C XMAT(13) = C : COEF De Schutter ft 128 : C XMAT(14) = D : COEF De Schutter Gft 129 : C___________________________________________________ 130 : C 131 : C NMATT = TAILLE DU VECTEUR DES PARAMETRES 132 : C MELE = NUMERO DE L ELEMENT FINI UTILISE 133 : C XE = VECTEUR DES COORDONNEES DE L ELEMENT 134 : C NBNN = TAILLE DU VECTEUR XE 135 : C IFOUR = INDICATEUR DU TYPE DE CALCUL 136 : C -3 : DEFORMATION PLANE GENERALISEE 137 : C -2 : CONTRAINTES PLANES 138 : C -1 : DEFORMATION PLANE 139 : C 0 : AXISYMETRIQUE 140 : C 1 : SERIE DE FOURIER 141 : C 2 : TRIDIMENSIONNEL 142 : C MFR = NUMERO DE FORMULATION 143 : C 1 : ELEMENTS MASSIFS 144 : C AUTRE: VOIR nummfr.eso 145 : C ITTHER = INDICATEUR DE CHARGEMENT THERMIQUE 146 : C 0 : PAS DE CHARGEMENT THERMIQUE 147 : C 1 : PRESENCE DE CHARGEMENT THERMIQUE 148 : C TINI = TEMPERATURE EN DEBUT DE PAS 149 : C TFIN = TEMPERATURE EN FIN DE PAS 150 : C TDT = PAS DE TEMPS 151 : C IB = NUMERO DE L ELEMENT CALCULE 152 : C IGAU = NUMERO DU POINT DE GAUSS 153 : C DEP(6,6) = MATRICE ELASTO-PLASTIQUE TANGENTE 154 : C DINI = ENDOMMAGEMENT EN DEBUT DE PAS 155 : C D = ENDOMMAGEMENT EN FIN DE PAS 156 : C SIG0 = VECTEUR DES CONTRAINTES INITIALES APPARENTES 157 : C SIGF = VECTEUR DES CONTRAINTES FINALES APPARENTES 158 : C SIGI = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES 159 : C EPSILO = VECTEUR DES DEFORMATIONS ELASTIQUES FIN PAS 160 : C____________________________________________________ 161 : C 162 : C VARIABLES EN SORTIES 163 : C____________________________________________________ 164 : C 165 : C SIGF(NSTRS) = VECTEUR DES CONTRAINTES FINALES 166 : C REMARQUE : 167 : C SIGF est dans un premier temps les contraintes effectives finales 168 : C puis les contraintes apparentes finales par multiplication par D 169 : C 170 : C DEPSFT = VECTEUR INCREMENT DES DEFORMATIONS DE FLUAGE 171 : C (EF11,EF22,EF33,EF12,EF31,EF23) 172 : C VARF(NVARI) = VECTEUR DES VARIABLES INTERNES FINALES 173 : C KERRE = INDICATEUR D ERREUR 174 : C 0 : PAS D ERREUR 175 : C 176 : C____________________________________________________ 177 : C____________________________________________________ 178 : C____________________________________________________ 179 : C____________________________________________________ 180 : C 181 : C 182 : IMPLICIT INTEGER(I-N) 183 : IMPLICIT REAL*8(A-H,O-Z) 184 : DIMENSION XMAT(NMATT),CMAT(NMAT) 185 : DIMENSION VAR0(NVARI),VARF(NVARI) 186 : DIMENSION DEPST(NSTRS),DEPSI(NSTRS) 187 : DIMENSION SIG0(NSTRS),SIGF(NSTRS),SIGI(NSTRS) 188 : DIMENSION DEP(NSTRS,NSTRS) 189 : DIMENSION EPSILO(NSTRS) 190 : DIMENSION EPSEQU(NSTRS),EPSFLU(NSTRS) 191 : DIMENSION DSIGT(NSTRS),SIGEF(NSTRS) 192 : DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS) 193 : C 194 : C TEST SUR L OPTION DE CALCUL 195 : C LE CALCUL N EST POSSIBLE QUE DANS LE CAS 196 : C CONTRAINTES PLANES : IFOUR = -2 197 : C DEFORMATION PLANE : IFOUR = -1 198 : C AXISYMETRIQUE : IFOUR = 0 199 : C TRIDIMENSIONEL : IFOUR = 2 200 : C 201 : IF ((IFOUR.NE.-2).AND.(IFOUR.NE.-1).AND.(IFOUR.NE.0) 202 : $ .AND.(IFOUR.NE.2)) THEN 203 : WRITE(*,*) 'L OPTION DE CALCUL N EST PAS IMPLEMENTE' 204 : WRITE(*,*) 'IL FAUT CHOISIR PARMIS LES OPTIONS SUIVANTES :' 205 : WRITE(*,*) ' CONTRAINTES PLANES' 206 : WRITE(*,*) ' DEFORMATIONS PLANES' 207 : WRITE(*,*) ' AXISYMETRIQUE' 208 : WRITE(*,*) ' TRIDIMENSIONEL' 209 : KERRE = 99 210 : RETURN 211 : ENDIF 212 : C 213 : C INITIALISATION DES MATRICES 214 : C 215 : XPRECIS = 1.D-6 215 : CALL ZERO(CMAT,NMAT,1) 215 : CALL ZERO(SIGF,NSTRS,1) 216 : C 217 : C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU 218 : C 219 : C 220 : C REMPLISSAGE D UN NOUVEAU VECTEUR DE MATERIAU 221 : C 222 : DO 300 I =1, NMATT 223 : CMAT(I)= XMAT(I) 224 : 300 CONTINUE 225 : C 226 : C AJOUT DU TEMPS A LA FIN 227 : C 228 : CMAT(NMAT) = TDT 229 : C 230 : C GESTION DES TEMPERATURES 231 : C 232 : CMAT(29) = TETA1 233 : CMAT(30) = TETA2 234 : C 235 : C RECUPERATION DES VARIABLES INTERNES INITIALES 236 : C 237 : C REMPLISSAGE DES VARIABLES INTERNES FINALES 238 : C 239 : DO 200 I =1, NVARI 240 : VARF(I)= VAR0(I) 241 : 200 CONTINUE 242 : C 243 : C_______________________________________________________________________ 244 : C 245 : C REDIMENSIONNEMENT DES VECTEURS DANS LE CAS DE CONTRAINTES PLANES 246 : C_______________________________________________________________________ 247 : C 248 : C 249 : C 250 : C CHANGEMENT DE NOM POUR TRAITEMENT 251 : C 252 : DO 100 I =1, NSTRS 253 : DEPSI(I)= DEPST(I) 254 : 100 CONTINUE 255 : C 256 : IPLAN = 0 257 : IF ((NSTRS.EQ.4).AND.(IFOUR.EQ.-2)) THEN 258 : IPLAN = 1 259 : DEPSI(3) = DEPST(4) 260 : DEPSI(4) = DEPST(3) 261 : IF (DABS(SIG0(3)).GE.1.D-3) THEN 262 : WRITE(*,*)'PROBLEME AVEC LES CONTRAINTES PLANES' 263 : KERRE=99 264 : RETURN 265 : ENDIF 266 : ENDIF 267 : C_______________________________________________________________________ 268 : C 269 : C MODELE MECANIQUE 270 : C_______________________________________________________________________ 271 : C 271 : CALL CAMECA(NSTRS,DEPSI,VAR0,NVARI,CMAT,NMAT, 272 : * IFOUR,IB,IGAU,SIGF,VARF,SIGI,DEP, 273 : * EPSILO,EPSEQU,EPSFLU,DSIGT,SIGEF, 274 : * AN,BN,CN,KERRE) 275 : C 276 : C 277 : C_______________________________________________________________________ 278 : C 279 : C REDIMENSIONNEMENT DES VECTEURS DANS LE CAS DE CONTRAINTES PLANES 280 : C_______________________________________________________________________ 281 : C 282 : IF (IPLAN.EQ.1) THEN 283 : temp = SIGF(4) 284 : SIGF(4)=SIGF(3) 285 : SIGF(3)=temp 286 : ENDIF 287 : C 288 : END 289 : 290 :
© Cast3M 2003 - All rights reserved.
Disclaimer