1 : C CURGCS SOURCE KICH 09/09/11 21:17:56 6505 2 : SUBROUTINE CURGCS(WRK52,WRK53,WRK54,WRK4,NSTRSS, 3 : 1 IFOUR,DDT,IB,IGAU,HCAR,IVIS,T0,TF) 4 : C--------------------------------------------------------------------- 5 : C PLASTICITE MODELE BETON 6 : C 7 : C ENTREES 8 : C SIG0(NSTRS) = CONTRAINTES INITIALES 9 : C NSTRS = NOMBRE DE CONTRAINTES 10 : C DEPST(NSTRS) = INCREMENT DE DEFORMATIONS TOTALES 11 : C VAR0(NVARI) = VARIABLES INTERNES DEBUT 12 : C VAR0( 1 ) = IFIS :Indicateur de fissuration (0 1 2) 13 : C VAR0( 2 ) = ANGL :Angle de fissuration 14 : C VAR0( 3 ) = IPLA :Indicateur d'etat en bicompres. (0 1 2 3) 15 : C VAR0( 4 ) = SIG1 :Force d écrouissage de traction 16 : C VAR0( 5 ) = SIG2 :Force d écrouissage de compression 17 : C VAR0( 6 ) = EPS1 :Variable d écrouissage de traction 18 : C VAR0( 7 ) = EPS2 :Variable d écrouissage de compression 19 : C VAR0( 8 ) = TDEF :Taux de déformation 20 : C VAR0( 9 ) = TCON :Taux de contrainte 21 : C VAR0( 10 ) = SIGP(1):Contrainte plastique en mode 22 : C viscoplastique 23 : C VAR0( 11 ) = SIGP(2): 24 : C VAR0( 12 ) = SIGP(3): 25 : C VAR0( 13 ) = SIGP(4): 26 : C VAR0( 14 ) = DPSTV1 :Variable d'écrouissage plastique 27 : C en mode viscoplastique 28 : C VAR0( 15 ) = DPSTV2 : 29 : C VAR0( 16 ) = SIGV1 :Force d ecrouissage de traction 30 : C en mode viscoplastique 31 : C VAR0( 17 ) = SIGV2 :Force d ecrouissage de compression 32 : C en mode viscoplastique 33 : C 34 : C XMAT(NCOMAT) = COMPOSANTES DE MATERIAU 35 : C NCOMAT = NOMBRE DE COMPOSANTES DE MATERIAU 36 : C SORTIES 37 : C SIGF(NSTRS) = CONTRAINTES FINALES 38 : C VARF(NVARI) = VARIABLES INTERNES FINALES 39 : C KERRE = 0 TOUT OK 40 : C--------------------------------------------------------------------- 41 : C 42 : C IFOUR INDICE DU TYPE DE PROBLEME 43 : C -2 CONTRAINTES PLANES 44 : C -1 DEFORMATIONS PLANES 45 : C 0 AXISYMETRIQUE 46 : C 1 SERIE DE FOURIER 47 : C 2 TRIDIMENSIONNEL 48 : C--------------------------------------------------------------------- 49 : C COMPOSANTES DE MATERIAU 50 : C===================================================================== 51 : C YOUN : Module d'Young 52 : C XNU : Coeficient de Poisson 53 : C RHO : Masse volumique (Facultatif) 54 : C ALPH : Coeficient de dilation thermique (Facultatif) 55 : C ALFA : Resis. tract. simple / resis. compr. simple 56 : C BETA : Resis. compr. biax. / resis. compr. simple (Modele de NADAI) 57 : C RB : Resis. compr. simple 58 : C GFC : Energie de rupt en compression 59 : C GFT : Energie de rupt en traction 60 : C ETA : parametre de viscosite 61 : C ICT : Choix de la courbe de traction (Calibrage) 62 : C ICC : Choix de la courbe de compression (Calibrage) 63 : C 64 : C IVIS : Modele visqueux ou non 65 : C 0 : non visqueux 66 : C 1 : viscoplastique 67 : C 2 : viscoelastoplastique 68 : C 3 : elastoplastique endommageable 69 : C 4 : elasto viscoplastique endommageable 70 : C 71 : C IMOD : Choix du modele Beton 72 : C 1 = Modele BETON_INSA ( Von MISES + Sigma Max en trac.) 73 : C Rheologie 2D 74 : C 2 = Modele BETON_INSA ( Von MISES + Sigma Max en trac.) 75 : C Rheologie 3D 76 : C 3 = Modele BETON_INSA ( Druck Prager + Sigma Max en trac.) 77 : C Rheologie 2D 78 : C 4 = Modele BETON_INSA ( Druck Prager + Sigma Max en trac.) 79 : C Rheologie 3D 80 : C ITER : nombre d'iterations internes admissibles 81 : C Pour modeles fluage 82 : C TZER = DATE MINI POUR APPLICATION DU PREMIER CHARGEMENT 83 : C PGDUR = CHOIX DE PROGRESSION DES PAS DE TEMPS DE DISCRETISATION DES COURBES DE RELAXATION 84 : C PGTZO = CHOIX DE PROGRESSION DES PAS DE TEMPS DE DISCRETISATION DES COURBES DE RELAXATION 85 : C ITYPE = CHOIX DU CODE DE CALCUL DE FLUAGE CEB-FIP, B3, EUROCODE2, ETC 86 : C TAU2 = SECOND TEMPS DE RELAXATION PERMETTANT LE CALCUL DE LA SUITE DES TEMPS DE RELAXATION 87 : C TAU1 = PREMIER TEMPS DE RELAXATION PERMETTANT LA REPONSE ELASTIQUE 88 : C DATSEC = DATE DE DEBUT DE SECHAGE 89 : C DATCUR = DATE DE FIN DE CURE DU BETON 90 : C NBRC = NOMBRE DE BRANCHE DE MAXWELL 91 : C DATCOU = DATE D APPLICATION DU CHARGEMENT 92 : C IFOR = 1 formulation solide 93 : C = 2 liquide 94 : C--------------------------------------------------------------------- 95 : C 96 : IMPLICIT INTEGER(I-N) 97 : IMPLICIT REAL*8(A-H,O-Z) 98 : -INC DECHE 99 : CHARACTER*40 FMT,TITRE 100 : DIMENSION SIR(9,4),EPST(4),SIGP(4) 101 : DIMENSION D(6,6),D1(6,6),STRN(4) 102 : DIMENSION SIGR(4),VART(200),VV(36),SIGMF(4) 103 : DIMENSION SIRL(8,4) 104 : C 105 : SEGMENT WRK4 106 : REAL*8 XE(3,NBNN) 107 : ENDSEGMENT 108 : 109 : SEGMENT BETJEF 110 : REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF, 111 : & TCON,DPSTF1,DPSTF2,TETA,PDT,TP00 112 : INTEGER ICT,ICC,IMOD,IVISS,ITER, 113 : & ISIM,IBB,IGAU1,IZON 114 : ENDSEGMENT 115 : SEGMENT VISCO 116 : REAL*8 DPSTV1,DPSTV2,SIGV1,SIGV2,ENDV 117 : ENDSEGMENT 118 : SEGMENT BETFLU 119 : REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2, 120 : & TP0,TZER 121 : INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR 122 : ENDSEGMENT 123 : SEGMENT NECH0 124 : REAL*8 DT,DC,ALFG,S0,ENDO 125 : ENDSEGMENT 126 : SEGMENT NECH1 127 : REAL*8 DLMT 128 : ENDSEGMENT 129 : SEGMENT NECH2 130 : REAL*8 ATR,GTR,ALPH0 131 : ENDSEGMENT 132 : SEGMENT NECH3 133 : REAL*8 RBT,ALFAT,YOUNT,GFCT,GFTT,ALPH 134 : ENDSEGMENT 135 : C 136 : C NCOMAT = nmatt 137 : C NBNN=XE(/2) 138 : NVARI=VAR0(/1) 139 : SEGINI BETJEF 140 : SEGINI VISCO 141 : SEGINI BETFLU 142 : SEGINI NECH0 143 : SEGINI NECH1 144 : SEGINI NECH2 145 : SEGINI NECH3 146 : PDT = DDT 147 : TP0=TEMP0 148 : TP00=TEMP0 149 : CAR = HCAR 149 : CALL ZERO(SIGP,4,1) 149 : CALL ZERO(SIR,9,4) 149 : CALL ZERO(SIRL,8,4) 150 : C 151 : C TEST DE CONSISTANCE DES DONNEES 152 : C 153 : YOUNT = 0.D0 154 : YOUN = 0.D0 155 : XNU = 0.D0 156 : RHO = 0.D0 157 : ALPH = 0.D0 158 : ALFAT = 0.D0 159 : ALFA = 0.D0 160 : BETA = 0.D0 161 : RBT = 0.D0 162 : RB = 0.D0 163 : GFCT = 0.D0 164 : GFTT = 0.D0 165 : GFC = 0.D0 166 : GFT = 0.D0 167 : ETA = 0.D0 168 : ITER = 0 169 : IMOD = 0 170 : ICT = 0 171 : ICC = 0 172 : DC =0.D0 173 : DT =0.D0 174 : IASO =0 175 : ICOU =0 176 : ENDO = 0.D0 177 : ENDV =0.D0 178 : C 179 : AA=0.D0 180 : BB=0.D0 181 : DK1=0.D0 182 : DK2=0.D0 183 : CNC1=0.D0 184 : CNC1=0.D0 185 : CNT1=0.D0 186 : CNT2=0.D0 187 : C 188 : YOUNT = XMAT( 1) 189 : XNU = XMAT( 2) 190 : RHO = XMAT( 3) 191 : ALPH = XMAT( 4) 192 : ALFAT = XMAT( 5) 193 : BETA = XMAT( 6) 194 : RBT = XMAT( 7) 195 : GFCT = XMAT( 8) 196 : GFTT = XMAT( 9) 197 : ITER = INT(SNGL(XMAT(10))) 198 : IMOD = INT(SNGL(XMAT(11))) 199 : IF (IVIS.EQ.0) THEN 200 : ICT = INT(SNGL(XMAT(12))) 201 : ICC = INT(SNGL(XMAT(13))) 202 : ENDIF 203 : IF (IVIS.EQ.1) THEN 204 : ICT = INT(SNGL(XMAT(12))) 205 : ICC = INT(SNGL(XMAT(13))) 206 : ETA = XMAT(14) 207 : ISIM = 1 208 : ENDIF 209 : IF (IVIS.EQ.2) THEN 210 : ICT = INT(SNGL(XMAT(12))) 211 : ICC = INT(SNGL(XMAT(13))) 212 : DATCOU = XMAT(14) 213 : NBRC = INT(SNGL(XMAT(15))) 214 : IFOR = INT(SNGL(XMAT(16))) 215 : ITYPE = INT(SNGL(XMAT(17))) 216 : DATCUR = XMAT(18) 217 : DATSEC = XMAT(19) 218 : E28= YOUN 219 : ENDIF 220 : IF (IVIS.EQ.3) THEN 221 : DC = XMAT(12) 222 : DT = XMAT(13) 223 : IASO = INT(SNGL(XMAT(14))) 224 : ICOU = INT(SNGL(XMAT(15))) 225 : ENDIF 226 : IF (IVIS.EQ.4) THEN 227 : DC = XMAT(12) 228 : DT = XMAT(13) 229 : IASO = INT(SNGL(XMAT(14))) 230 : ICOU = INT(SNGL(XMAT(15))) 231 : ETA = XMAT(16) 232 : ENDIF 233 : C 234 : C--------------------------------------------------------------------- 235 : C INITIALISATION DES VARIABLES DU MODELE D'ENDOMMAGEMENT 236 : C--------------------------------------------------------------------- 237 : C A: CALCUL ASSOCIE OU NON 238 : C 239 : IF(IASO.EQ.0) THEN 240 : ALFG = (BETA-1.D0)/(2*BETA-1.D0) 241 : ELSE 242 : ALFG = 0.2 243 : ENDIF 244 : C--------------------------------------------------------------------- 245 : C B: CALCUL COUPLE OU NON 246 : C 247 : IF(ICOU.EQ.0) THEN 248 : ATR=0.D0 249 : GTR=0.D0 250 : ELSE 251 : ATR=2.075D0 252 : GTR=0.235D0 253 : ENDIF 254 : C--------------------------------------------------------------------- 255 : C 256 : DLMT=0.9999D0 257 : C S0 = 0.25D0 258 : IF(IVIS.EQ.4)THEN 259 : C attention en comportement visqueux refermeture pose probleme 260 : S0 = 1.D0 261 : ELSE 262 : S0 = 0.25D0 263 : C S0 = 1.D0 264 : ENDIF 265 : C 266 : C--------------------------------------------------------------------- 267 : C CAS DU COUPLAGE THERMIQUE-MECANIQUE 268 : C--------------------------------------------------------------------- 269 : C IF IASO=0 Utilisation des lois f(T) locales 270 : C IF IASO=1 Utilisation des lois du fichier de donnée 271 : C 272 : ITHHER = 1 273 : IF(ITHHER.NE.0)THEN 273 : CALL EVOLTH(T0,TF,BETJEF,NECH0,NECH3) 274 : ELSE 275 : CC write(*,*)'ITHHER=',ITHHER 276 : RB = RBT 277 : ALFA = ALFAT 278 : YOUN = YOUNT 279 : GFC = GFCT 280 : GFT = GFTT 281 : ENDIF 282 : C 283 : C 284 : C--------------------------------------------------------------------- 285 : C 286 : C CORRESPONDANCE DES VARIABLES POUR URGCST 287 : C 287 : CALL ZERO(D,6,6) 287 : CALL ZERO(D1,6,6) 287 : CALL ZERO(VART,NVARI,1) 287 : CALL ZERO(SIGR,4,1) 288 : C 289 : IVISS = IVIS 290 : IBB = IB 291 : IGAU1 = IGAU 292 : EPAIST=1.D0 293 : NSTRS=NSTRSS 294 : IFOU=IFOUR 295 : * 296 : DO 1 I=1,NSTRS 297 : SIGR(I)=SIG0(I)/EPAIST 298 : * WRITE(*,*)'SIG0',I,'=',SIG0(I) 299 : STRN(I)=DEPST(I) 300 : * WRITE(*,*)'DEPST',I,'=',DEPST(I) 301 : 1 CONTINUE 302 : * 303 : IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN 304 : NSTRS=3 305 : IFOU=-2 306 : SIGR(3)=SIG0(4) 307 : STRN(3)=DEPST(4) 308 : STRN(4)=DEPST(3) 309 : ENDIF 310 : * 311 : IF((NSTRS.EQ.4).AND.IFOUR.EQ.-2.AND. 312 : * (IMOD.NE.2.AND.IMOD.NE.4)) THEN 313 : NSTRS=3 314 : SIGR(3)=SIG0(4) 315 : STRN(3)=DEPST(4) 316 : STRN(4)=DEPST(3) 317 : ENDIF 318 : * 319 : DO 22 II=1,NVARI 320 : VART(II)=VAR0(II) 321 : 22 CONTINUE 322 : C------------------------------------------------------ 323 : C INITIALISATION DES VARIABLES INTERNES 324 : C------------------------------------------------------ 325 : C A: MODELE PLASTIQUE 326 : C 327 : IFISU = INT(SNGL(VART(1))) 328 : ANGL = VART(2) 329 : IPLA = INT(SNGL(VART(3))) 330 : SIG1= VART(4) 331 : SIG2= VART(5) 332 : DPSTF1 = VART(6) 333 : DPSTF2 = VART(7) 334 : C 335 : C------------------------------------------------------ 336 : C B: MODELE VISCOPLASTIQUE CLASSIQUE 337 : C 338 : IF (IVIS.EQ.1) THEN 339 : TDEF = VART( 8) 340 : TCON = VART( 9) 341 : SIGP(1) = VART(10) 342 : SIGP(2) = VART(11) 343 : SIGP(3) = VART(12) 344 : SIGP(4) = VART(13) 345 : DPSTV1 = VART(14) 346 : DPSTV2 = VART(15) 347 : SIGV1 = VART(16) 348 : SIGV2 = VART(17) 349 : ENDIF 350 : C 351 : C------------------------------------------------------ 352 : C C: MODELE VISCO-ELASTO PLASTIQUE 353 : C 354 : IF (IVIS.EQ.2.AND.IFOR.EQ.1) THEN 355 : C modele visco elasto plastique solide de MAXWELL 356 : C------------------------------------------------- 357 : MC = NBRC + 1 358 : DO 101 I1 = 1,MC 359 : DO 102 J1 = 1,NSTRS 360 : NV = NSTRS * (I1 - 1) + J1 361 : SIR(I1,J1) = VART(7 + NV) 362 : 102 CONTINUE 363 : 101 CONTINUE 364 : C 365 : ELSE IF (IVIS.EQ.2.AND.IFOR.EQ.2) THEN 366 : C modele visco elasto plastique liquide de MAXWELL 367 : C------------------------------------------------- 368 : MC = NBRC + 1 369 : DO 131 I2 = 1,MC 370 : DO 132 J2 = 1,NSTRS 371 : NV = NSTRS * (I2 - 1) + J2 372 : SIRL(I2,J2) = VART(7 + NV) 373 : 132 CONTINUE 374 : 131 CONTINUE 375 : C 376 : ENDIF 377 : C 378 : C---------------------------------------------- 379 : C D: MODELE PLASTIQUE ENDOMMAGEABLE 380 : C 381 : IF (IVIS.EQ.3) THEN 382 : ENDO = VART( 8) 383 : DECO = VART( 9) 384 : ENDIF 385 : C 386 : C---------------------------------------------- 387 : C D: MODELE VISCOPLASTIQUE ENDOMMAGEABLE 388 : C 389 : IF (IVIS.EQ.4) THEN 390 : TDEF = VART( 8) 391 : TCON = VART( 9) 392 : SIGP(1) = VART(10) 393 : SIGP(2) = VART(11) 394 : SIGP(3) = VART(12) 395 : SIGP(4) = VART(13) 396 : DPSTV1 = VART(14) 397 : DPSTV2 = VART(15) 398 : SIGV1 = VART(16) 399 : SIGV2 = VART(17) 400 : ENDO = VART(18) 401 : ENDV = VART(19) 402 : ENDIF 403 : C ------------------------------------- 404 : TETA=ANGL 405 : C--------------------------------------------------------------------- 406 : C 407 : GO TO (10,10,30,40),NSTRS 408 : 10 CONTINUE 409 : KERRE=437 410 : WRITE(*,*) '!! ATTENTION DANS BETON NSTRS=',NSTRS 411 : STOP 412 : GO TO 1000 413 : C 414 : 40 CONTINUE 415 : C 416 : 30 CONTINUE 417 : C 418 : C--------------------------------------------------------------------- 419 : C 420 : IF(IMOD.GE.1.AND.IMOD.LE.4) THEN 421 : * *-----------------------* 422 : * * MODELE URGC ST * 423 : * *-----------------------* 424 : * 424 : CALL BONE(SIGR,SIGMF,STRN,IPLA,IFISU,SIG1,SIG2 425 : A ,NSTRS,D,D1,IFOU,SIGP,EPST,SIR,SIRL, 426 : B ENDO,ITHHER,T0,TF,BETJEF,VISCO,BETFLU, 427 : C NECH0,NECH1,NECH2,NECH3) 428 : C 429 : ELSE 430 : WRITE(*,*) '!! ATTENTION CE MODELE N EXISTE PAS (URGCST)' 431 : STOP 432 : RETURN 433 : ENDIF 434 : C--------------------------------------------------------------------- 435 : C 436 : C 437 : IF((NSTRS.EQ.3).AND.IFOU.EQ.-2.AND. 438 : * (IMOD.NE.2.AND.IMOD.NE.4)) THEN 439 : SIGMF(4)=SIGMF(3) 440 : SIGMF(3)=0.D0 441 : ENDIF 442 : C 443 : IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.3) THEN 444 : SIGMF(3)=0.D0 445 : SIGMF(4)=SIGMF(3) 446 : ENDIF 447 : C 448 : DO 2 I=1,NSTRSS 449 : SIGF(I)=SIGMF(I) 450 : C WRITE(*,*)'SIGMF(',I,')=',SIGMF(I) 451 : 2 CONTINUE 452 : C 453 : C--------------------------------------------------------------------- 454 : C 455 : C--------------------------------------------------------------------- 456 : C ACTUALISATION DES VARIABLES INTERNES 457 : C--------------------------------------------------------------------- 458 : C A: MODELE PLASTIQUE 459 : C 460 : VART( 1)=FLOAT(IFISU) 461 : VART( 2)=TETA 462 : VART( 3)=FLOAT(IPLA) 463 : VART( 4)=SIG1 464 : VART( 5)=SIG2 465 : VART( 6)=DPSTF1 466 : VART( 7)=DPSTF2 467 : K1 = 7 468 : C 469 : C--------------------------------------------------------------------- 470 : C B: MODELE VISCOPLASTIQUE CLASSIQUE 471 : C 472 : IF (IVIS.EQ.1) THEN 473 : VART(8)=TDEF 474 : VART(9)=TCON 475 : VART(10)=SIGP(1) 476 : VART(11)=SIGP(2) 477 : VART(12)=SIGP(3) 478 : VART(13)=SIGP(4) 479 : VART(14)=DPSTV1 480 : VART(15)=DPSTV2 481 : VART(16)=SIGV1 482 : VART(17)=SIGV2 483 : K1 = 17 484 : ENDIF 485 : C--------------------------------------------------------------------- 486 : C C: MODELE VISCO-ELASTO PLASTIQUE 487 : C 488 : IF (IVIS.EQ.2.AND.IFOR.EQ.1) THEN 489 : MC = NBRC + 1 490 : DO 103 I1 = 1,MC 491 : DO 104 J1 = 1,NSTRS 492 : NV = NSTRS * (I1 - 1) + J1 493 : VART(7 + NV) = SIR(I1,J1) 494 : 104 CONTINUE 495 : 103 CONTINUE 496 : C DO 123 I1 = 1,MC 497 : C DO 124 J1 = 1,MC 498 : C K = MC * (I1 - 1) + J1 499 : C VART(43 + K) = CODU(I1,J1) 500 : C 124 CONTINUE 501 : C 123 CONTINUE 502 : K1=43 503 : ENDIF 504 : C---------------------------------------- 505 : IF (IVIS.EQ.2.AND.IFOR.EQ.2) THEN 506 : DO 108 I1 = 1,NBRC 507 : DO 109 J1 = 1,NSTRS 508 : NV = NSTRS * (I1 - 1) + J1 509 : VART(7 + NV) = SIRL(I1,J1) 510 : 109 CONTINUE 511 : 108 CONTINUE 512 : C DO 143 I1 = 1,NBRC 513 : C DO 144 J1 = 1,NBRC 514 : C K = NBRC * (I1-1) + J1 515 : C VART(43 + K) = CODL(I1,J1) 516 : C 144 CONTINUE 517 : C 143 CONTINUE 518 : K1=43 519 : ENDIF 520 : C 521 : C--------------------------------------------------------------------- 522 : C D: MODELE PLASTIQUE ENDOMMAGEABLE 523 : C 524 : IF (IVIS.EQ.3) THEN 525 : VART(8) = ENDO 526 : VART(9) = MAX(VAR0(9),ENDO) 527 : K1= 9 528 : ENDIF 529 : C 530 : C--------------------------------------------------------------------- 531 : C B: MODELE VISCOPLASTIQUE ENDOMMAGEABLE 532 : C 533 : IF (IVIS.EQ.4) THEN 534 : VART(8)=TDEF 535 : VART(9)=TCON 536 : VART(10)=SIGP(1) 537 : VART(11)=SIGP(2) 538 : VART(12)=SIGP(3) 539 : VART(13)=SIGP(4) 540 : VART(14)=DPSTV1 541 : VART(15)=DPSTV2 542 : VART(16)=SIGV1 543 : VART(17)=SIGV2 544 : VART(18)=ENDO 545 : VART(19)=ENDV 546 : K1 = 19 547 : ENDIF 548 : C 549 : GO TO (11,11,31,41),NSTRS 550 : 11 CONTINUE 551 : KERRE=437 552 : WRITE(*,*) '!! ATTENTION DANS BETON NSTRS=',NSTRS 553 : STOP 554 : GO TO 1000 555 : C 556 : 41 CONTINUE 557 : C 558 : 31 CONTINUE 559 : C 560 : C CONTRAINTES OU DEFORMATIONS PLANES 561 : C SANS CISAILLEMENTS TRANSVERSAL 562 : C 563 : L1=0 564 : DO 66 I1 = 1,NSTRS 565 : DO 67 J1 = 1,NSTRS 566 : L1=L1+1 567 : VV(L1) = D(J1,I1) 568 : 67 CONTINUE 569 : 66 CONTINUE 570 : * 571 : L1=0 572 : DO 62 I1 = 1,NSTRS 573 : DO 63 J1 = 1,NSTRS 574 : K=K1+(I1-1)*NSTRS+J1 575 : L1=L1+1 576 : VART(K)=VV(L1) 577 : 63 CONTINUE 578 : 62 CONTINUE 579 : C 580 : DO 23 II=1,NVARI 581 : VARF(II)=VART(II) 582 : 23 CONTINUE 583 : C 584 : C--------------------------------------------------------------------- 585 : C 586 : 1000 CONTINUE 587 : SEGSUP BETJEF 588 : SEGSUP VISCO 589 : SEGSUP BETFLU 590 : SEGSUP NECH0 591 : SEGSUP NECH1 592 : SEGSUP NECH2 593 : SEGSUP NECH3 594 : C 595 : RETURN 596 : END 597 : 598 : 599 : 600 : 601 : 602 : 603 : 604 : 605 : 606 :
© Cast3M 2003 - All rights reserved.
Disclaimer