Esope source files

   1 : C CMAXOB    SOURCE    KICH      09/09/11    21:16:14     6505           
   2 :       SUBROUTINE CMAXOB(wrk52,wrk53,wrk54,WR12,NWA,NCHAIN,EPSFLU)
   3 : *
   4 : *      MODELE MAXOTT  : RECUPERATION DES VARIABLES INTERNES
   5 : *                           APRES ECOULEMENT
   6 : *
   7 :       IMPLICIT INTEGER(I-N)
   8 :       IMPLICIT REAL*8(A-H,O-Z)
   9 : *
  10 : -INC CCOPTIO
  11 : -INC SMEVOLL
  12 : -INC SMLREEL
  13 : -INC CCHAMP
  14 : -INC DECHE
  15 : 
  16 :       SEGMENT WR12
  17 :         REAL*8  EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  18 :         REAL*8  EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  19 :         REAL*8  EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  20 :         REAL*8  SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  21 :         REAL*8  SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  22 :         REAL*8  SM8(NSTRS)
  23 :       ENDSEGMENT
  24 : 
  25 :       DIMENSION NWA(9)
  26 :       DIMENSION EPSFLU(8)
  27 : 
  28 : *
  29 : *      ON CREE LES PROGS ASSOCIES AUX VARIABLES INTERNES
  30 : *
  31 : * Formulation tridimensionnelle massive
  32 : *
  33 :              IF(IFOUR.EQ.2.AND.MFR.EQ.1) THEN
  34 :                MXVREE = 22
  35 : *
  36 : * Formulation deformations planes / axisymetrique
  37 : *
  38 :              ELSE IF(IFOUR.EQ.-1.OR.IFOUR.EQ.-3
  39 :      &         .OR.IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  40 :                MXVREE = 17
  41 : *
  42 : * Formulation contraintes planes / 3D coques minces
  43 : *
  44 :              ELSE
  45 :                MXVREE = 14
  46 :              ENDIF
  47 :              DO JC=0,NCHAIN
  48 :                  JG=NSTRS
  49 :                  SEGINI MLREEL
  50 :                  VARF(JC+MXVREE)=DBLE(MLREEL)
  51 :                  IF(JC.EQ.0) THEN
  52 :                     DO JD=1,NSTRS
  53 :                        PROG(JD)=EPSFLU(JD)
  54 :                     END DO
  55 :                  ELSE IF(JC.EQ.1) THEN
  56 :                     DO JD=1,NSTRS
  57 :                        PROG(JD)=SM0(JD)
  58 :                     END DO
  59 :                  ELSE IF(JC.EQ.2) THEN
  60 :                     DO JD=1,NSTRS
  61 :                        PROG(JD)=SM1(JD)
  62 :                     END DO
  63 :                  ELSE IF(JC.EQ.3) THEN
  64 :                     DO JD=1,NSTRS
  65 :                        PROG(JD)=SM2(JD)
  66 :                     END DO
  67 :                  ELSE IF(JC.EQ.4) THEN
  68 :                     DO JD=1,NSTRS
  69 :                        PROG(JD)=SM3(JD)
  70 :                     END DO
  71 :                  ELSE IF(JC.EQ.5) THEN
  72 :                     DO JD=1,NSTRS
  73 :                        PROG(JD)=SM4(JD)
  74 :                     END DO
  75 :                  ELSE IF(JC.EQ.6) THEN
  76 :                     DO JD=1,NSTRS
  77 :                        PROG(JD)=SM5(JD)
  78 :                     END DO
  79 :                  ELSE IF(JC.EQ.7) THEN
  80 :                     DO JD=1,NSTRS
  81 :                        PROG(JD)=SM6(JD)
  82 :                     END DO
  83 :                  ELSE IF(JC.EQ.8) THEN
  84 :                     DO JD=1,NSTRS
  85 :                        PROG(JD)=SM7(JD)
  86 :                     END DO
  87 :                  ELSE IF(JC.EQ.9) THEN
  88 :                     DO JD=1,NSTRS
  89 :                        PROG(JD)=SM8(JD)
  90 :                     END DO
  91 :                  ENDIF
  92 :                  SEGDES MLREEL
  93 :             
  94 :              END DO
  95 :              DO JC=NCHAIN+1,9
  96 :                  JG=0
  97 :                  SEGINI MLREEL
  98 :                  VARF(JC+MXVREE)=DBLE(MLREEL)
  99 :                  SEGDES MLREEL
 100 :              END DO
 101 : *
 102 :       RETURN
 103 :       END
 104 : 
 105 :  
 106 :  
 107 :  
 108 :  
 109 :  
 110 :  
 111 :  
 112 :  
 113 :  
 114 :  
 115 :  
 116 :  
 117 :  
 118 :  
 119 :  
 120 :  
 121 :  
 122 : 
 123 : 
 124 : 
 125 : 
 126 :  
 127 :  
 128 :  
 129 :  
 130 :  
 131 :  
 132 :  
 133 :  
 134 :