OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat04.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_mat04 ../starter/source/materials/mat/mat004/hm_read_mat04.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_mat ../starter/source/materials/mat/hm_read_mat.F90
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!|| init_mat_keyword ../starter/source/materials/mat/init_mat_keyword.F
33!||--- uses -----------------------------------------------------
34!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.f
35!|| message_mod ../starter/share/message_module/message_mod.F
36!|| submodel_mod ../starter/share/modules1/submodel_mod.F
37!||====================================================================
38 SUBROUTINE hm_read_mat04(
39 . LSUBMODEL ,MTAG ,UNITAB ,IPM ,PM ,
40 . MAT_ID ,TITR ,ISRATE ,MAT_PARAM)
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE elbuftag_mod
45 USE message_mod
46 USE submodel_mod
47 USE matparam_def_mod
48 USE unitab_mod
49 USE names_and_titles_mod , ONLY : nchartitle
50C-----------------------------------------------
51C ROUTINE DESCRIPTION :
52C ===================
53C READ MAT LAW04 WITH HM READER
54C-----------------------------------------------
55C DUMMY ARGUMENTS DESCRIPTION:
56C ===================
57C UNITAB UNITS ARRAY
58C MAT_ID MATERIAL ID(INTEGER)
59C TITR MATERIAL TITLE
60C LSUBMODEL SUBMODEL STRUCTURE
61C-----------------------------------------------
62C I m p l i c i t T y p e s
63C-----------------------------------------------
64#include "implicit_f.inc"
65C-----------------------------------------------
66C C o m m o n B l o c k s
67C-----------------------------------------------
68#include "scr03_c.inc"
69#include "units_c.inc"
70#include "param_c.inc"
71C-----------------------------------------------
72C D u m m y A r g u m e n t s
73C-----------------------------------------------
74 INTEGER, INTENT(INOUT) :: MAT_ID
75 INTEGER, INTENT(INOUT) :: ISRATE
76 INTEGER, DIMENSION(NPROPMI) ,INTENT(INOUT) :: IPM
77 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
78 my_real, DIMENSION(NPROPM) ,INTENT(INOUT) :: pm
79 TYPE(unit_type_) ,INTENT(IN) :: UNITAB
80 TYPE(submodel_data), DIMENSION(NSUBMOD),INTENT(IN) :: LSUBMODEL
81 TYPE(mlaw_tag_), INTENT(INOUT) :: MTAG
82 TYPE(matparam_struct_),INTENT(INOUT) :: MAT_PARAM
83C-----------------------------------------------
84C L o c a l V a r i a b l e s
85C-----------------------------------------------
86 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED,IS_EOS, EOS_EMBEDDED
87 INTEGER JTUR, I, EMPTY_LINE,COUNT,IEOS,LAW6_OPT
88 INTEGER JTHE, IFLAGR
89 my_real
90 . young, anu, ca, cb, cn, epsm, sigm, bulk,
91 . pmin, cc, eps0, cm, tmelt, tmax, cs, t0, sph, as, bs,
92 . al, bl, g, tba, tbk, pco, fac_l, fac_t, fac_m, fac_c,
93 . bid2 ,bid3 ,bid4 ,bid5 ,
94 . bid2_,bid3_,bid4_,bid5_,
95 . c0,c1,c2,c3,c4,c5,e0,psh,ti,rho0,rhor
96C-----------------------------------------------
97C S o u r c e L i n e s
98C-----------------------------------------------
99 count = 0
100 is_eos=.false.
101 ieos = 0
102 eos_embedded=.false.
103 bulk=zero
104 israte = 0
105 is_encrypted = .false.
106 is_available = .false.
107!---
108 CALL hm_option_is_encrypted(is_encrypted)
109
110 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
111 CALL hm_get_floatv('Refer_Rho' ,rhor ,is_available, lsubmodel, unitab)
112
113 CALL hm_get_floatv('MAT_E' ,young ,is_available, lsubmodel, unitab)
114 CALL hm_get_floatv('MAT_NU' ,anu ,is_available, lsubmodel, unitab)
115
116 CALL hm_get_floatv('MAT_SIGY' ,ca ,is_available, lsubmodel, unitab)
117 CALL hm_get_floatv('MAT_BETA' ,cb ,is_available, lsubmodel, unitab)
118 CALL hm_get_floatv('MAT_HARD' ,cn ,is_available, lsubmodel, unitab)
119 CALL hm_get_floatv('MAT_EPS' ,epsm ,is_available, lsubmodel, unitab)
120 CALL hm_get_floatv('MAT_SIG' ,sigm ,is_available, lsubmodel, unitab)
121
122 CALL hm_get_floatv('MAT_PC' ,pmin ,is_available, lsubmodel, unitab)
123
124 CALL hm_get_floatv('MAT_SRC' ,cc ,is_available, lsubmodel, unitab)
125 CALL hm_get_floatv('MAT_SRP' ,eps0 ,is_available, lsubmodel, unitab)
126 CALL hm_get_floatv('MAT_M' ,cm ,is_available, lsubmodel, unitab)
127 CALL hm_get_floatv('MAT_TMELT' ,tmelt ,is_available, lsubmodel, unitab)
128 CALL hm_get_floatv('MAT_TMAX' ,tmax ,is_available, lsubmodel, unitab)
129
130 CALL hm_get_floatv('MAT_SPHEAT' ,cs ,is_available, lsubmodel, unitab)
131 CALL hm_get_floatv('MAT_T0' ,t0 ,is_available, lsubmodel, unitab)
132
133 CALL hm_get_intv ('Line_count' ,count ,is_available, lsubmodel)
134
135 IF(invers_init>=2018)THEN
136 eos_embedded=.false.
137 ELSEIF(invers_init>=110)THEN
138 IF(count==3)eos_embedded=.true.
139 ELSE
140 eos_embedded=.true.
141 ENDIF
142
143 !----------------------------------------------------
144 !CHECK IF OLD CARD FORMAT 100 IS USED (OBSOLETE)
145 ! MAKES IT COMPATIBLE + WARNING
146 !----------------------------------------------------
147 IF(eos_embedded)THEN ! will use POLYNOMIAL EOS
148 IF(invers_init>=140) THEN
149 CALL ancmsg(msgid=1072,
150 . msgtype=msgerror,
151 . anmode=aninfo,
152 . i1 = mat_id,
153 . c1 = titr)
154 RETURN
155 ENDIF
156 CALL hm_get_floatv('MAT_C0' ,c0 ,is_available, lsubmodel, unitab)
157 CALL hm_get_floatv('MAT_C1' ,c1 ,is_available, lsubmodel, unitab)
158 CALL hm_get_floatv('MAT_C2' ,c2 ,is_available, lsubmodel, unitab)
159 CALL hm_get_floatv('MAT_C3' ,c3 ,is_available, lsubmodel, unitab)
160 CALL hm_get_floatv('MAT_C4' ,c4 ,is_available, lsubmodel, unitab)
161 CALL hm_get_floatv('MAT_C5' ,c5 ,is_available, lsubmodel, unitab)
162 CALL hm_get_floatv('MAT_E0' ,e0 ,is_available, lsubmodel, unitab)
163 CALL hm_get_floatv('MAT_PSH' ,psh ,is_available, lsubmodel, unitab)
164 bulk = c1
165 pm(23)= e0
166 pm(31)= c0-psh
167 pm(32)= c1
168 pm(33)= c2
169 pm(34)= c3
170 pm(35)= c4
171 pm(36)= c5
172 pm(88)= psh
173 ieos = 1 ! Polynomial EOS
174 ELSE
175 ieos = 18 ! Linear EOS is used by default
176 ENDIF
177!
178 CALL init_mat_keyword(mat_param,"ELASTO_PLASTIC")
179 CALL init_mat_keyword(mat_param,"incremental")
180 CALL INIT_MAT_KEYWORD(MAT_PARAM,"large_strain")
181 MAT_PARAM%IEOS = IEOS
182 IPM(4) = IEOS ! keep this temporarily for output and fluid sections
183
184 IF (PMIN==ZERO) PMIN=-EP20
185 IF (RHOR==ZERO) RHOR=RHO0
186 PM(1) = RHOR
187 PM(89)= RHO0
188
189
190.AND. IF(CS == ZEROTMELT == ZERO)THEN
191 CALL ANCMSG(MSGID=593,
192 . MSGTYPE=MSGWARNING,
193 . ANMODE=ANINFO_BLIND_1,
194 . I1=MAT_ID,
195 . C1=TITR)
196 CS=EP20
197 ELSEIF(CS == ZERO)THEN
198 CALL ANCMSG(MSGID=594,
199 . MSGTYPE=MSGERROR,
200 . ANMODE=ANINFO_BLIND_1,
201 . I1=MAT_ID,
202 . C1=TITR)
203 END IF
204 SPH=CS
205 IF(PMIN == ZERO) PMIN =-EP20
206.OR. IF(CN == ZEROCN == ONE) CN = ONEP0001
207 IF(EPSM == ZERO) EPSM = EP20
208 IF(SIGM == ZERO) SIGM = EP20
209 IF(CC == ZERO) EPS0 = ONE
210 IF(CM == ZERO) CM = ONE
211 IF(TMELT == ZERO)TMELT = EP20
212 IF(TMAX == ZERO) TMAX = EP20
213 IF(T0<=ZERO) T0 = THREE100
214 G=YOUNG/(TWO*(ONE + ANU))
215.NOT. IF(EOS_EMBEDDED)BULK=YOUNG/(THREE*(ONE - TWO*ANU))
216C
217 PM(20)=YOUNG
218 PM(21)=ANU
219 PM(22)=G
220 PM(32)=BULK
221 PM(37)=PMIN
222 PM(38)=CA
223 PM(39)=CB
224 PM(40)=CN
225 PM(41)=EPSM
226 PM(42)=SIGM
227 PM(43)=CC
228 PM(44)=EPS0
229 PM(45)=CM
230 PM(46)=TMELT
231 PM(47)=TMAX
232 PM(48)=CS
233 PM(69)=SPH
234 PM(79)=T0
235 PM(80)=TMELT
236CC-----------
237C Formulation for solid elements time step computation.
238 IPM(252)= 2
239 PM(105) = (ONE -TWO*ANU)/(ONE - ANU) ! TWO*G/(BULK + FOUR_OVER_3*G)
240
241 !sizes for element buffer
242 MTAG%G_PLA = 1
243 MTAG%G_TEMP = 1
244 MTAG%G_EPSD = 1
245 !
246 MTAG%L_PLA = 1
247 MTAG%L_TEMP = 1
248 MTAG%L_EPSD = 1
249
250 ! Material compatibility with /EOS option
251 CALL INIT_MAT_KEYWORD(MAT_PARAM,"eos")
252
253 ! EOS/Thermo keyword for pressure treatment in elements
254 CALL INIT_MAT_KEYWORD(MAT_PARAM,"hydro_eos")
255
256 ! Properties compatibility
257 CALL INIT_MAT_KEYWORD(MAT_PARAM,"solid_isotropic")
258 CALL INIT_MAT_KEYWORD(MAT_PARAM,"sph")
259!
260 ! activate heat source calculation in material for /heat/mat
261
262 MAT_PARAM%HEAT_FLAG = 1
263
264C--------------------------------
265 WRITE(IOUT,2001) TITR,MAT_ID,4
266 WRITE(IOUT,1000)
267
268 IF(IS_ENCRYPTED)THEN
269 WRITE(IOUT,'(5X,A,//)')'CONFIDENTIAL DATA'
270 ELSE
271 WRITE(IOUT,2002)RHO0,RHOR
272 WRITE(IOUT,1300)YOUNG,ANU,G,BULK
273 WRITE(IOUT,1400)CA,CB,CN,EPSM,SIGM
274 WRITE(IOUT,1500)PMIN
275 WRITE(IOUT,1600)CC,EPS0,CM,TMELT,TMAX,CS,T0
276 ENDIF
277
278 !--------------------------------------------------!
279 ! OLD FORMAT 100 (OBSOLETE) !
280 !--------------------------------------------------!
281 IF(EOS_EMBEDDED)THEN !
282 WRITE(IOUT,2000) !
283 IF(IS_ENCRYPTED)THEN !
284 WRITE(IOUT,'(5X,A,//)')'CONFIDENTIAL DATA' !
285 ELSE !
286 WRITE(IOUT,2500)C0,C1,C2,C3,C4,C5,E0,PM(1),PSH !
287 ENDIF !
288 ENDIF !
289 !--------------------------------------------------!
290C
291 IF(EPS0 == ZERO) THEN
292 CALL ANCMSG(MSGID=298,MSGTYPE=MSGERROR,ANMODE=ANINFO,I1=4,I2=MAT_ID,C1=TITR)
293 ENDIF
294 RETURN
295 1000 FORMAT(
296 & 5X,' JOHNSON COOK LAW ',/,
297 & 5X,' ---------------- ',//)
298 1300 FORMAT(
299 & 5X,'YOUNG',1H','S MODULUS . . . . . . . . . . . .=',1PG20.13/,
300 & 5X,'POISSON',1H','S RATIO . . . . . . . . . . . .=',1PG20.13/,
301 & 5X,'SHEAR MODULUS . . . . . . . . . . . . .=',1PG20.13/,
302 & 5X,'BULK MODULUS . . . . . . . . . . . . .=',1PG20.13//)
303 1400 FORMAT(
304 & 5X,'YIELD COEFFICIENT CA. . . . . . . . . .=',1PG20.13/,
305 & 5X,'YIELD COEFFICIENT CB. . . . . . . . . .=',1PG20.13/,
306 & 5X,'YIELD COEFFICIENT CN. . . . . . . . . .=',1PG20.13/,
307 & 5X,'EPS-MAX . . . . . . . . . . . . . . . .=',1PG20.13/,
308 & 5X,'SIG-MAX . . . . . . . . . . . . . . . .=',1PG20.13//)
309 1500 FORMAT(
310 & 5X,'PRESSURE CUTOFF . . . . . . . . . . . .=',1PG20.13/)
311 1600 FORMAT(
312 & 5X,'STRAIN RATE COEFFICIENT CC. . . . . . .=',1PG20.13/,
313 & 5X,'REFERENCE STRAIN RATE . . . . . . . . .=',1PG20.13/,
314 & 5X,'TEMPERATURE EXPONENT. . . . . . . . . .=',1PG20.13/,
315 & 5X,'MELTING TEMPERATURE DEGREE K. . . . . .=',1PG20.13/,
316 & 5X,'THETA-MAX . . . . . . . . . . . . . . .=',1PG20.13/,
317 & 5X,'SPECIFIC HEAT . . . . . . . . . . . . .=',1PG20.13/,
318 & 5X,'ROOM TEMPERATURE . . . . . . . . . . . =',1PG20.13//)
319 2000 FORMAT(
320 & 5X,' POLYNOMIAL EOS ',/,
321 & 5X,' -------------- ',/,
322 & 5X,' INPUT 100 IS OBSOLETE ',/,
323 & 5X,' IT CAN BE ERASED BY /EOS CARD ',/)
324 2001 FORMAT(/
325 & 5X,A,/,
326 & 5X,' MATERIAL NUMBER . . . . . . . . . . . .=',I10/,
327 & 5X,' MATERIAL LAW. . . . . . . . . . . . . .=',I10/)
328 2002 FORMAT(
329 & 5X,' INITIAL DENSITY . . . . . . . . . . . .=',1PG20.13/,
330 & 5X,' REFERENCE DENSITY . . . . . . . . . . .=',1PG20.13/)
331 2500 FORMAT(
332 & 5X,'C0. . . . . . . . . . . . . . . . . . . .=',1PG20.13/,
333 & 5X,'C1. . . . . . . . . . . . . . . . . . . .=',1PG20.13/,
334 & 5X,'C2. . . . . . . . . . . . . . . . . . . .=',1PG20.13/,
335 & 5X,'C3. . . . . . . . . . . . . . . . . . . .=',1PG20.13/,
336 & 5X,'C4. . . . . . . . . . . . . . . . . . . .=',1PG20.13/,
337 & 5X,'C5. . . . . . . . . . . . . . . . . . . .=',1PG20.13/,
338 & 5X,'INITIAL INTERNAL ENERGY PER UNIT VOLUME .=',1PG20.13/,
339 & 5X,'REFERENCE DENSITY . . . . . . . . . . . .=',1PG20.13/,
340 & 5X,'PRESSURE SHIFT. . . . . . . . . . . . . .=',1PG20.13//)
341 RETURN
342 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_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_mat04
subroutine hm_read_mat04(lsubmodel, mtag, unitab, ipm, pm, mat_id, titr, israte, mat_param)
Definition hm_read_mat04.F:41
init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:150
elbuftag_mod
Definition elbuftag_mod.F:227
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
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
starter
program starter
Definition starter.F:39
elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:232
submodel_mod::submodel_data
Definition submodel_mod.F:723
unitab_mod::unit_type_
Definition unitab_mod.F:520