OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fail_johnson_ib.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!|| fail_johnson_ib ../engine/source/materials/fail/johnson_cook/fail_johnson_ib.F
25!||--- called by ------------------------------------------------------
26!|| fail_beam18 ../engine/source/elements/beam/fail_beam18.F
27!||====================================================================
28 SUBROUTINE fail_johnson_ib(
29 . NEL ,NGL ,IPT ,NUPARAM ,UPARAM ,
30 . TIME ,TSTAR ,SIGNXX ,SIGNXY ,SIGNZX ,
31 . DPLA ,EPSD ,OFF ,FOFF ,DFMAX ,
32 . TDEL ,IOUT ,ISTDO ,UELR , NPG)
33C-----------------------------------------------
34c Johnson-Cook failure model for integrated beams
35C-----------------------------------------------
36C I m p l i c i t T y p e s
37C-----------------------------------------------
38#include "implicit_f.inc"
39C-----------------------------------------------
40C G l o b a l P a r a m e t e r s
41C-----------------------------------------------
42#include "comlock.inc"
43C-----------------------------------------------
44C D u m m y A r g u m e n t s
45C-----------------------------------------------
46 INTEGER ,INTENT(IN) :: NEL ! size of element group
47 INTEGER ,INTENT(IN) :: NUPARAM ! size of parameter array
48 INTEGER ,INTENT(IN) :: IPT ! current integration point
49 INTEGER ,INTENT(IN) :: IOUT ! output file unit
50 INTEGER ,INTENT(IN) :: ISTDO ! output file unit
51 INTEGER ,DIMENSION(NEL) ,INTENT(IN) :: NGL ! table of element identifiers
52 my_real ,INTENT(IN) :: time ! current time
53 my_real ,DIMENSION(NUPARAM) ,INTENT(IN) :: uparam ! failure model parameter array
54 my_real ,DIMENSION(NEL) ,INTENT(IN) :: dpla ! plastic strain increment
55 my_real ,DIMENSION(NEL) ,INTENT(IN) :: epsd ! strain rate
56 my_real ,DIMENSION(NEL) ,INTENT(IN) :: tstar ! normalized temperature
57 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: off ! element deactivation flag
58 my_real ,DIMENSION(NEL) ,INTENT(IN) :: signxx ! stress component
59 my_real ,DIMENSION(NEL) ,INTENT(IN) :: signxy ! stress component
60 my_real ,DIMENSION(NEL) ,INTENT(IN) :: signzx ! stress component
61 INTEGER ,DIMENSION(NEL) ,INTENT(INOUT) :: FOFF ! integration point deactivation flag
62 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: dfmax ! maximum damage
63 my_real ,DIMENSION(NEL) ,INTENT(INOUT) :: tdel ! desactivation time
64 my_real, DIMENSION(NEL) ,INTENT(INOUT) :: uelr
65 INTEGER ,INTENT(IN) :: NPG
66C-----------------------------------------------
67C L o c a l V a r i a b l e s
68C-----------------------------------------------
69 INTEGER :: I,J,NINDX,FAILIP
70 INTEGER ,DIMENSION(NEL) :: INDX
71 my_real :: D1,D2,D3,D4,D5,EPSD0,P,EPSF,SVM,EPSF_MIN
72C=======================================================================
73 nindx = 0
74 d1 = uparam(1)
75 d2 = uparam(2)
76 d3 = uparam(3)
77 d4 = uparam(4)
78 d5 = uparam(5)
79 epsd0 = uparam(6)
80 epsf_min = uparam(12)
81 failip = min(nint(uparam(13)),npg)
82c-----------------------------
83 DO i=1,nel
84 IF (off(i) == one .and. foff(i) == 1 .and. dpla(i) > zero) THEN
85 p = third*signxx(i)
86 svm = signxx(i)**2 + three*(signxy(i)**2 + signzx(i)**2)
87 svm = sqrt(svm)
88 epsf = d3 * p / max(em20, svm)
89 epsf = d1 + d2 * exp(epsf)
90 IF (d4 /= zero) epsf = epsf * (one + d4*log(max(one,epsd(i)/epsd0)))
91 IF (d5 /= zero) epsf = epsf * (one + d5*tstar(i))
92 epsf = max(epsf,epsf_min)
93 IF (epsf > zero) dfmax(i) = dfmax(i) + dpla(i) / epsf
94 IF (dfmax(i) >= one) THEN ! total damage in integration point
95 nindx = nindx + 1
96 indx(nindx) = i
97 foff(i) = 0
98 tdel(i) = time
99 dfmax(i) = one
100 uelr(i) = uelr(i) + one
101 IF (nint(uelr(i)) >= failip) THEN
102 off(i) = four_over_5
103 ENDIF
104 ENDIF
105 ENDIF
106 ENDDO
107c------------------------
108 IF (nindx > 0) THEN
109 DO j=1,nindx
110 i = indx(j)
111#include "lockon.inc"
112 WRITE(iout, 1000) ngl(i),ipt,time
113 WRITE(istdo,1000) ngl(i),ipt,time
114#include "lockoff.inc"
115 IF (off(i) == four_over_5) THEN
116 WRITE(iout, 1111) ngl(i),time
117 WRITE(istdo,1111) ngl(i),time
118 ENDIF
119 END DO
120 END IF ! NINDX
121c------------------
122 1000 FORMAT(5x,'FAILURE (JC) OF BEAM ELEMENT ',i10,1x,',INTEGRATION PT',i5
123 . ,2x,'AT TIME :',1pe12.4)
124
125 1111 FORMAT(1x,'DELETED BEAM ELEMENT ',i10,1x,'at time :',1PE12.4)
126c------------------
127 RETURN
128 END
my_real
#define my_real
Definition cppsort.cpp:32
fail_johnson_ib
subroutine fail_johnson_ib(nel, ngl, ipt, nuparam, uparam, time, tstar, signxx, signxy, signzx, dpla, epsd, off, foff, dfmax, tdel, iout, istdo, uelr, npg)
Definition fail_johnson_ib.F:33
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21