1 : C INCRE4 SOURCE CHAT 05/01/13 00:35:20 5004 2 : SUBROUTINE INCRE4(SIG,VAR,EPSVPT,VARPT,XMAT,ALFA,NSTRS0, 3 : & NVARI,NCOMAT) 4 : IMPLICIT INTEGER(I-N) 5 : IMPLICIT REAL*8(A-H,O-Z) 6 : DIMENSION SIG(*),VAR(*),EPSVPT(*),VARPT(*),XMAT(*) 7 : DIMENSION SM(4),SB(4),XXM1(4),XXB1(4),XXM2(4),XXB2(4) 8 : DIMENSION GAMM(4),GABB(4),EPMM(4),EPBB(4),XXM(4),XXB(4) 9 : DIMENSION ANSM(4),ANSB(4),ANM(4),ANB(4) 10 : DIMENSION EPS0(8) 11 : PARAMETER (AMAX = 1.0E20) 12 : DETIER=2.0D0/3.0D0 13 : ROOT=SQRT(1.5D0) 14 : C-------------------------------------------------------------------| 15 : C STEP-1 SEPARATING MEMBRANE & BENDING PARTS | 16 : C-------------------------------------------------------------------I 17 : NSTRS = NSTRS0+2 18 : NS = NSTRS/2 19 : SM(1)=SIG(1) 20 : SM(2)=SIG(2) 21 : SM(3)=0.0D0 22 : IF (NS.EQ.4) SM(4)=SIG(3) 23 : SB(1)=SIG(NS) 24 : SB(2)=SIG(NS+1) 25 : SB(3)=0.0D0 26 : IF (NS.EQ.4) SB(4)=SIG(NS+2) 27 : DO 20 I=1,NS 28 : XXM1(I) = VAR( I) 29 : XXM2(I) = VAR(2*NS+I) 30 : XXB1(I) = VAR( NS+I) 31 : XXB2(I) = VAR(3*NS+I) 32 : GAMM(I) = VAR(4*NS+I) 33 : GABB(I) = VAR(5*NS+I) 34 : EPMM(I) = VAR(6*NS+I) 35 : EPBB(I) = VAR(7*NS+I) 36 : 20 CONTINUE 37 : C-------------------------------------------------------------------| 38 : C STEP-2: EVALUATION OF J2 SIGMA - X | 39 : C-------------------------------------------------------------------I 40 : XXM(1)= (2.D0*SM(1) -SM(2))/3.D0 -XXM1(1) -XXM2(1) 41 : XXM(2)= (2.D0*SM(2) -SM(1))/3.D0 -XXM1(2) -XXM2(2) 42 : XXM(3)= -1.D0*(SM(1)+SM(2))/3.D0 -XXM1(3) -XXM2(3) 43 : XXM(4)=0.0D0 44 : IF(NS.EQ.4) XXM(4)= SM(4) -XXM1(4) -XXM2(4) 45 : XXB(1)= (2.D0*SB(1) -SB(2))/3.D0 -XXB1(1) -XXB2(1) 46 : XXB(2)= (2.D0*SB(2) -SB(1))/3.D0 -XXB1(2) -XXB2(2) 47 : XXB(3)= -1.D0*(SB(1)+SB(2))/3.D0 -XXB1(3) -XXB2(3) 48 : XXB(4)=0.0D0 49 : IF(NS.EQ.4) XXB(4)= SB(4) -XXB1(4) -XXB2(4) 50 : * T1 = XXM(1)*XXM(1) + XXM(2)*XXM(2) + XXM(1)*XXM(2) 51 : * IF(NS.EQ.4) T1 = T1 + XXM(4)*XXM(4) 52 : T1 = PROCON(XXM,XXM,NS) 53 : AJM2 = 3.0D0*T1/2.D0 54 : * T1 = XXB(1)*XXB(1) + XXB(2)*XXB(2) + XXB(1)*XXB(2) 55 : * IF(NS.EQ.4) T1 = T1 + XXB(4)*XXB(4) 56 : T1 = PROCON(XXB,XXB,NS) 57 : AJB2 = 3.0D0*T1/2.D0 58 : AJ2 = SQRT( AJM2 + ALFA*ALFA*AJB2 ) 59 : C--------------------------------------------------------------------| 60 : C STEP-3: CALCUL OF SIGV/K | 61 : C -------------------------------------------------------------------I 62 : RR = VAR (4*NSTRS+2) 63 : RS = XMAT(10)*RR 64 : SK = (AJ2 - RS-XMAT(7))/( XMAT(8) + XMAT(9)*RR ) 65 : C--------------------------------------------------------------------| 66 : C STEP-4: CALCULATION OF EFFECTIVE INELASTIC STRAIN INCREMENT (P) | 67 : C -------------------------------------------------------------------I 68 : IF (SK.GT.0.0D0) THEN 69 : PPT1= XPUISS(SK,XMAT(6) ,AMAX) 70 : PPT2= XPUISS(SK,XMAT(6)+1,AMAX) 71 : ELSE 72 : PPT1= 0.0D0 73 : PPT2= 0.0D0 74 : ENDIF 75 : T1 = XMAT(5)*PPT2 76 : PPT2=1.0D20 77 : IF(DABS(T1).LT.40.0D0) PPT2= DEXP(T1) 78 : PPT = PPT1 * PPT2 79 : VARPT(4*NSTRS+1) = PPT 80 : C---------------------------------------------------------------------| 81 : C******* CALCULATION OF INELASTIC STRAIN INCREMENTS (EPSV0) | 82 : C---------------------------------------------------------------------I 83 : IF(AJ2.EQ.0.0D0) AJ2=1.0D-20 84 : T1 = 1.5D0*PPT/AJ2 85 : EPS0(1) = T1*XXM(1) 86 : EPS0(2) = T1*XXM(2) 87 : EPS0(3) = T1*XXM(3) 88 : IF(NS.EQ.4) EPS0(4)= T1*XXM(4) 89 : EPS0(NS+1) = T1*XXB(1) *ALFA *ALFA 90 : EPS0(NS+2) = T1*XXB(2) *ALFA *ALFA 91 : EPS0(NS+3) = T1*XXB(3) *ALFA *ALFA 92 : IF(NS.EQ.4) EPS0(NS+4)= T1*XXB(4) *ALFA*ALFA 93 : DO 30 I=1,NSTRS 94 : 30 VARPT(3*NSTRS+I)=EPS0(I) 95 : C 96 : C----------------------------------------------------------------| 97 : C ******* CALCULATION OF PI(P) | 98 : C----------------------------------------------------------------I 99 : P = VAR(4*NSTRS+1) 100 : T1 = XMAT(16)*P 101 : T2 = DEXP(-T1) 102 : PIP= XMAT(13)+(1.0D0-XMAT(13)) * T2 103 : C----------------------------------------------------------------| 104 : C******* CALCULATION OF XII | 105 : C----------------------------------------------------------------I 106 : C CORRECTIONS MLR 26/2/93 COMME DANS INCRE2 107 : C 108 : * T1 = 3.0D0*(XXM1(1)*XXM1(1)+XXM1(2)*XXM1(2)+XXM1(1)*XXM1(2)) 109 : * IF(NS.EQ.4) T1 = T1 + 3.0D0 * XXM1(4) * XXM1(4) 110 : T1=PROCON(XXM1,XXM1,NS) 111 : T1=3.D0*T1/2.D0 112 : * T2 = 3.0D0*(XXB1(1)*XXB1(1)+XXB1(2)*XXB1(2)+XXB1(1)*XXB1(2)) 113 : * IF(NS.EQ.4) T2 = T2 + 3.0D0 * XXB1(4) * XXB1(4) 114 : T2=PROCON(XXB1,XXB1,NS) 115 : T2=3.D0*T2/2.D0 116 : X1II = 0.D0 117 : IF(XMAT(17).NE.0.D0) 118 : . X1II = SQRT(T1+ALFA*ALFA*T2)/XMAT(17) 119 : * T1 = 3.0D0*(XXM2(1)*XXM2(1)+XXM2(2)*XXM2(2)+XXM2(1)*XXM2(2)) 120 : * IF(NS.EQ.4) T1 = T1 + 3.0D0 * XXM2(4) * XXM2(4) 121 : T1=PROCON(XXM2,XXM2,NS) 122 : T1=3.D0*T1/2.D0 123 : * T2 = 3.0D0*(XXB2(1)*XXB2(1)+XXB2(2)*XXB2(2)+XXB2(1)*XXB2(2)) 124 : * IF(NS.EQ.4) T2 = T2 + 3.0D0 * XXB2(4) * XXB2(4) 125 : T2=PROCON(XXB2,XXB2,NS) 126 : T2=3.D0*T2/2.D0 127 : X2II = 0.D0 128 : IF(XMAT(18).NE.0.D0) 129 : . X2II = SQRT(T1+ALFA*ALFA*T2)/XMAT(18) 130 : COX1 = XPUISS(X1II,XMAT(19)-1.0D0,AMAX) 131 : COX2 = XPUISS(X2II,XMAT(20)-1.0D0,AMAX) 132 : C ----------------------------------------------------------------| 133 : C******* CALCULATION OF DX1 AND DX2 | 134 : C-----------------------------------------------------------------I 135 : DO 73 I=1,NSTRS 136 : VARPT( I)=0.0D0 137 : 73 VARPT(NSTRS+I)=0.0D0 138 : DO 74 I=1,NSTRS 139 : T1 = DETIER * XMAT(11) * XMAT(12) * EPS0(I) 140 : T2 = XMAT(11) * PIP * VAR(I) * PPT 141 : T3 = COX1 * VAR(I) 142 : VARPT(I)= T1 - T2 -T3 143 : T1 = DETIER * XMAT(14) * XMAT(15) * EPS0(I) 144 : T2 = XMAT(14) * PIP * VAR(NSTRS+I) * PPT 145 : T3 = COX2 * VAR(NSTRS+I) 146 : 74 VARPT(NSTRS+I)= T1 - T2 - T3 147 : C-----------------------------------------------------------------| 148 : C******* CALCULATION OF DR | 149 : C-----------------------------------------------------------------I 150 : T1 = 1.0D0-(1.0D0-VAR(4*NSTRS+3)/XMAT(25))**2 151 : QR = VAR(4*NSTRS+3)-XMAT(26)*T1 152 : CO1= QR-VAR(4*NSTRS+2) 153 : CO0= ABS(CO1) 154 : CO = XPUISS(CO0,XMAT(22)-1.0D0,AMAX) 155 : T1= XMAT(16)*(VAR(4*NSTRS+3)-VAR(4*NSTRS+2))*PPT 156 : T2= XMAT(21)*CO*CO1 157 : DR= T1 + T2 158 : VARPT(4*NSTRS+2)=DR 159 : C----------------------------------------------------------------| 160 : C******* CALCULATION OF PROD (N X N*) | 161 : C----------------------------------------------------------------I 162 : TM1 = EPMM(1)-GAMM(1) 163 : TM2 = EPMM(2)-GAMM(2) 164 : TM3 = EPMM(3)-GAMM(3) 165 : TM4 = 0.0D0 166 : IF(NS.EQ.4) TM4 = EPMM(4)-GAMM(4) 167 : AIM2 = TM1*TM1+TM2*TM2+TM3*TM3+2.D0*TM4*TM4 168 : AIM2 = AIM2*3.D0/2.D0 169 : TB1 = EPBB(1)-GABB(1) 170 : TB2 = EPBB(2)-GABB(2) 171 : TB3 = EPBB(3)-GABB(3) 172 : TB4 = 0.0D0 173 : IF(NS.EQ.4) TB4 = EPBB(4)-GABB(4) 174 : AIB2 = TB1*TB1+TB2*TB2+TB3*TB3+2.D0*TB4*TB4 175 : AIB2 = AIB2*3.D0/2.D0 176 : AI2 = SQRT(AIM2+(AIB2/ALFA/ALFA)) 177 : IF (AI2.LT.1.D-10) THEN 178 : IF (AI2.EQ.0.D0) AI2=1.D-10 179 : AI20=AI2*1.D20 180 : TM1=TM1*1.D20 181 : TM2=TM2*1.D20 182 : TM3=TM3*1.D20 183 : TM4=TM4*1.D20 184 : TB1=TB1*1.D20 185 : TB2=TB2*1.D20 186 : TB3=TB3*1.D20 187 : TB4=TB4*1.D20 188 : ELSE 189 : AI20=AI2 190 : ENDIF 191 : ANSM(1) = ROOT*TM1/AI20 192 : ANSM(2) = ROOT*TM2/AI20 193 : ANSM(3) = ROOT*TM3/AI20 194 : ANSM(4) = 0.0D0 195 : IF(NS.EQ.4) ANSM(4) = ROOT*TM4/AI20 196 : ANSB(1) = ROOT*TB1/AI20 197 : ANSB(2) = ROOT*TB2/AI20 198 : ANSB(3) = ROOT*TB3/AI20 199 : ANSB(4) = 0.0D0 200 : IF(NS.EQ.4) ANSB(4) = ROOT*TB4/AI20 201 : ANM (1) = ROOT*XXM(1)/AJ2 202 : ANM (2) = ROOT*XXM(2)/AJ2 203 : ANM (3) = ROOT*XXM(3)/AJ2 204 : ANM (4) = 0.0D0 205 : IF(NS.EQ.4) ANM (4) = ROOT*XXM(4)/AJ2 206 : ANB (1) = ROOT*XXB(1)/AJ2 207 : ANB (2) = ROOT*XXB(2)/AJ2 208 : ANB (3) = ROOT*XXB(3)/AJ2 209 : ANB (4) = 0.0D0 210 : IF(NS.EQ.4) ANB (4) = ROOT*XXB(4)/AJ2 211 : 76 CONTINUE 212 : FF = DETIER*AI2-VAR(4*NSTRS+4) 213 : HF=0.0D0 214 : IF(FF.GT.0.0D0) HF=1.0D0 215 : ANAN1=PROCON(ANM,ANSM,NS) 216 : * ANAN1 = ANM(1)*ANSM(1)+ANM(2)*ANSM(2) 217 : * ANAN1 = ANAN1+2.D0*ANM(3)*ANSM(3) 218 : ANANSM=0.0D0 219 : IF(ANAN1.GT.0.0D0) ANANSM=ANAN1 220 : ANAN1=PROCON(ANB,ANSB,NS) 221 : * ANAN1 = ANB(1)*ANSB(1)+ANB(2)*ANSB(2) 222 : * ANAN1 = ANAN1+2.D0*ANB(3)*ANSB(3) 223 : ANANSB=0.0D0 224 : IF(ANAN1.GT.0.0D0) ANANSB=ANAN1 225 : C----------------------------------------------------------------| 226 : C******* CALCULATION OF DQ | 227 : C----------------------------------------------------------------I 228 : ANANS=ANANSM+ANANSB/ALFA/ALFA 229 : DQ = XMAT(23)*HF*ANANS*PPT 230 : VARPT(4*NSTRS+4)=DQ 231 : C----------------------------------------------------------------| 232 : C******* CALCULATION OF DGETA | 233 : C----------------------------------------------------------------I 234 : T1 = ROOT*(1.0D0-XMAT(23))*HF*ANANS*PPT 235 : DO 78 I=1,NS 236 : VARPT(2*NSTRS+I) = T1*ANSM(I) 237 : VARPT(2*NSTRS+NS+I)= T1*ANSB(I) 238 : 78 CONTINUE 239 : C----------------------------------------------------------------| 240 : C******* CALCULATION OF DQQ | 241 : C----------------------------------------------------------------I 242 : VARPT(4*NSTRS+3)=2.0D0*XMAT(24)*(XMAT(25)-VAR(4*NSTRS+3))*DQ 243 : C 244 : EPSVPT(1)=EPS0(1) 245 : EPSVPT(2)=EPS0(2) 246 : IF (NS.EQ.4) EPSVPT(3)=EPS0(4) 247 : EPSVPT(NS)=EPS0(NS+1) 248 : EPSVPT(NS+1)=EPS0(NS+2) 249 : IF (NS.EQ.4) EPSVPT(6)=EPS0(8) 250 : C 251 : RETURN 252 : END 253 : 254 : 255 : 256 :
© Cast3M 2003 - All rights reserved.
Disclaimer