1 : C CMAXOT SOURCE KICH 09/09/11 21:16:16 6505 2 : SUBROUTINE CMAXOT(wrk52,wrk53,wrk54,WR12,IB,IGAU,NCHAIN,DT,NWA, 3 : & EPSFLU) 4 : * 5 : * 6 : *============================================================== 7 : * modele maxott 8 : *============================================================== 9 : * 10 : * entrees 11 : * 12 : * WRK12 13 : * EMi(2,NWA(i+1)) i=0 a 8 14 : * = EVOLUTION DU MODULE D'YOUNG DE LA BRANCHE i 15 : * SMi(NSTRS) i=1 a 8 16 : * = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION 17 : * 18 : * IVALMA(NCXMAT) = INDEX CARACTERISTIQUES MECANIQUES DU MATERIAU 19 : * 20 : * 21 : * IB = NUMERO DE L'ELEMENT 22 : * IGAU = NUMERO DU POINT DE GAUSS 23 : * NCHAIN = NOMBRE DE BRANCHE DU MODELE DE MAXWELL 24 : * DT = INCREMENT DE TEMPS 25 : * NWA = DIMENSION DES EVOLUTIONS EMi 26 : * EPSFLU = DEFORMATION DE FLUAGE 27 : * 28 : * sorties 29 : * 30 : * SIGF(NSTRS)= CONTR. A LA FIN DU PAS D'INTEGRATION 31 : * DEFP(NSTRS)= INCREMENT DES DEFORM. PLASTIQUES A LA FIN DU PAS 32 : * D'INTEGRATION 33 : * EPINF(NSTRS) = DEFORMATION INELASTIQUE A LA FIN DU PAS D'INTEGRATION 34 : * 35 : *============================================================== 36 : * 37 : IMPLICIT INTEGER(I-N) 38 : IMPLICIT REAL*8(A-H,O-Z) 39 : -INC CCOPTIO 40 : -INC DECHE 41 : * 42 : * 43 : SEGMENT WR12 44 : REAL*8 EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3)) 45 : REAL*8 EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6)) 46 : REAL*8 EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9)) 47 : REAL*8 SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS) 48 : REAL*8 SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS) 49 : REAL*8 SM8(NSTRS) 50 : ENDSEGMENT 51 : * 52 : * 53 : DIMENSION NWA(9) 54 : DIMENSION EPSFLU(8) 55 : DIMENSION EPSOTT(8) 56 : DIMENSION DEPFIS(8) 57 : DIMENSION DEPSF(8) 58 : DIMENSION SFLU(8) 59 : DIMENSION SIGM(8) 60 : REAL*8 ETILDE,COE 61 : INTEGER IND 62 : 63 : NCXMAT = XMAT(/1) 64 : NSTRS = SIG0(/1) 65 : NVARI = VAR0(/1) 66 : ICARA=XCARB(/1) 67 : ANU=XMAT(2) 68 : UNANU=1.D0-ANU 69 : TPS1=TEMP0 70 : TPS2=TEMP0+DT 71 : SE1=0.D0 72 : SE2=0.D0 73 : ETILDE = 0.D0 74 : 75 : IF (IFOUR.EQ.2.AND.MFR.EQ.1) THEN 76 : IND = 21 77 : ELSE IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-3 78 : & .OR.IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN 79 : IND = 16 80 : ELSE 81 : IND = 13 82 : ENDIF 83 : 84 : *====================================================================== 85 : * Formulation TRIDIMENSIONNELLE MASSIVE 86 : * ===================================== 87 : * Formulation DEFORMATIONS PLANES /AXISYMETRIQUE 88 : * ============================================== 89 : *---------------------------------------------------------------------- 90 : * 91 : IF ((IFOUR.EQ.2.AND.MFR.EQ.1).OR. 92 : & (IFOMOD.EQ.-1.AND.IFOUR.NE.-2).OR. 93 : & (IFOMOD.EQ.0.OR.IFOMOD.EQ.1)) THEN 94 : CONS=1.D0/((1.D0+ANU)*(1.D0-2.D0*ANU)) 95 : DO M=0,NCHAIN-1 96 : * 97 : *---------------------------------------------------------------------- 98 : * module d'YOUNG apparent 99 : *---------------------------------------------------------------------- 100 : 100 : CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2) 101 : 102 : ETILDE = ETILDE + EFIC 103 : 104 : SE1=SE1+E1 105 : SE2=SE2+E2 106 : 107 : IF (E1.EQ.0.) THEN 108 : COE = 1. 109 : ELSE 110 : COE = E2 / E1 111 : ENDIF 112 : * 113 : *---------------------------------------------------------------------- 114 : * deformations de fluage 115 : *---------------------------------------------------------------------- 116 : * 117 : FLUA=EXP(-TR*DT) 118 : DO II=1,NSTRS 119 : IF (M.EQ.0) THEN 120 : SFLU(II)=SM0(II)*(COE*FLUA-1.0D0) 121 : ELSE IF (M.EQ.1) THEN 122 : SFLU(II)=SFLU(II)+SM1(II)*(COE*FLUA-1.0D0) 123 : ELSE IF (M.EQ.2) THEN 124 : SFLU(II)=SFLU(II)+SM2(II)*(COE*FLUA-1.0D0) 125 : ELSE IF (M.EQ.3) THEN 126 : SFLU(II)=SFLU(II)+SM3(II)*(COE*FLUA-1.0D0) 127 : ELSE IF (M.EQ.4) THEN 128 : SFLU(II)=SFLU(II)+SM4(II)*(COE*FLUA-1.0D0) 129 : ELSE IF (M.EQ.5) THEN 130 : SFLU(II)=SFLU(II)+SM5(II)*(COE*FLUA-1.0D0) 131 : ELSE IF (M.EQ.6) THEN 132 : SFLU(II)=SFLU(II)+SM6(II)*(COE*FLUA-1.0D0) 133 : ELSE IF (M.EQ.7) THEN 134 : SFLU(II)=SFLU(II)+SM7(II)*(COE*FLUA-1.0D0) 135 : ELSE IF (M.EQ.8) THEN 136 : SFLU(II)=SFLU(II)+SM8(II)*(COE*FLUA-1.0D0) 137 : ENDIF 138 : END DO 139 : END DO 140 : * 141 : EPSFLU(1)=(SFLU(1) - ANU*SFLU(2) - ANU*SFLU(3)) / ETILDE 142 : EPSFLU(2)=(SFLU(2) - ANU*SFLU(3) - ANU*SFLU(1)) / ETILDE 143 : EPSFLU(3)=(SFLU(3) - ANU*SFLU(1) - ANU*SFLU(2)) / ETILDE 144 : DO II=4,NSTRS 145 : EPSFLU(II)=(SFLU(II)*(1.0D0+ANU))/ETILDE 146 : END DO 147 : * 148 : *---------------------------------------------------------------------- 149 : * stockage du module d'elasticite 150 : *---------------------------------------------------------------------- 151 : * 152 : VARF(IND) = ETILDE 153 : XMAT(1) = ETILDE 154 : * 155 : *---------------------------------------------------------------------- 156 : * deformation corrigee pour Ottosen 157 : *---------------------------------------------------------------------- 158 : * 159 : DO I=1,NSTRS 160 : EPSOTT(I) = DEPST(I) + EPSFLU(I) 161 : END DO 162 : * 163 : *---------------------------------------------------------------------- 164 : * appel a Ottosen 165 : *---------------------------------------------------------------------- 166 : 167 : * 167 : CALL OTTOSE(INPLAS,SIG0,NSTRS,EPSOTT,VAR0,XMAT,IVALMA,NCXMAT, 168 : & XCARB,ICARA,NVARI,SIGF,VARF,DEPFIS,MFR,KERRE,IB,IGAU) 169 : * 170 : *---------------------------------------------------------------------- 171 : * deformation corrigee pour Maxwell 172 : *---------------------------------------------------------------------- 173 : * 174 : DO II=1,NSTRS 175 : DEPSF(II) = DEPST(II) - DEPFIS(II) 176 : END DO 177 : * 178 : *---------------------------------------------------------------------- 179 : * modele de Maxwell 180 : *---------------------------------------------------------------------- 181 : * 182 : DO M=0,NCHAIN-1 182 : CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2) 183 : IF (E1.EQ.0.) THEN 184 : COE = 1. 185 : ELSE 186 : COE = E2 / E1 187 : ENDIF 188 : FLUA=EXP(-TR*DT) 189 : IF (M.EQ.0) THEN 190 : SM0(1)=SM0(1)*COE*FLUA 191 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 192 : SM0(2)=SM0(2)*COE*FLUA 193 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 194 : SM0(3)=SM0(3)*COE*FLUA 195 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 196 : DO II = 4,NSTRS 197 : SM0(II)=SM0(II)*COE*FLUA 198 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 199 : END DO 200 : DO II=1,NSTRS 201 : SIGM(II)=SM0(II) 202 : END DO 203 : ELSE IF (M.EQ.1) THEN 204 : SM1(1)=SM1(1)*COE*FLUA 205 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 206 : SM1(2)=SM1(2)*COE*FLUA 207 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 208 : SM1(3)=SM1(3)*COE*FLUA 209 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 210 : DO II = 4,NSTRS 211 : SM1(II)=SM1(II)*COE*FLUA 212 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 213 : END DO 214 : DO II=1,NSTRS 215 : SIGM(II)=SIGM(II)+SM1(II) 216 : END DO 217 : ELSE IF (M.EQ.2) THEN 218 : SM2(1)=SM2(1)*COE*FLUA 219 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 220 : SM2(2)=SM2(2)*COE*FLUA 221 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 222 : SM2(3)=SM2(3)*COE*FLUA 223 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 224 : DO II = 4,NSTRS 225 : SM2(II)=SM2(II)*COE*FLUA 226 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 227 : END DO 228 : DO II=1,NSTRS 229 : SIGM(II)=SIGM(II)+SM2(II) 230 : END DO 231 : ELSE IF (M.EQ.3) THEN 232 : SM3(1)=SM3(1)*COE*FLUA 233 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 234 : SM3(2)=SM3(2)*COE*FLUA 235 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 236 : SM3(3)=SM3(3)*COE*FLUA 237 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 238 : DO II = 4,NSTRS 239 : SM3(II)=SM3(II)*COE*FLUA 240 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 241 : END DO 242 : DO II=1,NSTRS 243 : SIGM(II)=SIGM(II)+SM3(II) 244 : END DO 245 : ELSE IF (M.EQ.4) THEN 246 : SM4(1)=SM4(1)*COE*FLUA 247 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 248 : SM4(2)=SM4(2)*COE*FLUA 249 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 250 : SM4(3)=SM4(3)*COE*FLUA 251 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 252 : DO II = 4,NSTRS 253 : SM4(II)=SM4(II)*COE*FLUA 254 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 255 : END DO 256 : DO II=1,NSTRS 257 : SIGM(II)=SIGM(II)+SM4(II) 258 : END DO 259 : ELSE IF (M.EQ.5) THEN 260 : SM5(1)=SM5(1)*COE*FLUA 261 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 262 : SM5(2)=SM5(2)*COE*FLUA 263 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 264 : SM5(3)=SM5(3)*COE*FLUA 265 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 266 : DO II = 4,NSTRS 267 : SM5(II)=SM5(II)*COE*FLUA 268 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 269 : END DO 270 : DO II=1,NSTRS 271 : SIGM(II)=SIGM(II)+SM5(II) 272 : END DO 273 : ELSE IF (M.EQ.6) THEN 274 : SM6(1)=SM6(1)*COE*FLUA 275 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 276 : SM6(2)=SM6(2)*COE*FLUA 277 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 278 : SM6(3)=SM6(3)*COE*FLUA 279 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 280 : DO II = 4,NSTRS 281 : SM6(II)=SM6(II)*COE*FLUA 282 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 283 : END DO 284 : DO II=1,NSTRS 285 : SIGM(II)=SIGM(II)+SM6(II) 286 : END DO 287 : ELSE IF (M.EQ.7) THEN 288 : SM7(1)=SM7(1)*COE*FLUA 289 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 290 : SM7(2)=SM7(2)*COE*FLUA 291 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 292 : SM7(3)=SM7(3)*COE*FLUA 293 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 294 : DO II = 4,NSTRS 295 : SM7(II)=SM7(II)*COE*FLUA 296 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 297 : END DO 298 : DO II=1,NSTRS 299 : SIGM(II)=SIGM(II)+SM7(II) 300 : END DO 301 : ELSE IF (M.EQ.8) THEN 302 : SM8(1)=SM8(1)*COE*FLUA 303 : 1 +CONS*EFIC*((UNANU)*DEPSF(1)+ANU*(DEPSF(2)+DEPSF(3))) 304 : SM8(2)=SM8(2)*COE*FLUA 305 : 1 +CONS*EFIC*((UNANU)*DEPSF(2)+ANU*(DEPSF(1)+DEPSF(3))) 306 : SM8(3)=SM8(3)*COE*FLUA 307 : 1 +CONS*EFIC*((UNANU)*DEPSF(3)+ANU*(DEPSF(1)+DEPSF(2))) 308 : DO II = 4,NSTRS 309 : SM8(II)=SM8(II)*COE*FLUA 310 : 1 +CONS*EFIC*(1.D0-2.D0*ANU)*DEPSF(II) 311 : END DO 312 : DO II=1,NSTRS 313 : SIGM(II)=SIGM(II)+SM8(II) 314 : END DO 315 : ENDIF 316 : END DO 317 : DEFP(1)=DEPSF(1)-((SIGM(1)-ANU*(SIGM(2)+SIGM(3)))/SE2 318 : 1 -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1) 319 : DEFP(2)=DEPSF(2)-((SIGM(2)-ANU*(SIGM(1)+SIGM(3)))/SE2 320 : 1 -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1) 321 : DEFP(3)=DEPSF(3)-((SIGM(3)-ANU*(SIGM(2)+SIGM(1)))/SE2 322 : 1 -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1) 323 : DO II = 4,NSTRS 324 : DEFP(II)=DEPSF(II)-(1+ANU)*(SIGM(II)/SE2-SIG0(II)/SE1) 325 : END DO 326 : DO I=1,NSTRS 327 : EPINF(I)=EPIN0(I)+DEFP(I)+DEPFIS(I) 328 : END DO 329 : * 330 : *====================================================================== 331 : * Formulation CONTRAINTES PLANES 332 : * ============================== 333 : *---------------------------------------------------------------------- 334 : * 335 : ELSE IF (IFOMOD.EQ.-1.AND.IFOUR.EQ.-2) THEN 336 : CONS=1.D0/(1.D0-(ANU*ANU)) 337 : DO M=0,NCHAIN-1 338 : * 339 : *---------------------------------------------------------------------- 340 : * module d'YOUNG apparent 341 : *---------------------------------------------------------------------- 342 : 342 : CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2) 343 : ETILDE = ETILDE + EFIC 344 : SE1=SE1+E1 345 : SE2=SE2+E2 346 : IF (E1.EQ.0.) THEN 347 : COE = 1. 348 : ELSE 349 : COE = E2 / E1 350 : ENDIF 351 : * 352 : *---------------------------------------------------------------------- 353 : * deformations de fluage 354 : *---------------------------------------------------------------------- 355 : * 356 : FLUA=EXP(-TR*DT) 357 : DO II=1,NSTRS 358 : IF (M.EQ.0) THEN 359 : SFLU(II)=SM0(II)*(COE*FLUA-1.0D0) 360 : ELSE IF (M.EQ.1) THEN 361 : SFLU(II)=SFLU(II)+SM1(II)*(COE*FLUA-1.0D0) 362 : ELSE IF (M.EQ.2) THEN 363 : SFLU(II)=SFLU(II)+SM2(II)*(COE*FLUA-1.0D0) 364 : ELSE IF (M.EQ.3) THEN 365 : SFLU(II)=SFLU(II)+SM3(II)*(COE*FLUA-1.0D0) 366 : ELSE IF (M.EQ.4) THEN 367 : SFLU(II)=SFLU(II)+SM4(II)*(COE*FLUA-1.0D0) 368 : ELSE IF (M.EQ.5) THEN 369 : SFLU(II)=SFLU(II)+SM5(II)*(COE*FLUA-1.0D0) 370 : ELSE IF (M.EQ.6) THEN 371 : SFLU(II)=SFLU(II)+SM6(II)*(COE*FLUA-1.0D0) 372 : ELSE IF (M.EQ.7) THEN 373 : SFLU(II)=SFLU(II)+SM7(II)*(COE*FLUA-1.0D0) 374 : ELSE IF (M.EQ.8) THEN 375 : SFLU(II)=SFLU(II)+SM8(II)*(COE*FLUA-1.0D0) 376 : ENDIF 377 : END DO 378 : END DO 379 : * 380 : EPSFLU(1)=(SFLU(1) - ANU*SFLU(2) - ANU*SFLU(3)) / ETILDE 381 : EPSFLU(2)=(SFLU(2) - ANU*SFLU(3) - ANU*SFLU(1)) / ETILDE 382 : EPSFLU(3)=(SFLU(3) - ANU*SFLU(1) - ANU*SFLU(2)) / ETILDE 383 : EPSFLU(4)=(1.0D0+ANU) * SFLU(4) / ETILDE 384 : * 385 : *---------------------------------------------------------------------- 386 : * stockage du module d'elasticite 387 : *---------------------------------------------------------------------- 388 : * 389 : VARF(IND) = ETILDE 390 : XMAT(1) = ETILDE 391 : * 392 : *---------------------------------------------------------------------- 393 : * deformation corrigee pour Ottosen 394 : *---------------------------------------------------------------------- 395 : * 396 : DO I=1,NSTRS 397 : EPSOTT(I) = DEPST(I) + EPSFLU(I) 398 : END DO 399 : * 400 : *---------------------------------------------------------------------- 401 : * appel a Ottosen 402 : *---------------------------------------------------------------------- 403 : * 404 : 404 : CALL OTTOSE(INPLAS,SIG0,NSTRS,EPSOTT,VAR0,XMAT,IVALMA,NCXMAT, 405 : & XCARB,ICARA,NVARI,SIGF,VARF,DEPFIS,MFR,KERRE,IB,IGAU) 406 : * 407 : *---------------------------------------------------------------------- 408 : * deformation corrigee pour Maxwell 409 : *---------------------------------------------------------------------- 410 : * 411 : DO II=1,NSTRS 412 : DEPSF(II) = DEPST(II) - DEPFIS(II) 413 : END DO 414 : *---------------------------------------------------------------------- 415 : * modele de Maxwell 416 : *---------------------------------------------------------------------- 417 : DO M=0,NCHAIN-1 417 : CALL CFICMO(TPS1,TPS2,M,WR12,wrk52,wrk53,NWA,EFIC,TR,E1,E2) 418 : IF (E1.EQ.0.) THEN 419 : COE = 1. 420 : ELSE 421 : COE = E2 / E1 422 : ENDIF 423 : FLUA=EXP(-TR*DT) 424 : IF (M.EQ.0) THEN 425 : SM0(1)=SM0(1)*COE*FLUA 426 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 427 : SM0(2)=SM0(2)*COE*FLUA 428 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 429 : SM0(3)=0.D0 430 : SM0(4)=SM0(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 431 : DO II=1,NSTRS 432 : SIGM(II)=SM0(II) 433 : END DO 434 : ELSE IF (M.EQ.1) THEN 435 : SM1(1)=SM1(1)*COE*FLUA 436 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 437 : SM1(2)=SM1(2)*COE*FLUA 438 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 439 : SM1(3)=0.D0 440 : SM1(4)=SM1(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 441 : DO II=1,NSTRS 442 : SIGM(II)=SIGM(II)+SM1(II) 443 : END DO 444 : ELSE IF (M.EQ.2) THEN 445 : SM2(1)=SM2(1)*COE*FLUA 446 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 447 : SM2(2)=SM2(2)*COE*FLUA 448 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 449 : SM2(3)=0.D0 450 : SM2(4)=SM2(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 451 : DO II=1,NSTRS 452 : SIGM(II)=SIGM(II)+SM2(II) 453 : END DO 454 : ELSE IF (M.EQ.3) THEN 455 : SM3(1)=SM3(1)*COE*FLUA 456 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 457 : SM3(2)=SM3(2)*COE*FLUA 458 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 459 : SM3(3)=0.D0 460 : SM3(4)=SM3(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 461 : DO II=1,NSTRS 462 : SIGM(II)=SIGM(II)+SM3(II) 463 : END DO 464 : ELSE IF (M.EQ.4) THEN 465 : SM4(1)=SM4(1)*COE*FLUA 466 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 467 : SM4(2)=SM4(2)*COE*FLUA 468 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 469 : SM4(3)=0.D0 470 : SM4(4)=SM4(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 471 : DO II=1,NSTRS 472 : SIGM(II)=SIGM(II)+SM4(II) 473 : END DO 474 : ELSE IF (M.EQ.5) THEN 475 : SM5(1)=SM5(1)*COE*FLUA 476 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 477 : SM5(2)=SM5(2)*COE*FLUA 478 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 479 : SM5(3)=0.D0 480 : SM5(4)=SM5(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 481 : DO II=1,NSTRS 482 : SIGM(II)=SIGM(II)+SM5(II) 483 : END DO 484 : ELSE IF (M.EQ.6) THEN 485 : SM6(1)=SM6(1)*COE*FLUA 486 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 487 : SM6(2)=SM6(2)*COE*FLUA 488 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 489 : SM6(3)=0.D0 490 : SM6(4)=SM6(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 491 : DO II=1,NSTRS 492 : SIGM(II)=SIGM(II)+SM6(II) 493 : END DO 494 : ELSE IF (M.EQ.7) THEN 495 : SM7(1)=SM7(1)*COE*FLUA 496 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 497 : SM7(2)=SM7(2)*COE*FLUA 498 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 499 : SM7(3)=0.D0 500 : SM7(4)=SM7(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 501 : DO II=1,NSTRS 502 : SIGM(II)=SIGM(II)+SM7(II) 503 : END DO 504 : ELSE IF (M.EQ.8) THEN 505 : SM8(1)=SM8(1)*COE*FLUA 506 : 1 +CONS*EFIC*(DEPSF(1)+ANU*DEPSF(2)) 507 : SM8(2)=SM8(2)*COE*FLUA 508 : 1 +CONS*EFIC*(DEPSF(2)+ANU*DEPSF(1)) 509 : SM8(3)=0.D0 510 : SM8(4)=SM8(4)*COE*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPSF(4) 511 : DO II=1,NSTRS 512 : SIGM(II)=SIGM(II)+SM8(II) 513 : END DO 514 : END IF 515 : END DO 516 : DEFP(1)=DEPSF(1)-((SIGM(1)-ANU*(SIGM(2)+SIGM(3)))/SE2 517 : 1 -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1) 518 : DEFP(2)=DEPSF(2)-((SIGM(2)-ANU*(SIGM(1)+SIGM(3)))/SE2 519 : 1 -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1) 520 : DEFP(3)=DEPSF(3)-((SIGM(3)-ANU*(SIGM(2)+SIGM(1)))/SE2 521 : 1 -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1) 522 : DEFP(4)=DEPSF(4)-(1+ANU)*(SIGM(4)/SE2-SIG0(4)/SE1) 523 : DO I=1,NSTRS 524 : EPINF(I)=EPIN0(I)+DEFP(I)+DEPFIS(I) 525 : END DO 526 : ENDIF 527 : 528 : *====================================================================== 529 : * QUELQUES IMPRESSIONS 530 : *====================================================================== 531 : * 532 : IF (IIMPI.EQ.103) THEN 533 : WRITE(*,66770) IB,IGAU 534 : WRITE(*,*) dt,tps2 535 : WRITE(*,66775) XMAT(1) 536 : WRITE(*,66775) (XMAT(I),I=1,NCXMAT) 537 : WRITE(*,66776) (XCARB(I),I=1,ICARA) 538 : WRITE(*,66773) (VAR0(I),I=1,NVARI) 539 : WRITE(*,66774) (DEPST(I),I=1,NSTRS) 540 : WRITE(*,66777) (EPSFLU(I),I=1,NSTRS) 541 : WRITE(*,66778) (EPSOTT(I),I=1,NSTRS) 542 : WRITE(*,66779) (DEPFIS(I),I=1,NSTRS) 543 : WRITE(*,66771) (DEPSF(I),I=1,NSTRS) 544 : WRITE(*,66794) (DEFP(I),I=1,NSTRS) 545 : WRITE(*,66784) (SIGM(I),I=1,NSTRS) 546 : WRITE(*,66772) (SIG0(I),I=1,NSTRS) 547 : WRITE(*,66783) (SIGF(I),I=1,NSTRS) 548 : WRITE(*,66795) (EPIN0(I),I=1,NSTRS) 549 : WRITE(*,66782) (EPINF(I),I=1,NSTRS) 550 : ENDIF 551 : 66770 FORMAT(////////2X,'ELEMENT ',I6,2X,'POINT ',I3//) 552 : 66771 FORMAT(2X,' DEPSF '/(6(1X,1PE12.5))) 553 : 66772 FORMAT(2X,' SIG0 '/(6(1X,1PE12.5))) 554 : 66773 FORMAT(2X,' VAR0 '/(6(1X,1PE12.5))) 555 : 66774 FORMAT(2X,' DEPST '/(6(1X,1PE12.5))) 556 : 66775 FORMAT(2X,' XMAT '/(6(1X,1PE12.5))) 557 : 66776 FORMAT(2X,' XCARB '/(6(1X,1PE12.5))) 558 : 66777 FORMAT(2X,' EPSFLU '/(6(1X,1PE12.5))) 559 : 66778 FORMAT(2X,' EPSOTT '/(6(1X,1PE12.5))) 560 : 66779 FORMAT(2X,' DEPFIS '/(6(1X,1PE12.5))) 561 : 66782 FORMAT(2X,' EPINF '/(6(1X,1PE12.5))) 562 : 66783 FORMAT(2X,' SIGF '/(6(1X,1PE12.5))) 563 : 66784 FORMAT(2X,' SIGM '/(6(1X,1PE12.5))) 564 : 66794 FORMAT(2X,' DEFP '/(6(1X,1PE12.5))) 565 : 66795 FORMAT(2X,' EPIN0 '/(6(1X,1PE12.5))) 566 : 567 : *---------------------------------------------------------------------- 568 : RETURN 569 : END 570 : 571 : 572 : 573 : 574 : 575 : 576 : 577 : 578 : 579 : 580 : 581 : 582 : 583 : 584 : 585 : 586 : 587 : 588 : 589 : 590 : 591 : 592 : 593 : 594 : 595 : 596 : 597 : 598 :
© Cast3M 2003 - All rights reserved.
Disclaimer