OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sigeps70.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| sigeps70 ../starter/source/materials/mat/mat070/sigeps70.F
25!||--- called by ------------------------------------------------------
26!|| mulaw ../starter/source/materials/mat_share/mulaw.F
27!||--- calls -----------------------------------------------------
28!|| finter ../starter/source/tools/curve/finter.F
29!|| table_mat_vinterp ../starter/source/materials/tools/table_mat_vinterp.F
30!||--- uses -----------------------------------------------------
31!|| table_mat_vinterp_mod ../starter/source/materials/tools/table_mat_vinterp.F
32!||====================================================================
33 SUBROUTINE sigeps70(
34 1 NEL ,NUPARAM,NUVAR ,MFUNC ,KFUNC ,NPF ,
35 2 TF ,TIME ,TIMESTEP,UPARAM ,RHO0 ,RHO ,
36 3 VOLUME ,EINT ,NVARTMP ,VARTMP ,
37 4 DEPSXX ,DEPSYY ,DEPSZZ ,DEPSXY ,DEPSYZ ,DEPSZX ,
38 5 EPSXX ,EPSYY ,EPSZZ ,EPSXY ,EPSYZ ,EPSZX ,
39 6 SIG0XX ,SIG0YY ,SIG0ZZ ,SIG0XY ,SIG0YZ ,SIG0ZX ,
40 7 SIGNXX ,SIGNYY ,SIGNZZ ,SIGNXY ,SIGNYZ ,SIGNZX ,
41 8 SIGVXX ,SIGVYY ,SIGVZZ ,SIGVXY ,SIGVYZ ,SIGVZX ,
42 9 SOUNDSP,VISCMAX,UVAR ,OFF ,NGL ,
43 A PM ,IPM ,IMAT ,MATPARAM)
44C-----------------------------------------------
45C M o d u l e s
46C-----------------------------------------------
47 USE matparam_def_mod
48 USE table_mat_vinterp_mod
49C-----------------------------------------------
50C I m p l i c i t T y p e s
51C-----------------------------------------------
52#include "implicit_f.inc"
53C-----------------------------------------------
54C G l o b a l P a r a m e t e r s
55C-----------------------------------------------
56#include "mvsiz_p.inc"
57C---------+---------+---+---+--------------------------------------------
58C VAR | SIZE |TYP| RW| DEFINITION
59C---------+---------+---+---+--------------------------------------------
60C NEL | 1 | I | R | SIZE OF THE ELEMENT GROUP NEL
61C NUPARAM | 1 | I | R | SIZE OF THE USER PARAMETER ARRAY
62C NUVAR | 1 | I | R | NUMBER OF USER ELEMENT VARIABLES
63C---------+---------+---+---+--------------------------------------------
64C MFUNC | 1 | I | R | NUMBER FUNCTION USED FOR THIS USER LAW not used
65C KFUNC | NFUNC | I | R | FUNCTION INDEX not used
66C NPF | * | I | R | FUNCTION ARRAY
67C TF | * | F | R | FUNCTION ARRAY
68C---------+---------+---+---+--------------------------------------------
69C TIME | 1 | F | R | CURRENT TIME
70C TIMESTEP| 1 | F | R | CURRENT TIME STEP
71C UPARAM | NUPARAM | F | R | USER MATERIAL PARAMETER ARRAY
72C RHO0 | NEL | F | R | INITIAL DENSITY
73C RHO | NEL | F | R | DENSITY
74C VOLUME | NEL | F | R | VOLUME
75C EINT | NEL | F | R | TOTAL INTERNAL ENERGY
76C ... | | | |
77C DEPSXX | NEL | F | R | STRAIN INCREMENT XX
78C DEPSYY | NEL | F | R | STRAIN INCREMENT YY
79C ... | | | |
80C EPSXX | NEL | F | R | STRAIN XX
81C EPSYY | NEL | F | R | STRAIN YY
82C ... | | | |
83C SIG0XX | NEL | F | R | OLD ELASTO PLASTIC STRESS XX
84C SIG0YY | NEL | F | R | OLD ELASTO PLASTIC STRESS YY
85C ... | | | |
86C---------+---------+---+---+--------------------------------------------
87C SIGNXX | NEL | F | W | NEW ELASTO PLASTIC STRESS XX
88C SIGNYY | NEL | F | W | NEW ELASTO PLASTIC STRESS YY
89C ... | | | |
90C SIGVXX | NEL | F | W | VISCOUS STRESS XX
91C SIGVYY | NEL | F | W | VISCOUS STRESS YY
92C ... | | | |
93C SOUNDSP | NEL | F | W | SOUND SPEED (NEEDED FOR TIME STEP)
94C VISCMAX | NEL | F | W | MAXIMUM DAMPING MODULUS(NEEDED FOR TIME STEP)
95C---------+---------+---+---+--------------------------------------------
96C UVAR |NEL*NUVAR| F |R/W| USER ELEMENT VARIABLE ARRAY
97C OFF | NEL | F |R/W| DELETED ELEMENT FLAG (=1. ON, =0. OFF)
98C---------+---------+---+---+--------------------------------------------
99#include "param_c.inc"
100C
101 INTEGER ,INTENT(IN) :: NVARTMP,IMAT
102 INTEGER ,INTENT(INOUT) :: VARTMP(NEL,NVARTMP)
103 INTEGER NEL, NUPARAM, NUVAR,IPT,
104 . NGL(NEL),MAT(NEL),IPLA,IPM(NPROPMI,*)
105 my_real
106 . TIME,TIMESTEP,UPARAM(*),
107 . RHO(NEL),RHO0(NEL),VOLUME(NEL),EINT(NEL),
108 . DEPSXX(NEL),DEPSYY(NEL),DEPSZZ(NEL),
109 . DEPSXY(NEL),DEPSYZ(NEL),DEPSZX(NEL),
110 . EPSXX(NEL) ,EPSYY(NEL) ,EPSZZ(NEL) ,
111 . epsxy(nel) ,epsyz(nel) ,epszx(nel) ,
112 . sig0xx(nel),sig0yy(nel),sig0zz(nel),
113 . sig0xy(nel),sig0yz(nel),sig0zx(nel),
114 . pm(npropm,*)
115 TYPE(matparam_struct_) ,INTENT(INOUT) :: MATPARAM
116C-----------------------------------------------
117C O U T P U T A r g u m e n t s
118C-----------------------------------------------
119 my_real
120 . SIGNXX(NEL),SIGNYY(NEL),SIGNZZ(NEL),
121 . SIGNXY(NEL),SIGNYZ(NEL),SIGNZX(NEL),
122 . SIGVXX(NEL),SIGVYY(NEL),SIGVZZ(NEL),
123 . SIGVXY(NEL),SIGVYZ(NEL),SIGVZX(NEL),
124 . SOUNDSP(NEL),VISCMAX(NEL)
125C-----------------------------------------------
126C I N P U T O U T P U T A r g u m e n t s
127C-----------------------------------------------
128 my_real
129 . uvar(nel,nuvar), off(nel), pla(nel)
130C-----------------------------------------------
131C VARIABLES FOR FUNCTION INTERPOLATION
132C-----------------------------------------------
133 INTEGER NPF(*), MFUNC, KFUNC(MFUNC)
134 my_real
135 . TF(*),FINTER
136 EXTERNAL FINTER
137C EXTERNAL FINTER
138C Y = FINTER(IFUNC(J),X,NPF,TF,DYDX)
139C Y : y = f(x)
140C X : x
141C DYDX : f'(x) = dy/dx
142C IFUNC(J): FUNCTION INDEX
143C J : FIRST(J=1), SECOND(J=2) .. FUNCTION USED FOR THIS LAW
144C NPF,TF : FUNCTION PARAMETER
145C-----------------------------------------------
146C L o c a l V a r i a b l e s
147C-----------------------------------------------
148 INTEGER I,J,J1,J2,I1,I2,II,IADBUF,NC,NFUNC,ILOAD(MVSIZ),
149 . IFUNC(MVSIZ,100),NLOAD,NUNLOAD,
150 . IE_CST(MVSIZ),ILOAD0,IDAMAGE,IFLAG,IPFUNC,NUMTABL
151 INTEGER :: NINDEX_EPSSM,NINDEX_EPST
152 INTEGER, DIMENSION(MVSIZ) :: INDEX_EPSSM,INDEX_EPST
153 my_real R,FAC,E,DELTA,ALPHA,SVM2,P,YFACC,
154 . SVM,INTER_0,EMAX,EPSSMAX,EPS_MAX,DSIG,EXPO,HYS,DE,DF,MU
155 my_real YLDELAS(MVSIZ),YLDMIN(MVSIZ),YLDMAX(MVSIZ),EPSS(MVSIZ),
156 . E0(MVSIZ),G(MVSIZ),NU(MVSIZ),AA1(MVSIZ),AA2(MVSIZ),
157 . EPST(MVSIZ),AA(MVSIZ),RATE(MVSIZ,2),
158 . yfac(mvsiz,2),sig0(mvsiz),eps0(mvsiz)
159 INTEGER, DIMENSION(NEL) :: IPOS1
160 my_real, DIMENSION(NEL) :: XVEC1
161 INTEGER, DIMENSION(NEL,2) :: IPOS
162 my_real, DIMENSION(NEL,2) :: xvec
163 my_real, DIMENSION(MVSIZ) :: yld,dfc
164C-----------------------------------------------
165C USER VARIABLES INITIALIZATION
166C-----------------------------------------------
167C
168 nfunc = ipm(10,imat)
169 DO i= 1, nel
170 DO j=1,nfunc
171 ifunc(i,j)=ipm(10+j,imat)
172 ENDDO
173 ENDDO
174C
175 iadbuf = ipm(7,imat)-1
176 nload = uparam(iadbuf+7)
177 nunload = uparam(iadbuf+8)
178 inter_0 = uparam(2*nfunc + 15)
179 epssmax = uparam(iadbuf+4)
180 DO i=1,nel
181 e0(i) = uparam(iadbuf+2)
182 aa(i) = uparam(iadbuf+3)
183 g(i) = uparam(iadbuf+5)
184 nu(i) = uparam(iadbuf+6)
185 idamage = uparam(iadbuf+ 2*nfunc + 9)
186 expo = uparam(iadbuf+ 2*nfunc + 10)
187 hys = uparam(iadbuf+ 2*nfunc + 11)
188 emax = uparam(iadbuf+ 2*nfunc + 12)
189 iflag = uparam(iadbuf+ 2*nfunc + 13)
190 yfacc = uparam(iadbuf+ 8 + 2*nfunc )
191 ENDDO
192 numtabl = matparam%NTABLE
193C-----------------------------------------------
194 DO i=1,nel
195C-------------------
196C epst_spherique
197C-------------------
198 epst(i) = epsxx(i)**2+epsyy(i)**2 + epszz(i)**2 +
199 . half*(epsxy(i)**2+epsyz(i)**2+epszx(i)**2)
200 epst(i) = sqrt(epst(i))
201 epst(i) = min(one, epst(i))
202C-------------------
203C estimation elastique
204C-------------------
205 eps0(i) = uvar(i,1)
206 sig0(i) = uvar(i,2)
207 ENDDO
208C-------------------
209C CRITERE
210C-------------------
211 nindex_epst = 0
212 nindex_epssm = 0
213 DO i=1,nel
214c YLD_elast
215 IF (epst(i) >= epssmax) THEN
216 yldelas(i)= inter_0
217 yldelas(i)= yldelas(i) + emax*(epst(i) - epssmax)
218 nindex_epssm = nindex_epssm + 1
219 index_epssm(nindex_epssm) = i
220 ELSE
221 nindex_epst = nindex_epst + 1
222 index_epst(nindex_epst) = i
223 ENDIF
224 ENDDO
225c---------------------------------------------
226 IF (nindex_epst > 0) THEN
227 IF (matparam%TABLE(1)%NDIM == 1) THEN
228 DO ii=1,nindex_epst
229 i = index_epst(ii)
230 ipos1(i) = vartmp(i,1)
231 xvec1(i) = epst(i)
232 ENDDO
233c
234 CALL table_mat_vinterp(matparam%TABLE(1),nel,nindex_epst,ipos1,xvec1,yld,dfc)
235c
236 DO ii=1,nindex_epst
237 i = index_epst(ii)
238 yldelas(i) = yld(ii)
239 vartmp(i,1) = ipos1(ii)
240 ENDDO
241c
242 ELSE
243 DO ii=1,nindex_epst
244 i = index_epst(ii)
245 ipos(ii,1) = vartmp(i,1)
246 ipos(ii,2) = 1
247 xvec(ii,1) = epst(i)
248 xvec(ii,2) = uparam(iadbuf + 9) ! static strain rate
249 ENDDO
250c
251 CALL table_mat_vinterp(matparam%TABLE(1),nel,nindex_epst,ipos,xvec,yld,dfc)
252c
253 DO ii=1,nindex_epst
254 i = index_epst(ii)
255 yldelas(i) = yld(ii)
256 vartmp(i,1) = ipos(ii,1)
257 ENDDO
258 ENDIF
259c
260 ENDIF
261c---------------------------------------------
262c yld_max
263c---------------------------------------------
264c
265 IF (matparam%TABLE(1)%NDIM == 1) THEN
266 DO i=1,nel
267 ipos1(i) = vartmp(i,1)
268 xvec1(i) = epst(i)
269 ENDDO
270c
271 CALL table_mat_vinterp(matparam%TABLE(1),nel,nel,ipos1,xvec1,yldmax,dfc)
272c
273 ELSE
274 DO i=1,nel
275 ipos(i,1) = vartmp(i,1)
276 ipos(i,2) = 1
277 xvec(i,1) = epst(i)
278 xvec(i,2) = uparam(iadbuf + 9) ! static strain rate
279 ENDDO
280c
281 CALL table_mat_vinterp(matparam%TABLE(1),nel,nel,ipos,xvec,yldmax,dfc)
282c
283 END IF
284c
285c yld_min
286c
287 IF (nunload > 0) THEN
288 IF (matparam%TABLE(2)%NDIM == 1) THEN
289 DO i=1,nel
290 ipos1(i) = 1
291 xvec1(i) = epst(i)
292 ENDDO
293c
294 CALL table_mat_vinterp(matparam%TABLE(2),nel,nel,ipos1,xvec1,yldmin,dfc)
295c
296 ELSE
297 DO i=1,nel
298 ipos(i,1) = 1
299 ipos(i,2) = 1
300 xvec(i,1) = epst(i)
301 xvec(i,2) = uparam(iadbuf + 9 + nload) ! static strain rate
302 ENDDO
303c
304 CALL table_mat_vinterp(matparam%TABLE(2),nel,nel,ipos,xvec,yldmin,dfc)
305c
306 END IF
307 ELSE
308 yldmin(1:nel) = yldmax(1:nel)
309 END IF
310c-----------------------------------------------------
311 IF(iflag > 0) THEN
312 DO i=1,nel
313 alpha = one/max(em20,uvar(i,10))
314 sig0xx(i) = alpha*sig0xx(i)
315 sig0yy(i) = alpha*sig0yy(i)
316 sig0zz(i) = alpha*sig0zz(i)
317 sig0xy(i) = alpha*sig0xy(i)
318 sig0zx(i) = alpha*sig0zx(i)
319 sig0yz(i) = alpha*sig0yz(i)
320 ENDDO
321 ENDIF
322C
323C =====================================================
324 DO i = 1,nel
325 iload0 = uvar(i,5)
326 ie_cst(i)= 0
327 delta = epst(i) - uvar(i,4)
328 e = uvar(i,3)
329 uvar(i,9) = uvar(i,4)
330 IF(delta >= zero)THEN
331 yld(i) = yldmax(i)
332 iload(i) = 1
333 ELSEIF(delta < zero)THEN
334 yld(i) = yldmin(i)
335 iload(i) = -1
336 IF(idamage /= 0 )THEN
337 yld(i) = yldmax(i)
338 ENDIF
339 ENDIF
340C
341 epss(i) = epst(i) - yld(i)/ e
342 epss(i) = max(zero, epss(i))
343 de = aa(i)*(epss(i) - uvar(i,1))
344c
345 IF(iload(i) == 1) THEN
346 e = e + max(de, zero)
347 IF(iload0 == -1) e= uvar(i,3)
348 uvar(i,1) = max(uvar(i,1), epss(i))
349 ELSE
350 e = e + min(de ,zero)
351 IF(iload0 == 1) e= uvar(i,3)
352 uvar(i,1) = min(epss(i),uvar(i,1))
353 ENDIF
354 e = min(e, emax)
355 e = max(e, e0(i))
356 uvar(i,3) = e
357C==================================================
358 aa1(i) = e*(one-nu(i))/(one + nu(i))/(one - two*nu(i))
359 aa2(i) = aa1(i)*nu(i)/(one - nu(i))
360 g(i) =half*e/(one + nu(i))
361C ----
362 signxx(i)= aa1(i)*depsxx(i) + aa2(i)*(depsyy(i) + depszz(i))
363 signyy(i)= aa1(i)*depsyy(i) + aa2(i)*(depsxx(i) + depszz(i))
364 signzz(i)= aa1(i)*depszz(i) + aa2(i)*(depsxx(i) + depsyy(i))
365 signxy(i)= g(i) *depsxy(i)
366 signyz(i)= g(i) *depsyz(i)
367 signzx(i)= g(i) *depszx(i)
368 dsig = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
369 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
370 dsig =sqrt(dsig)
371C estimate stress
372 signxx(i)=sig0xx(i) + aa1(i)*depsxx(i)
373 . + aa2(i)*(depsyy(i) + depszz(i))
374 signyy(i)=sig0yy(i) + aa1(i)*depsyy(i)
375 . + aa2(i)*(depsxx(i) + depszz(i))
376 signzz(i)=sig0zz(i) + aa1(i)*depszz(i)
377 . + aa2(i)*(depsxx(i) + depsyy(i))
378 signxy(i)=sig0xy(i) + g(i) *depsxy(i)
379 signyz(i)=sig0yz(i) + g(i) *depsyz(i)
380 signzx(i)=sig0zx(i) + g(i) *depszx(i)
381
382 svm2 = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
383 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
384 svm =sqrt(svm2)
385C sound velocity
386 soundsp(i) = sqrt(aa1(i)/rho0(i))
387 viscmax(i) = zero
388C
389 IF(idamage == 0 ) THEN
390 IF(delta >= zero ) THEN
391 yld(i) = min(yldmax(i),svm)
392 ELSEIF(delta < zero) THEN
393 yld(i) = max(yldmin(i),svm)
394 yld(i) = min(yld(i),uvar(i,6))
395 ie_cst(i) = 1
396 IF( dsig > yldmin(i) .AND. dsig > svm)THEN
397 yld(i) = yldmin(i)
398 ie_cst(i) = 0
399 ENDIF
400 ENDIF
401 ELSE
402 yld(i) = yldmax(i)
403 IF(delta > zero .AND. svm < yldmax(i))yld(i)=svm
404 uvar(i,8) = uvar(i,8) +
405 . half*(yld(i) + uvar(i,6))*delta
406 uvar(i,8) = max(zero, uvar(i,8))
407 uvar(i,2) = max(uvar(i,2) , uvar(i,8))
408 ENDIF
409
410 ENDDO
411C
412C-------------------
413C projection spherique
414C-------------------
415
416 IF(idamage == 0 ) THEN
417 DO i=1,nel
418 signxx(i)= aa1(i)*epsxx(i) + aa2(i)*(epsyy(i) + epszz(i))
419 signyy(i)= aa1(i)*epsyy(i) + aa2(i)*(epsxx(i) + epszz(i))
420 signzz(i)= aa1(i)*epszz(i) + aa2(i)*(epsxx(i) + epsyy(i))
421 signxy(i)= g(i) *epsxy(i)
422 signyz(i)= g(i) *epsyz(i)
423 signzx(i)= g(i) *epszx(i)
424C
425 svm2 = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
426 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
427 svm =sqrt(svm2)
428 r = yld(i)/max(em20,svm)
429 signxx(i)=signxx(i)*r
430 signyy(i)=signyy(i)*r
431 signzz(i)=signzz(i)*r
432 signxy(i)=signxy(i)*r
433 signyz(i)=signyz(i)*r
434 signzx(i)=signzx(i)*r
435C
436 IF(ie_cst(i) == 1) THEN
437 iload0 = int(uvar(i,5))
438 IF(iload0 /= iload(i))THEN
439 iload(i) = iload0
440 uvar(i,1) = eps0(i)
441 ENDIF
442 ENDIF
443 uvar(i,4) = epst(i)
444 uvar(i,2) = svm*r
445 uvar(i,5) = iload(i)
446 uvar(i,6) = yld(i)
447 ENDDO
448 ELSEIF(idamage == 1) THEN
449 DO i=1,nel
450C
451 r = one
452 signxx(i)= aa1(i)*epsxx(i) + aa2(i)*(epsyy(i) + epszz(i))
453 signyy(i)= aa1(i)*epsyy(i) + aa2(i)*(epsxx(i) + epszz(i))
454 signzz(i)= aa1(i)*epszz(i) + aa2(i)*(epsxx(i) + epsyy(i))
455 signxy(i)= g(i) *epsxy(i)
456 signyz(i)= g(i) *epsyz(i)
457 signzx(i)= g(i) *epszx(i)
458C
459 svm2 = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
460 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
461 svm =sqrt(svm2)
462 r = (yld(i)/max(em20,svm))
463 signxx(i)=signxx(i)*r
464 signyy(i)=signyy(i)*r
465 signzz(i)=signzz(i)*r
466 signxy(i)=signxy(i)*r
467 signyz(i)=signyz(i)*r
468 signzx(i)=signzx(i)*r
469C
470 IF(iload(i) == -1) THEN
471 r = (yldmin(i)/max(em20,yldelas(i)))
472 p = third*(signxx(i) + signyy(i) + signzz(i))
473 signxx(i)=(signxx(i) - p)*r + p
474 signyy(i)=(signyy(i) - p)*r + p
475 signzz(i)=(signzz(i) - p)*r + p
476 signxy(i)=signxy(i)*r
477 signyz(i)=signyz(i)*r
478 signzx(i)=signzx(i)*r
479 ENDIF
480C
481 svm2 = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
482 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
483 svm =sqrt(svm2)
484 uvar(i,4) = epst(i)
485 uvar(i,5) = iload(i)
486 uvar(i,6) = yld(i)
487 ENDDO
488 ELSEIF(idamage == 2 ) THEN
489 DO i=1,nel
490C
491 signxx(i)= aa1(i)*epsxx(i) + aa2(i)*(epsyy(i) + epszz(i))
492 signyy(i)= aa1(i)*epsyy(i) + aa2(i)*(epsxx(i) + epszz(i))
493 signzz(i)= aa1(i)*epszz(i) + aa2(i)*(epsxx(i) + epsyy(i))
494 signxy(i)= g(i) *epsxy(i)
495 signyz(i)= g(i) *epsyz(i)
496 signzx(i)= g(i) *epszx(i)
497C
498 svm2 = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
499 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
500 svm =sqrt(svm2)
501 r = yld(i)/max(em20,svm)
502 signxx(i)=signxx(i)*r
503 signyy(i)=signyy(i)*r
504 signzz(i)=signzz(i)*r
505 signxy(i)=signxy(i)*r
506 signyz(i)=signyz(i)*r
507 signzx(i)=signzx(i)*r
508C
509 IF(iload(i) == -1) THEN
510 r = yldmin(i)/max(em20,yldelas(i))
511 signxx(i)=signxx(i)*r
512 signyy(i)=signyy(i)*r
513 signzz(i)=signzz(i)*r
514 signxy(i)=signxy(i)*r
515 signyz(i)=signyz(i)*r
516 signzx(i)=signzx(i)*r
517 ENDIF
518C
519 svm2 = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
520 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
521 svm =sqrt(svm2)
522 uvar(i,4) = epst(i)
523 uvar(i,5) = iload(i)
524 uvar(i,6) = yld(i)
525 ENDDO
526
527 ELSEIF(idamage == 3) THEN
528 DO i=1,nel
529C
530 signxx(i)= aa1(i)*epsxx(i) + aa2(i)*(epsyy(i) + epszz(i))
531 signyy(i)= aa1(i)*epsyy(i) + aa2(i)*(epsxx(i) + epszz(i))
532 signzz(i)= aa1(i)*epszz(i) + aa2(i)*(epsxx(i) + epsyy(i))
533 signxy(i)= g(i) *epsxy(i)
534 signyz(i)= g(i) *epsyz(i)
535 signzx(i)= g(i) *epszx(i)
536C
537 svm2 = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
538 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
539 svm =sqrt(svm2)
540 r = (yld(i)/max(em20,svm))
541 signxx(i)=signxx(i)*r
542 signyy(i)=signyy(i)*r
543 signzz(i)=signzz(i)*r
544 signxy(i)=signxy(i)*r
545 signyz(i)=signyz(i)*r
546 signzx(i)=signzx(i)*r
547C
548 IF(iload(i) == -1) THEN
549 r = one - (uvar(i,8)/max(em20,uvar(i,2)))**expo
550 r = one - (one - hys)*r
551 p = third*(signxx(i) + signyy(i) + signzz(i))
552 signxx(i)=(signxx(i) - p)*r + p
553 signyy(i)=(signyy(i) - p)*r + p
554 signzz(i)=(signzz(i) - p)*r + p
555 signxy(i)=signxy(i)*r
556 signyz(i)=signyz(i)*r
557 signzx(i)=signzx(i)*r
558 ENDIF
559C
560 uvar(i,4) = epst(i)
561 uvar(i,5) = iload(i)
562 uvar(i,6) = yld(i)
563 ENDDO
564 ELSEIF(idamage == 4) THEN
565 DO i=1,nel
566C
567 signxx(i)= aa1(i)*epsxx(i) + aa2(i)*(epsyy(i) + epszz(i))
568 signyy(i)= aa1(i)*epsyy(i) + aa2(i)*(epsxx(i) + epszz(i))
569 signzz(i)= aa1(i)*epszz(i) + aa2(i)*(epsxx(i) + epsyy(i))
570 signxy(i)= g(i) *epsxy(i)
571 signyz(i)= g(i) *epsyz(i)
572 signzx(i)= g(i) *epszx(i)
573C
574 svm2 = (signxx(i)**2 + signyy(i)**2 + signzz(i)**2)
575 . + two*(signxy(i)**2 + signzx(i)**2 + signyz(i)**2)
576 svm =sqrt(svm2)
577 r = (yld(i)/max(em20,svm))
578 signxx(i)=signxx(i)*r
579 signyy(i)=signyy(i)*r
580 signzz(i)=signzz(i)*r
581 signxy(i)=signxy(i)*r
582 signyz(i)=signyz(i)*r
583 signzx(i)=signzx(i)*r
584C
585 IF(iload(i) == -1) THEN
586 r = one - (uvar(i,8)/max(em20,uvar(i,2)))**expo
587 r = one - (one - hys)*r
588 signxx(i)=signxx(i)*r
589 signyy(i)=signyy(i)*r
590 signzz(i)=signzz(i)*r
591 signxy(i)=signxy(i)*r
592 signyz(i)=signyz(i)*r
593 signzx(i)=signzx(i)*r
594 ENDIF
595 uvar(i,4) = epst(i)
596 uvar(i,5) = iload(i)
597 uvar(i,6) = yld(i)
598 ENDDO
599 ENDIF
600C
601 IF(iflag > 0 ) THEN
602 ipfunc = ifunc(1,nfunc)
603 DO i=1,nel
604 mu = 1 - rho(i)/rho0(i)
605 alpha = one
606 IF(mu > 0) THEN
607 alpha = yfacc*finter(ipfunc,mu,npf,tf,df)
608 alpha = max(zero, alpha)
609 ENDIF
610 signxx(i) = alpha*signxx(i)
611 signyy(i) = alpha*signyy(i)
612 signzz(i) = alpha*signzz(i)
613 signxy(i) = alpha*signxy(i)
614 signzx(i) = alpha*signzx(i)
615 signyz(i) = alpha*signyz(i)
616 uvar(i,10) = alpha
617 ENDDO
618 ENDIF
619C------------------------------------
620 RETURN
621 END
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
table_mat_vinterp_mod
Definition table_mat_vinterp.F:32
table_mat_vinterp_mod::table_mat_vinterp
subroutine table_mat_vinterp(table, dimx, nel, ipos, xx, yy, dydx)
Definition table_mat_vinterp.F:45
sigeps70
subroutine sigeps70(nel, nuparam, nuvar, mfunc, kfunc, npf, tf, time, timestep, uparam, rho0, rho, volume, eint, nvartmp, vartmp, depsxx, depsyy, depszz, depsxy, depsyz, depszx, epsxx, epsyy, epszz, epsxy, epsyz, epszx, sig0xx, sig0yy, sig0zz, sig0xy, sig0yz, sig0zx, signxx, signyy, signzz, signxy, signyz, signzx, sigvxx, sigvyy, sigvzz, sigvxy, sigvyz, sigvzx, soundsp, viscmax, uvar, off, ngl, pm, ipm, imat, matparam)
Definition sigeps70.F:44