OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_fail_gene1.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_gene1 ../starter/source/materials/fail/gene1/hm_read_fail_gene1.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_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
31!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
32!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
33!||--- uses -----------------------------------------------------
34!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
35!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
36!|| message_mod ../starter/share/message_module/message_mod.F
37!|| submodel_mod ../starter/share/modules1/submodel_mod.F
38!||====================================================================
39 SUBROUTINE hm_read_fail_gene1(
40 . FAIL ,MAT_ID ,FAIL_ID ,IRUPT ,
41 . TITR ,LSUBMODEL,UNITAB ,FAIL_TAG )
42C-----------------------------------------------
43c ROUTINE DESCRIPTION :
44c Read generic tabulated failure model parameters
45C-----------------------------------------------
46C M o d u l e s
47C-----------------------------------------------
48 USE fail_param_mod
49 USE unitab_mod
50 USE message_mod
51 USE submodel_mod
52 USE hm_option_read_mod
53 USE elbuftag_mod
54C-----------------------------------------------
55C I m p l i c i t T y p e s
56C-----------------------------------------------
57#include "implicit_f.inc"
58C----------+---------+---+---+--------------------------------------------
59C C o m m o n B l o c k s
60C-----------------------------------------------
61#include "units_c.inc"
62C-----------------------------------------------
63C D u m m y A r g u m e n t s
64C-----------------------------------------------
65 INTEGER ,INTENT(IN) :: FAIL_ID !< failure model ID
66 INTEGER ,INTENT(IN) :: MAT_ID !< material law ID
67 INTEGER ,INTENT(IN) :: IRUPT !< failure model number
68 CHARACTER ,INTENT(IN) :: TITR*500 !< material model title
69 TYPE (UNIT_TYPE_) ,INTENT(IN) :: UNITAB !< table of input units
70 TYPE (SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(*) !< submodel table
71 TYPE (FAIL_PARAM_) ,INTENT(INOUT) :: FAIL !< failure model data structure
72 TYPE (FAIL_TAG_) ,INTENT(INOUT) :: FAIL_TAG !< failure model tag for buffer allocation
73C-----------------------------------------------
74C L o c a l V a r i a b l e s
75C-----------------------------------------------
76 INTEGER I,J,NSTEP,NCS,fct_IDsm,fct_IDps,fct_IDg12,fct_IDg13,
77 . fct_IDe1c,fct_IDel,Itab,tab_IDfld,Istrain,CRIT,
78 . ismooth,ncrit_max,nfunc,ntable,nuparam,nuvar,nmode,failip
79 INTEGER ,PARAMETER :: NSIZE = 6
80 integer ,DIMENSION(NSIZE) :: ifunc
81 my_real
82 . minpres,maxpres,sigp1,tmax,dtmin,epsdot_sm,sigvm,sigth,kf,
83 . epsdot_ps,maxeps,effeps,voleps,mineps,epssh,epsdot_fld,
84 . thin,volfrac,pthk,maxtemp,fscale_el,el_ref,freq_unit,
85 . pressure_unit,length_unit
86 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
87 CHARACTER(LEN=NCHARTITLE), DIMENSION(50) :: MODNAME
88C=======================================================================
89 is_encrypted = .false.
90 is_available = .false.
91 ! Initialization of the number of criteria
92 ncrit_max = 0
93 ! Initialization of the number of modes
94 nmode = 0
95 ! activated criteria
96 crit = 0
97C--------------------------------------------------
98C EXTRACT DATA (IS OPTION CRYPTED)
99C--------------------------------------------------
100 CALL hm_option_is_encrypted(is_encrypted)
101C======================================================================================
102C EXTRACT DATAS (REAL VALUES)
103C======================================================================================
104!---------------
105! -> Card1
106!---------------
107 CALL hm_get_floatv ('MAT_MINPRES' ,minpres ,is_available,lsubmodel,unitab)
108 IF (minpres /= zero) THEN
109 crit = ibset(crit,1)
110 nmode = nmode + 1
111 modname(nmode) = "Minimum hydrostatic pressure"
112 ENDIF
113 CALL hm_get_floatv ('MAT_MAXPRES' ,maxpres ,is_available,lsubmodel,unitab)
114 IF (maxpres /= zero) THEN
115 crit = ibset(crit,2)
116 nmode = nmode + 1
117 modname(nmode) = "Maximum hydrostatic pressure"
118 ENDIF
119 IF ((minpres /= zero).OR.(maxpres /= zero)) ncrit_max = ncrit_max + 1
120 CALL hm_get_floatv ('MAT_SIGP1' ,sigp1 ,is_available,lsubmodel,unitab)
121 IF (sigp1 /= zero) THEN
122 ncrit_max = ncrit_max + 1
123 crit = ibset(crit,3)
124 nmode = nmode + 1
125 modname(nmode) = "Maximum principal stress"
126 ENDIF
127 CALL hm_get_floatv ('MAT_TMAX' ,tmax ,is_available,lsubmodel,unitab)
128 IF (tmax /= zero) THEN
129 ncrit_max = ncrit_max + 1
130 crit = ibset(crit,4)
131 nmode = nmode + 1
132 modname(nmode) = "Failure time"
133 ENDIF
134 CALL hm_get_floatv ('MAT_DTMIN' ,dtmin ,is_available,lsubmodel,unitab)
135 IF (dtmin /= zero) THEN
136 ncrit_max = ncrit_max + 1
137 crit = ibset(crit,5)
138 nmode = nmode + 1
139 modname(nmode) = "Minimum timestep"
140 ENDIF
141!---------------
142! -> Card2
143!---------------
144 CALL hm_get_intv ('fct_IDsm' ,fct_idsm ,is_available,lsubmodel)
145 CALL hm_get_floatv ('MAT_EPSDSM' ,epsdot_sm ,is_available,lsubmodel,unitab)
146 CALL hm_get_floatv ('MAT_SIGVM' ,sigvm ,is_available,lsubmodel,unitab)
147 IF ((sigvm /= zero).OR.(fct_idsm /= 0)) THEN
148 ncrit_max = ncrit_max + 1
149 crit = ibset(crit,6)
150 nmode = nmode + 1
151 modname(nmode) = "Maximum equivalent stress"
152 ENDIF
153 CALL hm_get_floatv ('MAT_SIGTH' ,sigth ,is_available,lsubmodel,unitab)
154 CALL hm_get_floatv ('MAT_KF' ,kf ,is_available,lsubmodel,unitab)
155 IF ((sigth /= zero).OR.(kf /= zero)) THEN
156 ncrit_max = ncrit_max + 1
157 crit = ibset(crit,7)
158 nmode = nmode + 1
159 modname(nmode) = "Tuler-Butcher failure"
160 ENDIF
161!---------------
162! -> Card3
163!---------------
164 CALL hm_get_intv ('fct_IDps' ,fct_idps ,is_available,lsubmodel)
165 CALL hm_get_floatv ('MAT_EPSDPS' ,epsdot_ps ,is_available,lsubmodel,unitab)
166 CALL hm_get_floatv ('MAT_MAXEPS' ,maxeps ,is_available,lsubmodel,unitab)
167 IF ((maxeps /= zero).OR.(fct_idps /= 0)) THEN
168 ncrit_max = ncrit_max + 1
169 crit = ibset(crit,8)
170 nmode = nmode + 1
171 modname(nmode) = "Maximum principal strain"
172 ENDIF
173 CALL hm_get_floatv ('MAT_EFFEPS' ,effeps ,is_available,lsubmodel,unitab)
174 IF (effeps /= zero) THEN
175 ncrit_max = ncrit_max + 1
176 crit = ibset(crit,9)
177 nmode = nmode + 1
178 modname(nmode) = "Maximum effective strain"
179 ENDIF
180 CALL hm_get_floatv ('MAT_VOLEPS' ,voleps ,is_available,lsubmodel,unitab)
181 IF (voleps /= zero) THEN
182 ncrit_max = ncrit_max + 1
183 crit = ibset(crit,10)
184 nmode = nmode + 1
185 modname(nmode) = "Maximum volumetric strain"
186 ENDIF
187!---------------
188! -> Card4
189!---------------
190 CALL hm_get_floatv ('MAT_MINEPS' ,mineps ,is_available,lsubmodel,unitab)
191 IF (mineps /= zero) THEN
192 ncrit_max = ncrit_max + 1
193 crit = ibset(crit,11)
194 nmode = nmode + 1
195 modname(nmode) = "Minimum principal strain"
196 ENDIF
197 CALL hm_get_floatv ('MAT_EPSSH' ,epssh ,is_available,lsubmodel,unitab)
198 IF (epssh /= zero) THEN
199 ncrit_max = ncrit_max + 1
200 crit = ibset(crit,12)
201 nmode = nmode + 1
202 modname(nmode) = "Critical shear strain"
203 ENDIF
204 CALL hm_get_intv ('fct_IDg12' ,fct_idg12 ,is_available,lsubmodel)
205 IF (fct_idg12 /= 0) THEN
206 crit = ibset(crit,13)
207 nmode = nmode + 1
208 modname(nmode) = "Maximum in-plane shear strain"
209 ENDIF
210 CALL hm_get_intv ('fct_IDg13' ,fct_idg13 ,is_available,lsubmodel)
211 IF (fct_idg13 /= 0) THEN
212 crit = ibset(crit,14)
213 nmode = nmode + 1
214 modname(nmode) = "Maximum transversal shear strain"
215 ENDIF
216 CALL hm_get_intv ('fct_IDe1c' ,fct_ide1c ,is_available,lsubmodel)
217 IF (fct_ide1c /= 0) THEN
218 crit = ibset(crit,15)
219 nmode = nmode + 1
220 modname(nmode) = "Maximum in-plane major strain"
221 ENDIF
222 IF ((fct_idg13 /= 0).AND.(fct_ide1c /= 0)) THEN
223 ncrit_max = ncrit_max + 2
224 ELSEIF ((fct_idg12 /= 0).AND.((fct_idg13 /= 0).OR.(fct_ide1c /= 0))) THEN
225 ncrit_max = ncrit_max + 1
226 ELSEIF ((fct_idg12 == 0).AND.(fct_idg13 == 0).AND.(fct_ide1c == 0)) THEN
227 ELSE
228 ncrit_max = ncrit_max + 1
229 ENDIF
230!---------------
231! -> Card5
232!---------------
233 CALL hm_get_intv ('tab_IDfld' ,tab_idfld ,is_available,lsubmodel)
234 IF (tab_idfld /= 0) THEN
235 ncrit_max = ncrit_max + 1
236 crit = ibset(crit,16)
237 nmode = nmode + 1
238 modname(nmode) = "Forming Limit Diagram"
239 ENDIF
240 CALL hm_get_intv ('Itab' ,itab ,is_available,lsubmodel)
241 CALL hm_get_floatv ('MAT_EPSDFLD' ,epsdot_fld,is_available,lsubmodel,unitab)
242 CALL hm_get_intv ('MAT_NSTEP' ,nstep ,is_available,lsubmodel)
243 CALL hm_get_intv ('Ismooth' ,ismooth ,is_available,lsubmodel)
244 CALL hm_get_intv ('Istrain' ,istrain ,is_available,lsubmodel)
245 CALL hm_get_floatv ('MAT_THIN' ,thin ,is_available,lsubmodel,unitab)
246 IF (thin /= zero) THEN
247 ncrit_max = ncrit_max + 1
248 crit = ibset(crit,17)
249 nmode = nmode + 1
250 modname(nmode) = "Thinning failure"
251 ENDIF
252!---------------
253! -> Card6
254!---------------
255 CALL hm_get_intv ('FAILIP' ,failip ,is_available,lsubmodel)
256 IF (failip == 0) failip = 1
257 CALL hm_get_floatv ('MAT_VOLFRAC' ,volfrac ,is_available,lsubmodel,unitab)
258 CALL hm_get_floatv ('Pthk' ,pthk ,is_available,lsubmodel,unitab)
259 CALL hm_get_intv ('MAT_NCS' ,ncs ,is_available,lsubmodel)
260 CALL hm_get_floatv ('MAT_MAXTEMP' ,maxtemp ,is_available,lsubmodel,unitab)
261 IF (maxtemp /= zero) THEN
262 ncrit_max = ncrit_max + 1
263 crit = ibset(crit,18)
264 nmode = nmode + 1
265 modname(nmode) = "Maximum temperature"
266 ENDIF
267!---------------
268! -> Card7
269!---------------
270 CALL hm_get_intv ('fct_IDel' ,fct_idel ,is_available,lsubmodel)
271 CALL hm_get_floatv ('Fscale_el' ,fscale_el ,is_available,lsubmodel,unitab)
272 CALL hm_get_floatv ('El_ref' ,el_ref ,is_available,lsubmodel,unitab)
273C======================================================================================
274C DEFAULT VALUES
275C======================================================================================
276 ! Card 1
277 minpres = -abs(minpres)
278 IF (minpres == zero) minpres = -infinity
279 maxpres = abs(maxpres)
280 IF (maxpres == zero) maxpres = infinity
281 IF (sigp1 == zero) sigp1 = infinity
282 IF (tmax == zero) tmax = infinity
283 IF (dtmin == zero) dtmin = -infinity
284 ! Card 2
285 IF (fct_idsm > 0) THEN
286 IF (epsdot_sm == zero) THEN
287 CALL hm_get_floatv_dim('MAT_EPSDSM' ,freq_unit,is_available, lsubmodel, unitab)
288 epsdot_sm = one*freq_unit
289 ENDIF
290 IF (sigvm == zero) THEN
291 CALL hm_get_floatv_dim('MAT_SIGVM' ,pressure_unit,is_available, lsubmodel, unitab)
292 sigvm = one*pressure_unit
293 ENDIF
294 ELSE
295 IF (sigvm == zero) sigvm = infinity
296 ENDIF
297 IF (sigth == zero) sigth = infinity
298 IF (kf == zero) kf = infinity
299 ! Card 3
300 IF (fct_idps > 0) THEN
301 IF (epsdot_ps == zero) THEN
302 CALL hm_get_floatv_dim('MAT_EPSDPS' ,freq_unit,is_available, lsubmodel, unitab)
303 epsdot_ps = one*freq_unit
304 ENDIF
305 IF (maxeps == zero) THEN
306 maxeps = one
307 ENDIF
308 ELSE
309 IF (maxeps == zero) maxeps = two*infinity
310 ENDIF
311 IF (effeps == zero) effeps = infinity
312 IF (voleps == zero) voleps = infinity
313 ! Card 4
314 mineps = -abs(mineps)
315 IF (mineps == zero) mineps = -infinity
316 IF (epssh == zero) epssh = infinity
317 ! Card 5
318 IF (tab_idfld > 0) THEN
319 IF (epsdot_fld == zero) THEN
320 CALL hm_get_floatv_dim('MAT_EPSDFLD' ,freq_unit,is_available, lsubmodel, unitab)
321 epsdot_fld = one*freq_unit
322 ENDIF
323 ENDIF
324 IF (nstep == 0) nstep = 10
325 IF (ismooth == 0) ismooth = 1
326 IF (thin == zero) thin = -infinity
327 ! Card 6
328 IF (volfrac == zero) volfrac = half
329 volfrac = min(volfrac,one)
330 volfrac = max(zero,volfrac)
331 IF (pthk == zero) pthk = one-em06
332 pthk = min(pthk, one)
333 pthk = max(pthk,-one)
334 IF (ncs == 0 ) ncs = 1
335 IF (ncs > ncrit_max) THEN
336 IF (ncrit_max > 0) THEN
337 CALL ancmsg(msgid=2043,msgtype=msgwarning,
338 . anmode=aninfo_blind_1,i1=mat_id,
339 . c1=titr,i2=ncs,i3=ncrit_max,i4=ncrit_max)
340 ncs = ncrit_max
341 ELSE
342 CALL ancmsg(msgid=3040,msgtype=msgwarning,
343 . anmode=aninfo_blind_1,i1=mat_id,
344 . c1=titr)
345 ENDIF
346 ENDIF
347 ncs = max(1,ncs)
348 IF (maxtemp == zero) maxtemp = infinity
349 ! Card 7
350 IF (fscale_el == zero) fscale_el = one
351 IF (((fct_idel > 0).OR.(fct_idg12 > 0).OR.(fct_idg13 > 0).OR.(fct_ide1c > 0)
352 . .OR.(tab_idfld > 0.AND.itab == 2)).AND.(el_ref == zero)) THEN
353 CALL hm_get_floatv_dim('El_ref' ,length_unit,is_available, lsubmodel, unitab)
354 el_ref = one*length_unit
355 ENDIF
356C======================================================================================
357c Filling buffer tables
358c-----------------------------------------------------
359 nuvar = 8
360 nuparam = 23
361 nfunc = 6
362 ntable = 0
363 IF (tab_idfld > 0) ntable = 1
364c
365 fail%KEYWORD = 'GENE1'
366 fail%IRUPT = irupt
367 fail%FAIL_ID = fail_id
368 fail%NUPARAM = nuparam
369 fail%NIPARAM = 6
370 fail%NUVAR = nuvar
371 fail%NFUNC = nfunc
372 fail%NTABLE = ntable
373 fail%NMOD = nmode
374 fail%PTHK = pthk
375c
376 ALLOCATE (fail%UPARAM(fail%NUPARAM))
377 ALLOCATE (fail%IPARAM(fail%NIPARAM))
378 ALLOCATE (fail%IFUNC (fail%NFUNC))
379 ALLOCATE (fail%TABLE (fail%NTABLE))
380c
381 ! Modes of failure
382 fail_tag%LF_DAMMX = fail_tag%LF_DAMMX + fail%NMOD
383 ALLOCATE (fail%MODE(fail%NMOD))
384 DO i = 1,fail%NMOD
385 fail%MODE(i) = modname(i)
386 ENDDO
387c
388 fail%IFUNC(1) = fct_idsm
389 fail%IFUNC(2) = fct_idps
390 fail%IFUNC(3) = fct_idg12
391 fail%IFUNC(4) = fct_idg13
392 fail%IFUNC(5) = fct_ide1c
393 fail%IFUNC(6) = fct_idel
394 IF (ntable == 1) fail%TABLE(1) = tab_idfld
395c
396 fail%IPARAM(1) = crit
397 fail%IPARAM(2) = itab
398 fail%IPARAM(3) = nstep
399 fail%IPARAM(4) = ncs
400 fail%IPARAM(5) = ismooth
401 fail%IPARAM(6) = istrain
402c
403 fail%UPARAM(1) = minpres
404 fail%UPARAM(2) = maxpres
405 fail%UPARAM(3) = sigp1
406 fail%UPARAM(4) = tmax
407 fail%UPARAM(5) = dtmin
408 fail%UPARAM(6) = epsdot_sm
409 fail%UPARAM(7) = sigvm
410 fail%UPARAM(8) = sigth
411 fail%UPARAM(9) = kf
412 fail%UPARAM(10) = epsdot_ps
413 fail%UPARAM(11) = maxeps
414 fail%UPARAM(12) = effeps
415 fail%UPARAM(13) = voleps
416 fail%UPARAM(14) = mineps
417 fail%UPARAM(15) = epssh
418 fail%UPARAM(16) = epsdot_fld
419 fail%UPARAM(17) = thin
420 fail%UPARAM(18) = volfrac
421 fail%UPARAM(19) = pthk
422 fail%UPARAM(20) = maxtemp
423 fail%UPARAM(21) = fscale_el
424 fail%UPARAM(22) = el_ref
425 fail%UPARAM(23)= failip
426c--------------------------
427c Printout data
428c--------------------------
429 IF (is_encrypted) THEN
430 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
431 ELSE
432 ! Title of the failure criterion
433 WRITE(iout,1000)
434 ! Number of criteria setup in the input
435 WRITE(iout,1001) ncrit_max
436 ! Number of criteria to reach for element deletion
437 WRITE(iout,1002) ncs
438 ! Definition of the criteria
439 WRITE(iout,1003)
440 ! Card 1
441 IF (btest(crit,1)) WRITE(iout,1100) minpres
442 IF (btest(crit,2)) WRITE(iout,1101) maxpres
443 IF (btest(crit,1) .OR. btest(crit,2)) WRITE(iout,1105)
444 IF (btest(crit,3)) WRITE(iout,1102) sigp1
445 IF (btest(crit,4)) WRITE(iout,1103) tmax
446 IF (btest(crit,5)) WRITE(iout,1104) dtmin
447 ! Card 2
448 IF (btest(crit,6)) THEN
449 IF (fct_idsm > 0) THEN
450 WRITE(iout,2000) fct_idsm,epsdot_sm,sigvm
451 ELSE
452 WRITE(iout,2100) sigvm
453 ENDIF
454 ENDIF
455 IF (btest(crit,7)) WRITE(iout,2001) sigth,kf
456 ! Card 3
457 IF (btest(crit,8)) THEN
458 IF (fct_idps > 0) THEN
459 WRITE(iout,3000) fct_idps,epsdot_ps,maxeps
460 ELSE
461 WRITE(iout,3100) maxeps
462 ENDIF
463 ENDIF
464 IF (btest(crit,9)) WRITE(iout,3001) effeps
465 IF (btest(crit,10)) WRITE(iout,3002) voleps
466 ! Card 4
467 IF (btest(crit,11)) WRITE(iout,4000) mineps
468 IF (btest(crit,12)) WRITE(iout,4001) epssh
469 IF (btest(crit,13)) WRITE(iout,4100) fct_idg12
470 IF (btest(crit,14)) WRITE(iout,4200) fct_idg13
471 IF (btest(crit,15)) WRITE(iout,4300) fct_ide1c
472 ! Card 5
473 IF (btest(crit,16)) THEN
474 WRITE(iout,5000) tab_idfld,itab,epsdot_fld,
475 . ismooth,istrain
476 ENDIF
477 IF (btest(crit,17)) WRITE(iout,5002) thin
478 ! Card 6
479 IF (btest(crit,18)) WRITE(iout,6001) maxtemp
480 ! Card 7
481 IF (fct_idel > 0) WRITE(iout,7000) fct_idel,fscale_el
482 IF (((fct_idel > 0).OR.(fct_idg12 > 0).OR.(fct_idg13 > 0).OR.(fct_ide1c > 0)
483 . .OR.(tab_idfld > 0.AND.itab == 2))) WRITE(iout,7100) el_ref
484 ! Element deletion properties
485 WRITE(iout,7101)
486 WRITE(iout,6000) volfrac,pthk,failip
487 WRITE(iout,5001) nstep
488 WRITE(iout,8000)
489 ENDIF
490c-----------------------------------------------------
491 1000 FORMAT(
492 & 5x,' ----------------------------------------------------',/,
493 & 5x,' FAILURE CRITERION : GENE1 ',/,
494 & 5x,' ----------------------------------------------------',/)
495 1001 FORMAT(
496 & 5x,'NUMBER OF CRITERIA SPECIFIED IN THE INPUT. . . . . .=',i10,/)
497 1002 FORMAT(
498 & 5x,'NUMBER OF CRITERIA TO REACH FOR ELEMENT DELETION . .=',i10,/)
499 1003 FORMAT(
500 & 5x,'SPECIFIED ELEMENT DELETION CRITERIA:',/,
501 & 5x,'------------------------------------',/)
502 1100 FORMAT(
503 & 5x,'MINIMUM PRESSURE . . . . . . . . . . . . . . . . . .=',1pg20.13,/)
504 1101 FORMAT(
505 & 5x,'MAXIMUM PRESSURE . . . . . . . . . . . . . . . . . .=',1pg20.13,/)
506 1102 FORMAT(
507 & 5x,'MAXIMUM PRINCIPAL STRESS . . . . . . . . . . . . . .=',1pg20.13,/
508 & 5x,' < 0.0 : RESTRICTED TO POSITIVE STRESS TRIAXIALITIES',/,
509 & 5x,' > 0.0 : UNRESTRICTED ',/)
510 1103 FORMAT(
511 & 5x,'FAILURE TIME . . . . . . . . . . . . . . . . . . . .=',1pg20.13,/)
512 1104 FORMAT(
513 & 5x,'MINIMUM TIME STEP . . . . . . . . . . . . . . . . .=',1pg20.13,/)
514 1105 FORMAT(
515 & 5x,'NOTE : PRESSURE IS DEFINED SO THAT IT IS POSITIVE IN COMPRESSION',/)
516 2000 FORMAT(
517 & 5x,'FUNCTION ID FOR EQ. STRESS VS STRAIN-RATE . . . . .=',i10,/,
518 & 5x,'REFERENCE STRAIN-RATE . . . . . . . . . . . . . . .=',1pg20.13,/
519 & 5x,'EQUIVALENT STRESS SCALE FACTOR . . . . . . . . . . .=',1pg20.13,/)
520 2001 FORMAT(
521 & 5x,'INITIAL STRESS FOR TULER-BUTCHER CRITERION . . . . .=',1pg20.13,/
522 & 5x,'CRITICAL TULER-BUTCHER INTEGRAL VALUE . . . . . . .=',1pg20.13,/)
523 2100 FORMAT(
524 & 5x,'MAXIMUM EQUIVALENT STRESS . . . . . . . . . . . . .=',1pg20.13,/)
525 3000 FORMAT(
526 & 5x,'FUNCTION ID FOR PRINC. STRAIN VS STRAIN-RATE . . . .=',i10,/,
527 & 5x,'REFERENCE STRAIN-RATE . . . . . . . . . . . . . . .=',1pg20.13,/
528 & 5x,'PRINCIPAL STRAIN SCALE FACTOR . . . . . . . . . . .=',1pg20.13,/)
529 3001 FORMAT(
530 & 5x,'MAXIMUM EFFECTIVE STRAIN . . . . . . . . . . . . . .=',1pg20.13,/)
531 3002 FORMAT(
532 & 5x,'MAXIMUM VOLUMETRIC STRAIN . . . . . . . . . . . . .=',1pg20.13,/)
533 3100 FORMAT(
534 & 5x,'MAXIMUM PRINCIPAL STRAIN . . . . . . . . . . . . . .=',1pg20.13,/)
535 4000 FORMAT(
536 & 5x,'MINIMUM PRINCIPAL STRAIN . . . . . . . . . . . . . .=',1pg20.13,/)
537 4001 FORMAT(
538 & 5x,'TENSORIAL SHEAR STRAIN AT FAILURE. . . . . . . . . .=',1pg20.13,/)
539 4100 FORMAT(
540 & 5x,'FUNCTION ID FOR MAX.IN-PLANE SHEAR VS EL.SIZE . . .=',i10,/)
541 4200 FORMAT(
542 & 5x,'FUNCTION ID FOR MAX.TRANSVERSE SHEAR VS EL.SIZE . .=',i10,/)
543 4300 FORMAT(
544 & 5x,'FUNCTION ID FOR IN-PLANE MAJOR STRAIN VS EL.SIZE . .=',i10,/)
545 5000 FORMAT(
546 & 5x,'TABLE ID FOR FORMING LIMIT DIAGRAM . . . . . . . . .=',i10,/,
547 & 5x,'flag for table dependency . . . . . . . . . . . . .=',I10,/,
548 & 5X,' = 1 : forming limit diagram vs strain-rate ',/,
549 & 5X,' = 2 : forming limit diagram vs element SIZE ',/,
550 & 5X,'reference strain-rate . . . . . . . . . . . . . . .=',1PG20.13,/,
551 & 5X,'tabulated FUNCTION interpolation flag . . . . . . .=',I10,/,
552 & 5X,'engineering / true input strain flag . . . . . . . .=',I10,/)
553 5001 FORMAT(
554 & 5X,'number of cycles before element deletion . . . . . .=',I10,/)
555 5002 FORMAT(
556 & 5X,'shell thinning at failure . . . . . . . . . . . . .=',1PG20.13,/)
557 6000 FORMAT(
558 & 5X,'volume fraction for solid element deletion . . . . .=',1PG20.13,/,
559 & 5X,'shell element deletion parameter pthickfail. . . . .=',1PG20.13,/,
560 & 5X,' > 0.0 : fraction of failed thickness ',/,
561 & 5X,' < 0.0 : fraction of failed intg. points or layers ',/,
562 & 5X,'number of failed intg. points prior to elem deletion .=',I10/)
563 6001 FORMAT(
564 & 5X,'maximum temperature . . . . . . . . . . . . . . . .=',1PG20.13,/)
565 7000 FORMAT(
566 & 5X,'regularization function id . . . . . . . . . . . . .=',I10,/,
567 & 5X,'scale factor for regularization function . . . . . .=',1PG20.13,/)
568 7100 FORMAT(
569 & 5X,'reference element length . . . . . . . . . . . . . .=',1PG20.13,/)
570 7101 FORMAT(
571 & 5X,'element deletion properties:',/,
572 & 5X,'----------------------------',/)
573 8000 FORMAT(
574 & 5X,' ----------------------------------------------------',//)
575c-----------
576 END
my_real
#define my_real
Definition cppsort.cpp:32
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
hm_get_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:159
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:284
hm_read_fail_gene1
subroutine hm_read_fail_gene1(fail, mat_id, fail_id, irupt, titr, lsubmodel, unitab, fail_tag)
Definition hm_read_fail_gene1.F:42
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
id
initmumps id
Definition multiplerhs_example.m:4
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
elbuftag_mod
Definition elbuftag_mod.F:227
hm_option_read_mod
Definition hm_option_read_mod.F:408
message_mod
Definition message_mod.F:1249
submodel_mod
Definition submodel_mod.F:712
unitab_mod
Definition unitab_mod.F:513
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:889
tabulated
subroutine tabulated(iflag, nel, pm, off, eint, mu, espe, dvol, df, vnew, mat, psh, pnew, dpdm, dpde, npf, tf)
Definition tabulated.F:32