Sources Esope | Cast3M

cmaxob
C CMAXOB    SOURCE    PV        17/12/08    21:15:55     9660           
      SUBROUTINE CMAXOB(wrk52,wrk53,wrk54,WR12,NWA,NCHAIN,EPSFLU)
*
*      MODELE MAXOTT  : RECUPERATION DES VARIABLES INTERNES
*                           APRES ECOULEMENT
*
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
*
 
-INC PPARAM
-INC CCOPTIO
-INC SMEVOLL
-INC SMLREEL
-INC CCHAMP
-INC DECHE
 
      SEGMENT WR12
        REAL*8  EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
        REAL*8  EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
        REAL*8  EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
        REAL*8  SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
        REAL*8  SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
        REAL*8  SM8(NSTRS)
      ENDSEGMENT
 
      DIMENSION NWA(9)
      DIMENSION EPSFLU(8)
 
*
*      ON CREE LES PROGS ASSOCIES AUX VARIABLES INTERNES
*
* Formulation tridimensionnelle massive
*
             IF(IFOUR.EQ.2.AND.MFR.EQ.1) THEN
               MXVREE = 22
*
* Formulation deformations planes / axisymetrique
*
             ELSE IF(IFOUR.EQ.-1.OR.IFOUR.EQ.-3
     &         .OR.IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
               MXVREE = 17
*
* Formulation contraintes planes / 3D coques minces
*
             ELSE
               MXVREE = 14
             ENDIF
             DO JC=0,NCHAIN
                 JG=NSTRS
                 SEGINI MLREEL
                 VARF(JC+MXVREE)=DBLE(MLREEL)
                 IF(JC.EQ.0) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=EPSFLU(JD)
                    END DO
                 ELSE IF(JC.EQ.1) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM0(JD)
                    END DO
                 ELSE IF(JC.EQ.2) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM1(JD)
                    END DO
                 ELSE IF(JC.EQ.3) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM2(JD)
                    END DO
                 ELSE IF(JC.EQ.4) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM3(JD)
                    END DO
                 ELSE IF(JC.EQ.5) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM4(JD)
                    END DO
                 ELSE IF(JC.EQ.6) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM5(JD)
                    END DO
                 ELSE IF(JC.EQ.7) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM6(JD)
                    END DO
                 ELSE IF(JC.EQ.8) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM7(JD)
                    END DO
                 ELSE IF(JC.EQ.9) THEN
                    DO JD=1,NSTRS
                       PROG(JD)=SM8(JD)
                    END DO
                 ENDIF
                 SEGDES MLREEL
 
             END DO
             DO JC=NCHAIN+1,9
                 JG=0
                 SEGINI MLREEL
                 VARF(JC+MXVREE)=DBLE(MLREEL)
                 SEGDES MLREEL
             END DO
*
      RETURN
      END