OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_fail_tab1.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_tab1 ../starter/source/materials/fail/tabulated/hm_read_fail_tab1.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!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
35!|| message_mod ../starter/share/message_module/message_mod.F
36!|| submodel_mod ../starter/share/modules1/submodel_mod.F
37!|| table_mod ../starter/share/modules1/table_mod.F
38!||====================================================================
39 SUBROUTINE hm_read_fail_tab1(FAIL ,
40 . MAT_ID ,FAIL_ID ,IRUPT ,IXFEM ,
41 . LSUBMODEL,UNITAB )
42C-----------------------------------------------
43C ROUTINE DESCRIPTION :
44C ===================
45C READ TABULATED FAILURE MODEL (/FAIL/TAB1)
46C-----------------------------------------------
47C M o d u l e s
48C-----------------------------------------------
49 USE fail_param_mod
50 USE unitab_mod
51 USE message_mod
52 USE submodel_mod
53 USE hm_option_read_mod
54 USE table_mod
55C-----------------------------------------------
56C I m p l i c i t T y p e s
57C-----------------------------------------------
58#include "implicit_f.inc"
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62#include "units_c.inc"
63C-----------------------------------------------
64C D u m m y A r g u m e n t s
65C-----------------------------------------------
66 INTEGER ,INTENT(IN) :: FAIL_ID ! failure model ID
67 INTEGER ,INTENT(IN) :: MAT_ID ! material law ID
68 INTEGER ,INTENT(IN) :: IRUPT ! failure model number
69 TYPE (UNIT_TYPE_) ,INTENT(IN) :: UNITAB ! table of input units
70 TYPE (SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(*)
71 INTEGER ,INTENT(INOUT) :: IXFEM ! XFEM activation flag
72 TYPE (FAIL_PARAM_) ,INTENT(INOUT) :: FAIL
73C-----------------------------------------------
74C L o c a l V a r i a b l e s
75C-----------------------------------------------
76 INTEGER :: IFAIL_SH,ISOLID,DMG_FLAG,INST_FLAG,
77 . IFUN_DMG,ITAB_EPSF,ITAB_INST,IFUN_SIZE,IFUN_TEMP
78 my_real :: p_thick,pthkf,p_thinnfail,scale_temp,scale_el,el_ref,
79 . y1scale,x1scale,y2scale,x2scale,dcrit,dd,dn,dadv,
80 . ecrit,fade_expo,fscal_unit,shrf,biaxf
81C-----------------------------------------------
82 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
83C-----------------------------------------------
84c UVAR storage:
85C 1 = DAMAGE
86C 2 = initial shell thickness
87C 3 = DCrit_NS --> instability starts
88C 4 = percent from instability to failure
89C 5 = initial characteristic el. length
90C 6 = IPOS 1 for Table
91C 7 = IPOS 2 for Table
92C 8 = IPOS 1 for vinter
93C=======================================================================
94 is_encrypted = .false.
95 is_available = .false.
96 dadv = zero
97 ecrit = zero
98C--------------------------------------------------
99C EXTRACT DATA (IS OPTION CRYPTED)
100C--------------------------------------------------
101 CALL hm_option_is_encrypted(is_encrypted)
102C--------------------------------------------------
103C EXTRACT INPUT DATA
104C--------------------------------------------------
105Card1
106 CALL hm_get_intv ('Ifail_sh' ,ifail_sh ,is_available,lsubmodel)
107 CALL hm_get_intv ('Ifail_so' ,isolid ,is_available,lsubmodel)
108 CALL hm_get_floatv ('P_thickfail' ,p_thick ,is_available,lsubmodel,unitab)
109 CALL hm_get_floatv ('P_thinfail' ,p_thinnfail ,is_available,lsubmodel,unitab)
110 CALL hm_get_intv ('Ixfem' ,ixfem ,is_available,lsubmodel)
111Card2
112 CALL hm_get_floatv ('Dcrit' ,dcrit ,is_available,lsubmodel,unitab)
113 CALL hm_get_floatv ('D' ,dd ,is_available,lsubmodel,unitab)
114 CALL hm_get_floatv ('n' ,dn ,is_available,lsubmodel,unitab)
115 CALL hm_get_floatv ('Dadv' ,dadv ,is_available,lsubmodel,unitab)
116 CALL hm_get_intv ('fct_IDd' ,ifun_dmg ,is_available,lsubmodel)
117Card3
118 CALL hm_get_intv ('table1_ID' ,itab_epsf ,is_available,lsubmodel)
119 CALL hm_get_floatv ('Xscale1' ,y1scale ,is_available,lsubmodel,unitab)
120 CALL hm_get_floatv ('Xscale2' ,x1scale ,is_available,lsubmodel,unitab)
121 CALL hm_get_intv ('table2_ID' ,itab_inst ,is_available,lsubmodel)
122 CALL hm_get_floatv ('Xscale3' ,y2scale ,is_available,lsubmodel,unitab)
123 CALL hm_get_floatv ('Xscale4' ,x2scale ,is_available,lsubmodel,unitab)
124Card4
125 CALL hm_get_intv ('fct_IDel' ,ifun_size ,is_available,lsubmodel)
126 CALL hm_get_floatv ('Fscale_el' ,scale_el ,is_available,lsubmodel,unitab)
127 CALL hm_get_floatv ('EI_ref' ,el_ref ,is_available,lsubmodel,unitab)
128 CALL hm_get_floatv ('Inst_start' ,ecrit ,is_available,lsubmodel,unitab)
129 CALL hm_get_floatv ('Fad_exp' ,fade_expo ,is_available,lsubmodel,unitab)
130 CALL hm_get_intv ('Ch_i_f' ,inst_flag ,is_available,lsubmodel)
131Card5
132 CALL hm_get_intv ('fct_IDt' ,ifun_temp ,is_available,lsubmodel)
133 CALL hm_get_floatv ('FscaleT' ,scale_temp ,is_available,lsubmodel,unitab)
134 CALL hm_get_floatv ('Shear_limit' ,shrf ,is_available,lsubmodel,unitab)
135 CALL hm_get_floatv ('Biax_limit' ,biaxf ,is_available,lsubmodel,unitab)
136c-----------------------------------------------------------------------
137 ! Error massage: 'table1_ID' is mandatory:
138 IF (itab_epsf == 0) THEN
139 CALL ancmsg(msgid=2068, msgtype=msgerror, anmode=aninfo_blind,
140 . i1=mat_id)
141 ENDIF
142c-----------------------------------------------------------------------
143c Set default parameter values
144c-----------------------------------------------------------------------
145 IF (dcrit == zero) dcrit = one
146 IF (dadv == zero) dadv = dcrit
147 IF (dadv > dcrit) THEN
148 dadv = dcrit
149 CALL ancmsg(msgid=974, msgtype=msgwarning, anmode=aninfo,
150 . i1=mat_id)
151 ENDIF
152 IF (dd == one ) dd = 0.999
153 IF (itab_inst > 0) THEN
154 ecrit = zero
155 ELSEIF (ecrit == zero) THEN
156 ecrit = dd
157 ENDIF
158 IF (dn == zero) dn = one
159 IF (ifail_sh == 0) ifail_sh = 1
160 IF (isolid == 0) isolid = 1
161 IF (ixfem /= 1 .AND. ixfem /= 2) ixfem = 0
162 IF (ixfem > 0) isolid = 0
163 IF (y1scale == zero) y1scale = one
164 IF (y2scale == zero) y2scale = one
165 IF (scale_el == zero) scale_el = one
166 IF (scale_temp == zero) scale_temp = one
167 IF (shrf == zero) shrf =-one
168 IF (biaxf == zero) biaxf = one
169c
170 IF (fade_expo > zero .or. ecrit /= zero) THEN
171 dmg_flag = 1
172 ELSE
173 dmg_flag = 0
174 ENDIF
175 IF (inst_flag == 0 .OR. inst_flag > 3) inst_flag = 1
176 IF (p_thick == one ) p_thick = p_thick - em06
177 IF (p_thick == zero) p_thick = one-em06
178 p_thick = min(p_thick, one)
179 p_thick = max(p_thick,-one)
180c---------------------------
181 IF (p_thick > zero .and. ifail_sh > 1) THEN
182 pthkf = p_thick
183 ELSEIF (ifail_sh == 1) THEN
184 pthkf = em06
185 ELSEIF (ifail_sh == 2) THEN
186 pthkf = one - em06
187 ENDIF
188c-----------------------------------------------------------------------
189 CALL hm_get_floatv_dim('Xscale2' ,fscal_unit ,is_available ,lsubmodel ,unitab)
190 IF (x1scale == zero) x1scale = one*fscal_unit
191c
192 CALL hm_get_floatv_dim('Xscale4' ,fscal_unit ,is_available ,lsubmodel ,unitab)
193 IF (x2scale == zero) x2scale = one*fscal_unit
194c
195 CALL hm_get_floatv_dim('EI_ref' ,fscal_unit ,is_available ,lsubmodel ,unitab)
196 IF (el_ref == zero) el_ref = one*fscal_unit
197c-----------------------------------------------------------------------
198 fail%KEYWORD = 'TAB1'
199 fail%IRUPT = irupt
200 fail%FAIL_ID = fail_id
201 fail%NUPARAM = 22
202 fail%NIPARAM = 0
203 fail%NUVAR = 8
204 fail%NFUNC = 4
205 fail%NTABLE = 2
206 fail%NMOD = 0
207 fail%PTHK = pthkf
208c
209 ALLOCATE (fail%UPARAM(fail%NUPARAM))
210 ALLOCATE (fail%IPARAM(fail%NIPARAM))
211 ALLOCATE (fail%IFUNC (fail%NFUNC))
212 ALLOCATE (fail%TABLE (fail%NTABLE))
213c
214 fail%UPARAM(1) = isolid
215 fail%UPARAM(2) = ifail_sh
216 fail%UPARAM(3) = 0 ! not used (P_THICK)
217 fail%UPARAM(4) = dcrit
218 fail%UPARAM(5) = dd
219 fail%UPARAM(6) = dn
220 fail%UPARAM(7) = scale_temp
221 fail%UPARAM(8) = scale_el
222 fail%UPARAM(9) = el_ref
223 fail%UPARAM(10)= 0 ! not used (xfem)
224 fail%UPARAM(11)= dadv
225 fail%UPARAM(12)= y1scale
226 fail%UPARAM(13)= one / x1scale
227 fail%UPARAM(14)= y2scale
228 fail%UPARAM(15)= one / x2scale
229 fail%UPARAM(16)= p_thinnfail
230 fail%UPARAM(17)= ecrit
231 fail%UPARAM(18)= fade_expo
232 fail%UPARAM(19)= dmg_flag
233 fail%UPARAM(20)= inst_flag
234 fail%UPARAM(21)= shrf
235 fail%UPARAM(22)= biaxf
236c
237 fail%TABLE(1) = itab_epsf
238 fail%TABLE(2) = itab_inst
239 fail%IFUNC(1) = ifun_size
240 fail%IFUNC(2) = ifun_temp
241 fail%IFUNC(3) = ifun_dmg
242 IF (fade_expo < zero) THEN
243 fail%IFUNC(4) = int(abs(fade_expo))
244 ELSE
245 fail%IFUNC(4) = 0
246 ENDIF
247c-----------------------------------------------------------------------
248 IF (is_encrypted)THEN
249 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
250 ELSE
251 WRITE(iout,100) mat_id,irupt,fail_id
252
253 WRITE(iout, 1001) itab_epsf, y1scale, x1scale
254 IF (itab_inst /= 0) THEN
255 WRITE(iout, 1002) itab_inst, y2scale, x2scale
256 ENDIF
257 IF (ixfem > 0) WRITE(iout, 1003) ixfem,dadv
258 WRITE(iout, 1004) p_thick,p_thinnfail
259 IF (ifun_dmg > 0) THEN
260 WRITE(iout, 1009) dcrit,ifun_dmg,ecrit
261 ELSE
262 WRITE(iout, 1005) dcrit,dd,dn,ecrit
263 ENDIF
264 IF (fade_expo >= zero) THEN
265 WRITE(iout, 1006) fade_expo
266 ELSE
267 WRITE(iout, 1007) int(abs(fade_expo))
268 ENDIF
269 WRITE(iout, 1008) dmg_flag,
270 . ifun_temp,scale_temp,
271 . ifun_size,scale_el,el_ref,shrf,biaxf,inst_flag
272 WRITE(iout, 1008) dmg_flag,ifun_temp,scale_temp,ifun_size,scale_el,
273 . el_ref,shrf,biaxf,inst_flag
274c for shell
275 IF (ixfem == 0) THEN
276 IF(ifail_sh == 1) THEN
277 WRITE(iout, 1100)
278 ELSEIF (ifail_sh == 2) THEN
279 WRITE(iout, 1200)
280 ELSEIF (ifail_sh == 3) THEN
281 WRITE(iout, 1300)
282 ENDIF
283 ELSE IF (ixfem == 1) THEN
284 WRITE(iout, 1400)
285 END IF
286c for solid
287 IF (isolid == 1) THEN
288 WRITE(iout, 2100)
289 ELSEIF(isolid == 2) THEN
290 WRITE(iout, 2200)
291 ENDIF
292C
293 ENDIF ! IS_ENCRYPTED
294c-----------
295 RETURN
296c-----------------------------------------------------------------------
297 100 FORMAT(//
298 & 5x,'MAT_ID . . . . . . . . . . . . . . .=',i10/
299 & 5x,'FAILURE MODEL. . . . . . . . . . . .=',i10/
300 & 5x,'FAIL_ID. . . . . . . . . . . . . . .=',i10/)
301 1001 FORMAT(
302 & 5x,'TABULATED FAILURE CRITERIA WITH DAMAGE',//,
303 & 5x,'STRAIN TABLE ID . . . . . . . . . . . . . . . . . . .=',i10/
304 & 5x,' SCALE FACTOR FOR FAILURE STRAIN . . . . . . . . .=',1pg20.13/
305 & 5x,' STRAIN RATE FACTOR FOR FAILURE STRAIN . . . . . .=',1pg20.13)
306 1002 FORMAT(
307 & 5x,'NECKING TABLE ID . . . . . . . . . . . . . . . .=',i10/
308 & 5x,' SCALE FACTOR FOR NECKING STRAIN. . . . . . . . .=',1pg20.13/
309 & 5x,' STRAIN RATE FACTOR FOR NECKING STRAIN. . . . . .=',1pg20.13)
310 1003 FORMAT(
311 & 5x,'FLAG XFEM. . . . . . . .. . . . . . . . . . . . . . .=',i10/
312 & 5x,' XFEM ADVANCEMENT SOFTENING FACTOR. . . . . . . .=',1pg20.13)
313 1004 FORMAT(
314 & 5x,'SHELL ELEMENT DELETION PARAMETER PTHICKFAIL . . . . .=',1pg20.13,/,
315 & 5x,' > 0.0 : FRACTION OF FAILED THICKNESS ',/,
316 & 5x,' < 0.0 : FRACTION OF FAILED INTG. POINTS OR LAYERS',/,
317 & 5x,'SHELL FAILURE DUE TO THINNING . . . . . . . . . . . .=',1pg20.13)
318 1005 FORMAT(
319 & 5x,'CRITICAL DAMAGE VALUE . . . . . . . . . . . . . . . .=',1pg20.13/
320 & 5x,'DAMAGE PARAMETER D. . . . . . . . . . . . . . . . . .=',1pg20.13/
321 & 5x,'DAMAGE PARAMETER N. . . . . . . . . . . . . . . . . .=',1pg20.13/
322 & 5x,'INSTABILITY STRAIN. . . . . . . . . . . . . . . . . .=',1pg20.13)
323 1009 FORMAT(
324 & 5x,'CRITICAL DAMAGE VALUE . . . . . . . . . . . . . . . .=',1pg20.13/
325 & 5x,'DAMAGE EVOLUTION FUNCTION . . . . . . . . . . . . . .=',i10/
326 & 5x,'INSTABILITY STRAIN. . . . . . . . . . . . . . . . . .=',1pg20.13)
327 1006 FORMAT(
328 & 5x,'FADE PARAMETER. . . . . . . . . . . . . . . . . . . .=',1pg20.13)
329 1007 FORMAT(
330 & 5x,'FADE FUNCTION . . . . . . . . . . . . . . . . . . . .=',i10)
331 1008 FORMAT(
332 & 5x,'DAMAGE FLAG . . . . . . . . . . . . . . . . . . . . .=',i10/
333 & 5x,'TEMPERATURE SCALE FUNCTION. . . . . . . . . . . . . .=',i10/
334 & 5x,'SCALE FACTOR OF TEMPERATURE FUNCTION. . . . . . . . .=',1pg20.13/
335 & 5x,'ELEMENT LENGTH FUNCTION . . . . . . . . . . . . . . .=',i10/
336 & 5x,'SCALE FACTOR OF LENGTH FUNCTION . . . . . . . . . . .=',1pg20.13/
337 & 5x,'REFERENCE ELEMENT LENGTH. . . . . . . . . . . . . . .=',1pg20.13/
338 & 5x,'SHEAR TRIAXIALITY LIMIT FOR ELEMENT SIZE SCALING. . .=',1pg20.13/
339 & 5x,'BI-TRACT TRIAXIALITY LIMIT FOR ELEMENT SIZE SCALING .=',1pg20.13/
340 & 5x,'REGULARIZATION FLAG . . . . . . . . . . . . . . . . .=',i10/
341 & 5x,'FAILURE OPTION:')
342 1100 FORMAT(
343 & 5x,' SHELL ELEMENT DELETION AFTER FAILURE OF ONE LAYER')
344 1200 FORMAT(
345 & 5x,' STRESS TENSOR IN SHELL LAYER SET TO ZERO AFTER FAILURE')
346 1300 FORMAT(
347 & 5x,' SHELL ELEMENT DELETION AFTER FAILURE OF ALL LAYERS')
348 1400 FORMAT(
349 & 5x,' SHELL ELEMENT CRACKING AFTER FAILURE')
350 2100 FORMAT(
351 & 5x,' SOLID ELEMENT DELETION AFTER FAILURE')
352 2200 FORMAT(
353 & 5x,' DEVIATORIC STRESS IN SOLID WILL VANISH AFTER FAILURE')
354c-----------------------------------------------------------------------
355 END
my_real
#define my_real
Definition cppsort.cpp:32
ecrit
subroutine ecrit(timers, partsav, ms, v, in, r, dmas, weight, enintot, ekintot, a, ar, fxbipm, fxbrpm, monvol, xmom_sms, sensors, qfricint, ipari, weight_md, wfexth, iflag, ms_2d, multi_fvm, mas_nd, kend, h3d_data, dynain_data, usreint, output)
Definition ecrit.F:51
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418
hm_get_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:161
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_tab1
subroutine hm_read_fail_tab1(fail, mat_id, fail_id, irupt, ixfem, lsubmodel, unitab)
Definition hm_read_fail_tab1.F:42
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
submodel_mod
Definition submodel_mod.F:717
table_mod
Definition table_mod.F:113
unitab_mod
Definition unitab_mod.F:517
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