Esope source files

   1 : C ELOCRAK1  SOURCE    CHAT      12/08/06    21:15:21     7479           
   2 :       SUBROUTINE ELOCRAK1(WRK52,WRK53,WRK54,NVARI,Iecou)
   3 : C
   4 : C====&===1=========2=========3=========4=========5=========6=========7==
   5 : C Commentaires :   Subroutine permettant de mettre en oeuvre le 
   6 : C                          modele ELOCRACK (CEA) pour representer
   7 : C                  le comportement cyclique du bÃ©ton sous                 
   8 : C                  sollicitations sismiques 
   9 : C
  10 : C Traits       : - 
  11 : C                - 
  12 : C                - 
  13 : C                - 
  14 : C                - 
  15 : C
  16 : C Auteur       : B. Richard (Dr - Ing.) - CEA/DEN/DANS/DM2S/SEMT/EMSI
  17 : C====&===1=========2=========3=========4=========5=========6=========7==
  18 : C
  19 : C----DECLARATION GENERALES----------------------------------------------
  20 : C
  21 :       IMPLICIT INTEGER(I-N)
  22 :       IMPLICIT REAL*8(A-H,O-Z)
  23 : C
  24 : C----LES INCLUDES-------------------------------------------------------
  25 : C      
  26 : -INC CCOPTIO
  27 : -INC DECHE
  28 : C
  29 : C----DECLARATION PARAMETRES---------------------------------------------
  30 : C 
  31 :        SEGMENT IECOU
  32 :        INTEGER icow1,icow2,icow3,icow4,icow5,icow6,icow7,
  33 :      1 icow8,icow9,icow10,icow11,icow12,icow13,icow14,icow15,icow16,
  34 :      2 icow17,icow18,icow19,icow20,icow21,icow22,icow23,icow24,
  35 :      3 icow25,icow26,icow27,icow28,icow29,icow30,ICARA,
  36 :      4 icow32,icow33,NSTRS1,MFR1,icow36,icow37,icow38,
  37 :      5 icow39,icow40,icow41,icow42,icow43,icow44
  38 :      
  39 :        INTEGER icow45,icow46,icow47,icow48,icow49,icow50,
  40 :      1 icow51,icow52,icow53,icow54,icow55,icow56
  41 :      2 icow57,icow58
  42 :        ENDSEGMENT
  43 : 
  44 :       SEGMENT WRKK2
  45 :          REAL*8 EPSILI(NSTRSV)
  46 :       END SEGMENT
  47 : 
  48 :       SEGMENT WRK6
  49 :         REAL*8 SIG0S(NSTRS),DEPSTS(NSTRS)
  50 :       END SEGMENT
  51 : 
  52 :       INTEGER NVARI,KCAS,IRTD,ISTRS
  53 :       
  54 :       REAL*8 PREC,EPAI,FAC,AUX,AUX1,AUX2,YOUN,XNU,UN
  55 :       
  56 :       PARAMETER (UN=1.D0)
  57 : 
  58 :       KERRE=0
  59 : C
  60 : C----MISE EN DONNEES----------------------------------------------------
  61 : C
  62 :       YOUN = XMAT(1)
  63 :       XNU  = XMAT(2)
  64 : C
  65 : C----CAS DE LA FORMULATION COQUES MINCES--------------------------------
  66 : C
  67 :       IF (MFR1 .EQ. 9) THEN
  68 :          
  69 :          EPAI=XCARB(1)
  70 :          SEGINI WRK6
  71 :          
  72 :          DO 101 ISTRS=1,NSTRS1
  73 :             SIG0S(ISTRS)=SIG0(ISTRS)
  74 :             DEPSTS(ISTRS)=DEPST(ISTRS)
  75 :   101    CONTINUE
  76 : 
  77 :          NSTRSV=4
  78 :          IFOUR2=-2
  79 :          JDIM =2
  80 : 
  81 :          DO 102 ISTRS=1,2
  82 :             SIG0(ISTRS)=SIG0S(ISTRS)/EPAI
  83 :   102    CONTINUE
  84 : 
  85 :          DEPST(3)=0.D0
  86 :          DEPST(4)=DEPSTS(3)
  87 : 
  88 :          SIG0(3)=0.D0
  89 :          SIG0(4)=SIG0S(3)/EPAI
  90 : C
  91 : C----CAS DE LA FORMULATION MASSIVE--------------------------------------
  92 : C
  93 :       ELSE IF (MFR1 .EQ. 1) THEN
  94 :       
  95 :          NSTRSV=NSTRS1
  96 :          IFOUR2=IFOUR
  97 :          JDIM = IDIM
  98 :       
  99 :       ELSE
 100 :       
 101 :          PRINT*,'MFR1=',MFR1
 102 :          KERRE=57
 103 :          RETURN
 104 :       
 105 :       END IF
 106 : 
 107 :       SEGINI WRKK2
 108 : 
 109 :       CMATE = 'ISOTROPE'
 110 :       KCAS=1
 111 :       
 111 :       CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,KCAS,DDHOOK,IRTD)
 112 : 
 113 :       IF ( IRTD .EQ. 1) THEN
 114 :          
 114 :          CALL MATVE1 (DDHOOK,DEPST,NSTRSV,NSTRSV,DSIGT,1)
 115 :          
 116 :          PREC=1.D-08
 117 :          
 117 :          CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRSV,2,DDHOOK,IRTD) 
 118 :          
 118 :          CALL INVALM(DDHOOK,NSTRSV,NSTRSV,IRTD,PREC)
 119 :          
 120 :          IF (IRTD.EQ.0)THEN
 121 : 
 121 :          CALL MATVE1 (DDHOOK,SIG0,NSTRSV,NSTRSV,EPSILI,1)
 122 : 
 123 :          icarbi=icara
 124 : C
 125 : C----CALCUL DES CONTRAINTES DE CAUCHY DANS LE CAS 3D--------------------
 126 : C         
 126 :          CALL  ELOCRAK2(WRK52,WRK53,WRK54,WRKK2,NSTRSV,NVARI,
 127 :      1                  ICARbi,JDIM,IFOUR2)
 128 :          
 129 :          icara=icarbi
 130 : C
 131 : C----CALCUL DES CONTRAINTES VRAIES DANS LE CAS DES COQUES MINCES--------
 132 : C          
 133 :          IF (MFR1 .EQ. 9) THEN
 134 :                               
 135 :             FAC=(EPAI**3)/12.D0
 136 :             AUX =FAC*YOUN/(1.D0-XNU*XNU)
 137 :             AUX1=FAC*YOUN*.5D0/(1.D0+XNU)
 138 :             AUX2=EPAI*YOUN*.5D0/(1.d0+XNU)/1.2d0
 139 : 
 140 :             DO 103 ISTRS=1,NSTRS1
 141 :                SIG0(ISTRS)=SIG0S(ISTRS)
 142 :                DEPST(ISTRS)=DEPSTS(ISTRS)
 143 :  103        CONTINUE
 144 : 
 145 :             DO 104 ISTRS=1,2
 146 :                SIGF(ISTRS)=SIGF(ISTRS)*EPAI
 147 :  104        CONTINUE
 148 : 
 149 :             SIGF(3)=SIGF(4)*EPAI
 150 :             SIGF(4)=SIG0(4)+AUX*(DEPST(4)+XNU*DEPST(5))
 151 :             SIGF(5)=SIG0(5)+AUX*(DEPST(5)+XNU*DEPST(4))
 152 :             SIGF(6)=SIG0(6)+AUX1*DEPST(6)
 153 :             SIGF(7)=SIG0(7)+AUX2*DEPST(7)
 154 :             SIGF(8)=SIG0(8)+AUX2*DEPST(8)
 155 :             SEGSUP WRK6
 156 :          END IF
 157 :          
 158 :          ELSE
 159 :             PRINT*,'ERREUR DANS INVALM'
 160 :             KERRE=56
 161 :          END IF
 162 :       ELSE
 163 :          PRINT*,'ERREUR DANS DOHMAS'
 164 :          KERRE=56
 165 :       END IF
 166 :       SEGSUP WRKK2
 167 :       RETURN
 168 :       END
 169 :  
 170 :  
 171 : 
 172 :  
 173 :  
 174 :  
 175 :  
 176 :  
 177 :  
 178 :  
 179 :