1 : C ELOCRAK2 SOURCE CHAT 12/08/06 21:15:22 7479 2 : SUBROUTINE ELOCRAK2 (WRK52,WRK53,WRK54,WRKK2,NSTRS1,NVARI, 3 : 1 ICARA,JDIM,IFOUR2) 4 : C 5 : C====&===1=========2=========3=========4=========5=========6=========7== 6 : C Commentaires : Subroutine permettant de mettre en oeuvre le 7 : C modele ELOCRACK (CEA) pour representer 8 : C le comportement cyclique du béton sous 9 : C sollicitations sismiques 10 : C 11 : C Traits : - 12 : C - 13 : C - 14 : C - 15 : C - 16 : C 17 : C Auteur : B. Richard (Dr - Ing.) - CEA/DEN/DANS/DM2S/SEMT/EMSI 18 : C====&===1=========2=========3=========4=========5=========6=========7== 19 : C 20 : C----DECLARATION GENERALES---------------------------------------------- 21 : C 22 : IMPLICIT INTEGER(I-N) 23 : IMPLICIT REAL*8(A-H,O-Z) 24 : C 25 : C----LES INCLUDES------------------------------------------------------- 26 : C 27 : -INC CCOPTIO 28 : -INC DECHE 29 : C 30 : C----DECLARATION PARAMETRES--------------------------------------------- 31 : C 32 : SEGMENT WRKK2 33 : REAL*8 EPSILI(NSTRS1) 34 : ENDSEGMENT 35 : 36 : SEGMENT WRK3 37 : REAL*8 EPSILO(NSTRS1) 38 : ENDSEGMENT 39 : 40 : INTEGER NSTRS1,NVARI,ISTRS,JSTRS,KCAS,IRTD 41 : 42 : REAL*8 EPS33(3,3),EPSIPP(3),VALP33(3,3) 43 : REAL*8 YOUN,XNU,XYDT,XDELTA,EPSIP1(3) 44 : REAL*8 XZERO,UN,DEUX,XPETIT 45 : 46 : PARAMETER (XZERO=0.D0 , UN=1.D0 , DEUX=2.D0, XPETIT=1.D-12) 47 : C 48 : C----MISE EN DONNEE----------------------------------------------------- 49 : C 50 : N=NSTRS1 51 : 52 : CMATE = 'ISOTROPE' 53 : 54 : C Module d Young et coefficient de Poisson 55 : YOUN = XMAT(1) 56 : XNU = XMAT(2) 57 : 58 : C Resistance et parametre endommagement traction 59 : XFT = XMAT(5) 60 : XEPSUT = XMAT(6) 61 : XY0T = XFT/YOUN 62 : 63 : C Resistance et parametre endommagement compression 64 : XFC = XMAT(7) 65 : XEPSUC = XMAT(8) 66 : XY0C = XFC/YOUN 67 : C 68 : C----VARIABLES INTERNES------------------------------------------------- 69 : C 70 : C Deformation totale 71 : EPSILI(1) = VAR0(1) 72 : EPSILI(2) = VAR0(2) 73 : EPSILI(3) = VAR0(3) 74 : EPSILI(4) = VAR0(4) 75 : 76 : C Endommagement et ecrouissage traction 77 : XDT = VAR0(7) 78 : XZT = VAR0(8) 79 : 80 : C Endommagement et ecrouissage traction 81 : XDC = VAR0(9) 82 : XZC = VAR0(10) 83 : 84 : C 85 : C----ACTUALISATION DEFORMATION------------------------------------------ 86 : C 87 : SEGINI WRK3 88 : DO 100 ISTRS=1,NSTRS1 89 : EPSILO(ISTRS)=EPSILI(ISTRS)+DEPST(ISTRS) 90 : 100 CONTINUE 91 : 92 : C 93 : C----PASSAGE EN MATRICIEL----------------------------------------------- 94 : C 94 : CALL ENDOCB (EPSILO,EPS33,2,IFOUR2) 95 : C 96 : C----CALCUL DE L ENDOMMAGEMENT------------------------------------------ 97 : C 98 : 99 : C Indicateur pour le sens du chargement 100 : XINDQ = SIG0(1)+SIG0(2) 101 : 102 : C Diagonalisation des deformations de membrane 103 : C ! La diagonalisation est analytique car absence de ! 104 : C ! convergence ... voire resultats faux ... ! 105 : C ! en cas de distorsion pure ! 106 : XDELTA = (EPS33(1,1)-EPS33(2,2))**2.0D0+EPS33(1,2)**2.0D0 107 : XDELTA = DSQRT(XDELTA) 108 : 109 : EPSIPP(1) = (EPS33(1,1) + EPS33(2,2) + XDELTA)/2.0D0 110 : EPSIPP(2) = (EPS33(1,1) + EPS33(2,2) - XDELTA)/2.0D0 111 : EPSIPP(3) = -(EPSIPP(1) + EPSIPP(2))*XNU / (UN-XNU) 112 : 113 : C Energie due a l endommagement traction 114 : DO I = 1,3 115 : EPSIP1(I) = 0.5D0*(DABS(EPSIPP(I))+EPSIPP(I)) 116 : ENDDO 117 : 118 : XYDT = 0.0D0 119 : DO I = 1,3 120 : XYDT = XYDT + (EPSIP1(I))**2.0D0 121 : ENDDO 122 : 123 : XYDT = DSQRT(XYDT) 124 : 125 : C Seuil d'endommagement traction 126 : IF (XINDQ.GE.0.0D0) THEN 127 : XINDQT = 1.0D0 128 : ELSE 129 : XINDQT = 0.0D0 130 : ENDIF 131 : 132 : XFDT = XYDT*XINDQT -(XY0T + XZT) 133 : 134 : C Calcul de l'endommagement en traction 135 : IF (XFDT.GT.0.0D0) THEN 136 : XDT = 1.0D0-(XY0T/XYDT)*DEXP(XEPSUT*(XY0T-XYDT)*XFT) 137 : XZT = XYDT - XY0T 138 : ENDIF 139 : 140 : IF (XDT.GE.0.999D0) THEN 141 : XDT = 0.999D0 142 : ENDIF 143 : 144 : C Seuil d'endommagement compression 145 : IF (XINDQ.LT.0.0D0) THEN 146 : XINDQC = 1.0D0 147 : ELSE 148 : XINDQC = 0.0D0 149 : ENDIF 150 : 151 : XFDC = XYDT*XINDQC -(XY0C + XZC) 152 : 153 : IF (XFDC.GT.0.0D0) THEN 154 : XDC = 1.0D0-(XY0C/XYDT)*DEXP(XEPSUC*(XY0C-XYDT)*XFC) 155 : XZC = XYDT -XY0C 156 : ENDIF 157 : C 158 : C----MATRICE DE HOOKE--------------------------------------------------- 159 : C 160 : CMATE = 'ISOTROPE' 161 : KCAS=1 162 : 162 : CALL DOHMAS(XMAT,CMATE,IFOUR2,NSTRS1,KCAS,DDHOOK,IRTD) 163 : 163 : CALL MATVE1 (DDHOOK,EPSILO,NSTRS1,NSTRS1,SIGF,2) 164 : C 165 : C----SORTIE DES CONTRAINTES--------------------------------------------- 166 : C 167 : DO 500 ISTRS=1,NSTRS1 168 : SIGF(ISTRS)=SIGF(ISTRS)*(1.0D0-XDT*XINDQT)*(1.0D0-XDC*XINDQC) 169 : 500 CONTINUE 170 : C 171 : C----SORTIE DES VARIABLES INTERNES-------------------------------------- 172 : C 173 : 174 : C Deformation totale 175 : VARF(1) = EPSILO(1) 176 : VARF(2) = EPSILO(2) 177 : VARF(3) = EPSILO(3) 178 : VARF(4) = EPSILO(4) 179 : 180 : C Endommagement et ecrouissage traction 181 : VARF(7) = XDT 182 : VARF(8) = XZT 183 : 184 : C Endommagement et ecrouissage traction 185 : VARF(9) = XDC 186 : VARF(10) = XZC 187 : 188 : SEGSUP WRK3 189 : RETURN 190 : END 191 : 192 : 193 : 194 : 195 : 196 : 197 : 198 :
© Cast3M 2003 - All rights reserved.
Disclaimer