OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_fail_johnson.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!|| hm_read_fail_johnson ../starter/source/materials/fail/johnson_cook/hm_read_fail_johnson.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_fail ../starter/source/materials/fail/hm_read_fail.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
30!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
31!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
32!||--- uses -----------------------------------------------------
33!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
34!|| message_mod ../starter/share/message_module/message_mod.F
35!|| submodel_mod ../starter/share/modules1/submodel_mod.F
36!||====================================================================
37 SUBROUTINE hm_read_fail_johnson(
38 . FAIL ,MAT_ID ,FAIL_ID ,ILAW ,IRUPT ,IXFEM ,
39 . TITR ,LSUBMODEL,UNITAB )
40C-----------------------------------------------
41c ROUTINE DESCRIPTION :
42c Read Johnoson-Cook failure model parameters
43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE fail_param_mod
47 USE unitab_mod
48 USE message_mod
49 USE submodel_mod
50 USE hm_option_read_mod
51 USE names_and_titles_mod , ONLY : nchartitle
52C-----------------------------------------------
53C I m p l i c i t T y p e s
54C-----------------------------------------------
55#include "implicit_f.inc"
56C-----------------------------------------------
57#include "units_c.inc"
58C-----------------------------------------------
59C D u m m y A r g u m e n t s
60C-----------------------------------------------
61 INTEGER ,INTENT(IN) :: FAIL_ID ! failure model ID
62 INTEGER ,INTENT(IN) :: MAT_ID ! material law ID
63 INTEGER ,INTENT(IN) :: ILAW ! material law type number
64 INTEGER ,INTENT(IN) :: IRUPT ! failure model type number
65 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR ! material model title
66 TYPE(unit_type_) ,INTENT(IN) :: UNITAB ! table of input units
67 TYPE(submodel_data),INTENT(IN) :: LSUBMODEL(*) ! submodel table
68 INTEGER ,INTENT(INOUT) :: IXFEM ! XFEM activation flag
69 TYPE(fail_param_) ,INTENT(INOUT) :: FAIL ! failure model data structure
70C-----------------------------------------------
71C L o c a l V a r i a b l e s
72C-----------------------------------------------
73 INTEGER IFAIL_SH,ISOLID,I,FAILIP
74 my_real D1, D2, D3, D4, D5, EPSP0, DADV, PTHKF, EPSF_MIN
75 my_real triax,epsf,depsf_dtriax
76 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
77C=======================================================================
78 is_encrypted = .false.
79 is_available = .false.
80C--------------------------------------------------
81C Check encryption
82C--------------------------------------------------
83 CALL hm_option_is_encrypted(is_encrypted)
84C--------------------------------------------------
85C EXTRACT PARAMETERS
86C--------------------------------------------------
87 CALL hm_get_floatv ('D1' ,d1 ,is_available,lsubmodel,unitab)
88 CALL hm_get_floatv ('D2' ,d2 ,is_available,lsubmodel,unitab)
89 CALL hm_get_floatv ('D3' ,d3 ,is_available,lsubmodel,unitab)
90 CALL hm_get_floatv ('D4' ,d4 ,is_available,lsubmodel,unitab)
91 CALL hm_get_floatv ('D5' ,d5 ,is_available,lsubmodel,unitab)
92 CALL hm_get_floatv ('Epsilon_Dot_0',epsp0 ,is_available,lsubmodel,unitab)
93 CALL hm_get_intv ('Ifail_sh' ,ifail_sh,is_available,lsubmodel)
94 CALL hm_get_intv ('Ifail_so' ,isolid ,is_available,lsubmodel)
95 CALL hm_get_floatv ('EPSF_MIN' ,epsf_min,is_available,lsubmodel,unitab)
96 CALL hm_get_floatv ('Dadv' ,dadv ,is_available,lsubmodel,unitab)
97 CALL hm_get_intv ('Ixfem' ,ixfem ,is_available,lsubmodel)
98 CALL hm_get_intv ('FAILIP' ,failip ,is_available,lsubmodel)
99 IF (failip == 0) failip = 1
100!
101 ! CHECK COMPATIBILITY WITH MATERIAL LAW
102 IF (ilaw == 25) THEN
103 CALL ancmsg(msgid=3001, msgtype=msgwarning, anmode=aninfo_blind,
104 . i1=mat_id,
105 . c1=titr)
106 ENDIF
107!
108 ! Check the Johnson-Cook Failure criterion
109 triax = one
110 epsf = d1 + d2*exp(d3*triax)
111 ! -> Check if criterion reaches negative values
112 IF ((epsf <= zero).AND.(epsf_min == zero)) THEN
113 DO i = 1,10
114 depsf_dtriax = d2*d3*exp(d3*triax)
115 triax = triax - epsf/depsf_dtriax
116 triax = min(triax, one)
117 triax = max(triax,-one)
118 epsf = d1 + d2*exp(d3*triax)
119 ENDDO
120 ! -> Print warning message
121 IF ((abs(epsf) < em06) .AND. (triax < one) .AND. (triax > -one)) THEN
122 CALL ancmsg(msgid=3058,
123 . msgtype=msgwarning,
124 . anmode=aninfo_blind,
125 . i1=mat_id,
126 . c1=titr,
127 . r1=triax)
128 ENDIF
129 ENDIF
130!
131 IF (epsp0 == zero) epsp0=em20
132 IF (ifail_sh == 0) ifail_sh=1
133 IF (isolid == 0) isolid=1
134 IF (ixfem /= 1 .AND. ixfem /= 2) ixfem = 0
135 IF (ixfem > 0) isolid = 0
136 IF (dadv == zero .OR. dadv > one) dadv = one
137c-----------------------------
138 IF (ifail_sh == 1) THEN
139 pthkf = em06
140 ELSEIF (ifail_sh == 2) THEN
141 pthkf = one
142 ENDIF
143c----------------------------------
144 fail%KEYWORD = 'JOHNSON-COOK'
145 fail%IRUPT = irupt
146 fail%FAIL_ID = fail_id
147 fail%NUPARAM = 13
148 fail%NIPARAM = 0
149 fail%NUVAR = 0
150 fail%NFUNC = 0
151 fail%NTABLE = 0
152 fail%NMOD = 0
153 fail%PTHK = pthkf
154c
155 ALLOCATE (fail%UPARAM(fail%NUPARAM))
156 ALLOCATE (fail%IPARAM(fail%NIPARAM))
157 ALLOCATE (fail%IFUNC (fail%NFUNC))
158 ALLOCATE (fail%TABLE (fail%NTABLE))
159c
160 fail%UPARAM(1) = d1
161 fail%UPARAM(2) = d2
162 fail%UPARAM(3) = d3
163 fail%UPARAM(4) = d4
164 fail%UPARAM(5) = d5
165 fail%UPARAM(6) = epsp0
166 fail%UPARAM(7) = ifail_sh
167 fail%UPARAM(8) = isolid
168 fail%UPARAM(9) = 0 ! not used
169 fail%UPARAM(10)= dadv
170 fail%UPARAM(11)= ixfem
171 fail%UPARAM(12)= epsf_min
172 fail%UPARAM(13)= failip
173c---------------------------
174c OUTPUT
175c---------------------------
176 IF (is_encrypted) THEN
177 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
178 ELSE
179
180 WRITE(iout,500) mat_id,irupt,fail_id
181
182C---
183C for shell
184C---
185 IF (ixfem == 0)THEN
186 WRITE(iout, 1000)d1,d2,d3,d4,d5,epsf_min,failip
187 IF (ifail_sh == 1)THEN
188 WRITE(iout, 1100)
189 ELSEIF(ifail_sh == 2)THEN
190 WRITE(iout, 1200)
191 ENDIF
192 ELSE
193 WRITE(iout, 1010)d1,d2,d3,d4,d5,ixfem,dadv,epsf_min
194 WRITE(iout, 1400)
195 END IF
196C---
197C for solid
198C---
199 IF (ixfem == 0) THEN ! XFEM is not available for solids
200 IF(isolid == 1) THEN
201 WRITE(iout, 2100)
202 ELSEIF (isolid == 2) THEN
203 WRITE(iout, 2200)
204 ELSEIF (isolid == 3) THEN
205 WRITE(iout, 2300)
206 ELSEIF (isolid == 4) THEN
207 WRITE(iout, 2400)
208 ENDIF
209 ENDIF
210C
211 IF (d3 > zero)THEN
212 CALL ancmsg(msgid=831, msgtype=msgwarning, anmode=aninfo,
213 . i1=mat_id)
214 ENDIF
215 ENDIF ! IF (IS_ENCRYPTED)
216c-----------------------------------------------------------
217 500 FORMAT(//
218 & 5x,'MAT_ID . . . . . . . . . . . . . . .=',i10/
219 & 5x,'failure model. . . . . . . . . . . .=',I10/
220 & 5X,'fail_id. . . . . . . . . . . . . . .=',I10/)
221 1000 FORMAT(
222 & 5X,' ---------------------------------------------------- ',/
223 & 5X,' failure criterion : johnson-cook ',/,
224 & 5X,' ---------------------------------------------------- ',/
225 & 5X,'first failure parameter(d1). . . . . . . . . . . . . .=',1PG20.13/
226 & 5X,'second failure parameter(d2). . . . . . . . . . . . . .=',1PG20.13/
227 & 5X,'third failure parameter(d3). . . . . . . . . . . . . .=',1PG20.13/
228 & 5X,'forth failure parameter(d4). . . . . . . . . . . . . .=',1PG20.13/
229 & 5X,'fifth failure parameter(d5). . . . . . . . . . . . . .=',1PG20.13/
230 & 5X,'lower plastic strain at failure(epsf_min). . . . . . . =',1PG20.13/
231 & 5X,'number of failed intg. points prior to elem deletion .=',I10/)
232 1100 FORMAT(
233 & 5X,' shell element deletion after failure')
234 2100 FORMAT(
235 & 5X,' solid element deletion after failure')
236
237 1200 FORMAT(
238 & 5X,' stress tensor in shell layer set to zero after failure')
239 1400 FORMAT(
240 & 5X,' shell element cracking after failure')
241 2200 FORMAT(
242 & 5X,' deviatoric stress tensor in solid will vanish after failure')
243 2300 FORMAT(
244 & 5X,' deviatoric stress tensor in solid will vanish after failure',/,
245 & 5X,' pressure can still be positive or null')
246 2400 FORMAT(
247 & 5X,' the stress tensor in solid will vanish after failure')
248C
249 1010 FORMAT(
250 & 5X,' ---------------------------------------------------- ',/
251 & 5X,' failure criterion : x-fem johnson-cook ',/,
252 & 5X,' ---------------------------------------------------- ',/
253 & 5X,'first failure parameter(d1). . . . . . . . . . . . . .=',1PG20.13/
254 & 5X,'second failure parameter(d2). . . . . . . . . . . . . .=',1PG20.13/
255 & 5X,'third failure parameter(d3). . . . . . . . . . . . . .=',1PG20.13/
256 & 5X,'forth failure parameter(d4). . . . . . . . . . . . . .=',1PG20.13/
257 & 5X,'fifth failure parameter(d5). . . . . . . . . . . . . .=',1PG20.13/
258 & 5X,'flag xfem. . . . . . . . . . . . . . . . . . . . . . . .=',I10/
259 & 5X,'critical advancement VALUE . . . . . . . . . . . . . . .=',1PG20.13/
260 & 5X,'lower plastic strain at failure(epsf_min). . . . . . . =',1PG20.13//)
261C-----------
262 RETURN
263 END
the
end diagonal values have been computed in the(sparse) matrix id.SOL
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288
hm_read_fail_johnson
subroutine hm_read_fail_johnson(fail, mat_id, fail_id, ilaw, irupt, ixfem, titr, lsubmodel, unitab)
Definition hm_read_fail_johnson.F:40
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
hm_option_read_mod
Definition hm_option_read_mod.F:412
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
submodel_mod
Definition submodel_mod.F:717
unitab_mod
Definition unitab_mod.F:517
second
real function second()
SECOND Using ETIME
Definition second_EXT_ETIME.f:35
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895
submodel_mod::submodel_data
Definition submodel_mod.F:728
unitab_mod::unit_type_
Definition unitab_mod.F:524