1 : C ARED3M SOURCE CHAT 05/01/12 21:23:28 5004 2 : SUBROUTINE ARED3M(SIGT,EPST,ENDT,XLIMT, 3 : * SIGD,EPSD,ENDD,XLIMD, 4 : * DEPS,VDEP,VDEP1,NPA,VAL,VP,EPSPLS,VAL1,VP1, 5 : * VAL2,VAL3,VP2,VP3,EPSTT,FDDD,D,SIFF,VDEII, 6 : * VDEP2,EPSC,PPLUS,U,PROD,ENDOPL,VDEI,KM,AAA) 7 : C====================================================================== 8 : C TEST D ENDOMMAGEMENT INITIAL POUR LES ELEMENTS TRIDIMENSIONNELS 9 : C====================================================================== 10 : C 11 : C 12 : C MODIF. : F.CORMERY 13 : C E.N.S.M.A - LMPM 14 : C JUIN 1993 15 : C 16 : C====================================================================== 17 : C ENTREES : SIGT,EPST,ENDT, VALEURS AU PAS PRECEDENT 18 : C DSIG,DEPS ACCROISSEMENTS INITIAUX 19 : C SORTIES SIGD,EPSD,ENDD, VALEURS A LA FIN DU PAS 20 : C DEPS ACCROISSEMENTS 21 : C VDEP MATRICE TANGENTE 22 : C====================================================================== 23 : IMPLICIT INTEGER(I-N) 24 : IMPLICIT REAL*8(A-H,O-Z) 25 : REAL*8 LAMB,MUT 26 : C---------------------------------------------------------------------- 27 : C COMMON 28 : C---------------------------------------------------------------------- 29 : COMMON /APENTE/ LAMB,MUT,G,ALP1,BETA 30 : C---------------------------------------------------------------------- 31 : C DATA ET DIMENSIONS 32 : C---------------------------------------------------------------------- 33 : C N9 34 : CC DIMENSION SIGD(6),EPSD(6),ENDD(6),D(3,3),VAL1(3,3),VP1(3), 35 : C N15 N51 36 : CC * SIGT(6),EPST(6),ENDT(6),SIFF(6),VDEII(6,6), 37 : C 38 : CC * DEPS(6),VDEP(6,6),VDEP1(6,6),VAL(3,3),VP(3), 39 : C N87 40 : CC * VAL2(3,3),VAL3(3,3),VP3(3),VP2(3),VDEP2(6,6) 41 : C N93 N174 N255 42 : CC DIMENSION EPSC(6),EPSPLS(6),PPLUS(3,3,3,3),U(3,3,3,3) 43 : C N261 N267 N303 N309 44 : CC DIMENSION PROD(6),ENDOPL(6),VDEI(6,6),KM(6),EPSTT(6) 45 : C 46 : C 47 : DIMENSION AAA(*) 48 : DIMENSION SIGD(*),EPSD(*),ENDD(*),D(3,*),VAL1(3,*),VP1(*), 49 : C N15 N51 50 : * SIGT(*),EPST(*),ENDT(*),SIFF(*),VDEII(6,*), 51 : C 52 : * DEPS(*),VDEP(6,*),VDEP1(6,*),VAL(3,*),VP(*), 53 : C N87 54 : * VAL2(3,*),VAL3(3,*),VP3(*),VP2(*),VDEP2(6,*) 55 : C N93 N174 N255 56 : DIMENSION EPSC(*),EPSPLS(*),PPLUS(3,3,3,*),U(3,3,3,*) 57 : C N261 N267 N303 N309 58 : DIMENSION PROD(*),ENDOPL(*),VDEI(6,*),KM(*),EPSTT(*) 59 : DIMENSION NII(6),MII(6),NIII(6),MIII(6) 60 : DATA NII/1,2,3,3,3,2/ 61 : DATA MII/1,2,3,2,1,1/ 62 : DATA NIII/1,2,3,2,3,3/ 63 : DATA MIII/1,2,3,1,2,1/ 64 : DATA ZERO/0.D0/ 65 : DATA IMA/6/,DEUX/2.D0/,UD/0.5D0/ 66 : C---------------------------------------------------------------------- 67 : MT=6 68 : IKK=1 69 : C---------------------------------------------------------------------- 70 : C TENSEUR DE DEFORMATION ACTUALISEE 71 : C---------------------------------------------------------------------- 72 : MMI=1 73 : GOTO 121 74 : C------------ 75 : 3 CONTINUE 76 : C------------ 77 : DO 144 I=1,6 78 : M= MII(I) 79 : N=NII(I) 80 : DO 144 J=1,6 81 : K=MII(J) 82 : L=NII(J) 83 : VDEI(I,J)=U(M,N,K,L) 84 : VDEII(I,J)=U(M,N,K,L) 85 : 144 CONTINUE 86 : C---------------------------Calcul de la deformation initiale---------- 87 : IF(NPA.EQ.1)THEN 88 : C------------ 88 : CALL INVERE(VDEI,6,6,KM,DET) 89 : C------------ 90 : DO 4 I=1,6 91 : EPST(I)=0.D0 92 : DO 4 J=1,6 93 : EPST(I)=EPST(I)+VDEI(I,J)*(-G*ENDT(J)) 94 : 4 CONTINUE 95 : C------------ 96 : MMI=0 97 : ENDIF 98 : C---------------------------------------------------------------------- 99 : 234 CONTINUE 100 : DO 10 I=1,IMA 101 : EPSD(I)=EPST(I)+DEPS(I) 102 : EPSC(I)=EPSD(I) 103 : 10 CONTINUE 104 : C---------------------------------------------------------------------- 105 : C CALCUL DU TENSEUR DE DEFORMATION POSITIVE 106 : C---------------------------------------------------------------------- 106 : 7 CALL EPPLUM(EPSC,PPLUS,1,VAL1,VP1,AAA(442), 107 : . AAA(478),AAA(514),AAA(517),AAA(526), 108 : . AAA(532)) 109 : C---------------- 110 : EPSD(4)=EPSD(4)/2 111 : EPSD(5)=EPSD(5)/2 112 : EPSD(6)=EPSD(6)/2 113 : C---------------- 113 : CALL T4CT2(PPLUS,EPSD,EPSPLS,6) 114 : C---------------- 115 : DO 324 I=1,6 116 : IF(ABS(EPSPLS(I)).LE.1E-15)EPSPLS(I)=0.D0 117 : 324 CONTINUE 118 : C---------------- 119 : EPSD(4)=EPSD(4)*2 120 : EPSD(5)=EPSD(5)*2 121 : EPSD(6)=EPSD(6)*2 122 : EPSPLS(4)=EPSPLS(4)*2 123 : EPSPLS(5)=EPSPLS(5)*2 124 : EPSPLS(6)=EPSPLS(6)*2 125 : C---------------------------------------------------------------------- 126 : C TEST D ENDOMMAGEMENT 127 : 128 : C---------------------------------------------------------------------- 128 : CALL AVENDO(EPSPLS,ENDT,DEQ) 129 : CRIT=DEQ-XLIMT 130 : IF(CRIT.LE.ZERO)GO TO 100 131 : C---------------------------------------------------------------------- 132 : C CALCUL DU MULTIPLICATEUR 133 : C---------------------------------------------------------------------- 134 : XLMD=ZERO 134 : CALL AMULTM(EPSPLS,ENDT,XLMD,XLIMT) 135 : C---------------------------------------------------------------------- 136 : C ACTUALISATION DES GRANDEURS 137 : C---------------------------------------------------------------------- 137 : CALL AACTUM(SIGT,EPST,ENDT,XLIMT,SIGD,EPSD,ENDD,XLIMD, 138 : * VDEP,EPSC,EPSPLS,XLMD,PPLUS, 139 : * VDEP1,SIFF,VDEP2,FDDD,AAA(442),AAA(448), 140 : * AAA(454),AAA(460),AAA(466),AAA(472), 141 : * AAA(478),AAA(484),AAA(493),AAA(502), 142 : * AAA(511),AAA(520),AAA(529),AAA(538), 143 : * AAA(547),AAA(628),AAA(709),AAA(790), 144 : * AAA(871)) 145 : ISS=2 146 : GO TO 899 147 : C********************************************************************** 148 : C CAS ELASTIQUE 149 : C********************************************************************** 150 : 100 CONTINUE 151 : ISS=1 152 : C---------------------------------------------------------------------- 153 : DO 110 I=1,IMA 154 : ENDD(I)=ENDT(I) 155 : 110 CONTINUE 156 : XLIMD=XLIMT 157 : C--------------------------------CALCUL PPLUS:D=D+ ------------ 158 : * MODEE1=1 159 : * IF(MODEE1.EQ.1)THEN 160 : DO 123 I=1,6 161 : 123 ENDOPL(I)=ENDD(I) 162 : * GOTO 108 163 : * ENDIF 163 : * CALL T4CT2(PPLUS,ENDD,ENDOPL,6) 164 : C--------------------------------CALCUL DES CONTRAINTES---------------- 165 : *108 TRAC2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3) 166 : TRAC2=EPSD(1)*ENDD(1)+EPSD(2)*ENDD(2)+EPSD(3)*ENDD(3) 167 : * +EPSD(4)*ENDD(4)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6) 168 : TRAC3=EPSD(1)+EPSD(2)+EPSD(3) 169 : PROD(1)=DEUX*EPSD(1)*ENDD(1)+EPSD(5)*ENDD(5)+EPSD(6)*ENDD(6) 170 : PROD(2)=DEUX*EPSD(2)*ENDD(2)+EPSD(4)*ENDD(4)+EPSD(6)*ENDD(6) 171 : PROD(3)=DEUX*EPSD(3)*ENDD(3)+EPSD(5)*ENDD(5)+EPSD(4)*ENDD(4) 172 : PROD(4)=UD*EPSD(6)*ENDD(5)+EPSD(2)*ENDD(4)+UD*EPSD(4)*ENDD(3)+ 173 : * UD*EPSD(5)*ENDD(6)+UD*EPSD(4)*ENDD(2)+EPSD(3)*ENDD(4) 174 : PROD(5)=UD*EPSD(5)*ENDD(1)+UD*EPSD(4)*ENDD(6)+EPSD(3)*ENDD(5)+ 175 : * EPSD(1)*ENDD(5)+UD*EPSD(6)*ENDD(4)+UD*EPSD(5)*ENDD(3) 176 : PROD(6)=UD*EPSD(6)*ENDD(1)+EPSD(2)*ENDD(6)+UD*EPSD(4)*ENDD(5)+ 177 : * EPSD(1)*ENDD(6)+UD*EPSD(6)*ENDD(2)+UD*EPSD(5)*ENDD(4) 178 : C--------------------- 179 : DO 115 I=1,3 180 : SIGD(I)=G*ENDOPL(I)+LAMB*TRAC3+ALP1*(TRAC2+TRAC3*ENDD(I)) 181 : * +DEUX*BETA*PROD(I)+DEUX*MUT*EPSD(I) 182 : 115 CONTINUE 183 : C--------------------- 184 : DO 116 I=4,6 185 : SIGD(I)=G*ENDOPL(I)+ALP1*TRAC3*ENDD(I) 186 : * +DEUX*BETA*PROD(I)+MUT*EPSD(I) 187 : 116 CONTINUE 188 : C------------------ 189 : DO 17 I=1,6 190 : SIFF(I)=0 191 : DO 17 J=1,6 192 : SIFF(I)=SIFF(I)+VDEII(I,J)*EPSC(J) 193 : 17 CONTINUE 194 : DO 18 I=1,6 195 : SIFF(I)=SIFF(I)+G*ENDOPL(I) 196 : 18 CONTINUE 197 : C---------------------------------------------------------------------- 198 : MMI=0 199 : C---------------------------------------------------------------------- 200 : C CALCUL MATRICE TANGENTE ELASTIQUE ANISOTROPE 201 : C---------------------------------------------------------------------- 202 : C 203 : C---------------------------------Passage en notation indicielle------- 204 : 121 CONTINUE 205 : DO 42 I=1,3 206 : D(I,I)=ENDT(I) 207 : DO 43 J=I+1,3 208 : K = 9 - I -J 209 : * IF ((I.EQ.1).AND.(J.EQ.2)) K=6 210 : * IF ((I.EQ.1).AND.(J.EQ.3)) K=5 211 : * IF ((I.EQ.2).AND.(J.EQ.3)) K=4 212 : D(I,J)=ENDT(K) 213 : D(J,I)=ENDT(K) 214 : 43 CONTINUE 215 : 42 CONTINUE 216 : 217 : C----------------------------Calcul de la matrice tangente------------- 218 : DO 51 I=1,3 219 : DO 51 J=1,3 220 : DO 51 K=1,3 221 : DO 51 L=1,3 222 : C------------------------- 223 : U(I,J,K,L)=LAMB*DELT(I,J)*DELT(K,L)+MUT*(DELT(I,K)*DELT(J,L) 224 : * +DELT(I,L)*DELT(J,K))+ALP1*(DELT(I,J)*D(K,L)+D(I,J)*DELT(K,L)) 225 : * +BETA*(D(I,K)*DELT(J,L)+D(L,J)*DELT(I,K)+DELT(I,L) 226 : * *D(J,K)+D(I,L)*DELT(J,K)) 227 : C------------------------- 228 : 51 CONTINUE 229 : C-------------------------------Notation de VOIGT-------------------- 230 : IF(MMI.EQ.1)GOTO 3 231 : DO 134 I=1,6 232 : M=MIII(I) 233 : N=NIII(I) 234 : DO 134 J=1,6 235 : K= MIII(J) 236 : L= NIII(J) 237 : VDEP(I,J)=U(M,N,K,L) 238 : 134 CONTINUE 239 : C------------------------------------------------------------------- 240 : DO 142 I=1,6 241 : DO 142 J=1,6 242 : VDEP1(I,J)=VDEP(I,J) 243 : VDEP2(I,J)=VDEP(I,J) 244 : 142 CONTINUE 245 : C******************************************************************** 246 : C 247 : C******************************************************************** 248 : 899 DO 321 I=1,6 249 : IF(ABS(SIGD(I)).LE.1E-15)SIGD(I)=0.D0 250 : 321 CONTINUE 251 : C-------------------------------------------------------------------- 251 : CALL VECT(SIGD,VAL,VP,AAA(442),AAA(451),AAA(457)) 252 : DO 332 I=4,6 253 : 332 EPSPLS(I)=EPSPLS(I)/2 253 : CALL VECT(EPSPLS,VAL2,VP2,AAA(442),AAA(451),AAA(457)) 254 : DO 333 I=4,6 255 : 333 EPSPLS(I)=EPSPLS(I)*2 255 : CALL VECT(ENDD,VAL3,VP3,AAA(442),AAA(451),AAA(457)) 256 : C-------------------------------------------------------------------- 256 : CALL INVERE(VDEP2,6,6,KM,DET) 257 : C------------ 258 : DO 401 I=1,6 259 : EPSTT(I)=0.D0 260 : DO 401 J=1,6 261 : EPSTT(I)=EPSTT(I)+VDEP2(I,J)*(-G*ENDD(J)) 262 : 401 CONTINUE 263 : RETURN 264 : END 265 : 266 : 267 :
© Cast3M 2003 - All rights reserved.
Disclaimer