OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat03.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_mat03 ../starter/source/materials/mat/mat003/hm_read_mat03.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_mat03(
39 . LSUBMODEL ,MTAG ,UNITAB ,IPM ,PM ,
40 . MAT_ID ,TITR ,ISRATE ,MATPARAM)
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE elbuftag_mod !MAXEOS
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 LAW03 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) :: MATPARAM
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 COUNT, IEOS
88 my_real
89 . young, anu, ca, cb, cn, epsm, sigm, bulk,
90 . pmin, e0, g,
91 .
92 .
93 . c0,c1,c2,c3,c4,c5,psh,rho0,rhor
94C-----------------------------------------------
95C S o u r c e L i n e s
96C-----------------------------------------------
97 count = 0
98 is_eos=.false.
99 ieos=0
100 eos_embedded=.false.
101 bulk=zero
102 is_encrypted = .false.
103 is_available = .false.
104!---
105 CALL hm_option_is_encrypted(is_encrypted)
106
107 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
108 CALL hm_get_floatv('Refer_Rho' ,rhor ,is_available, lsubmodel, unitab)
109
110 CALL hm_get_floatv('MAT_E' ,young ,is_available, lsubmodel, unitab)
111 CALL hm_get_floatv('MAT_NU' ,anu ,is_available, lsubmodel, unitab)
112
113 CALL hm_get_floatv('MAT_SIGY' ,ca ,is_available, lsubmodel, unitab)
114 CALL hm_get_floatv('MAT_BETA' ,cb ,is_available, lsubmodel, unitab)
115 CALL hm_get_floatv('MAT_HARD' ,cn ,is_available, lsubmodel, unitab)
116 CALL hm_get_floatv('MAT_EPS' ,epsm ,is_available, lsubmodel, unitab)
117 CALL hm_get_floatv('MAT_SIG' ,sigm ,is_available, lsubmodel, unitab)
118
119 CALL hm_get_floatv('MAT_PC' ,pmin ,is_available, lsubmodel, unitab)
120
121 CALL hm_get_intv ('Line_count' ,count ,is_available, lsubmodel)
122
123 IF(invers_init>=2018)THEN
124 eos_embedded=.false.
125 ELSEIF(invers_init>=110)THEN
126 IF(count==3)eos_embedded=.true.
127 ELSE
128 eos_embedded=.true.
129 ENDIF
130
131 !----------------------------------------------------
132 !CHECK IF OLD CARD FORMAT 100 IS USED (OBSOLETE)
133 ! MAKES IT COMPATIBLE + WARNING
134 !----------------------------------------------------
135 IF(eos_embedded)THEN ! POLYNOMIAL EOS
136 IF(invers_init>=140) THEN
137 CALL ancmsg(msgid=1072,
138 . msgtype=msgerror,
139 . anmode=aninfo,
140 . i1 = mat_id,
141 . c1 = titr)
142 RETURN
143 ENDIF
144 CALL hm_get_floatv('MAT_C0' ,c0 ,is_available, lsubmodel, unitab)
145 CALL hm_get_floatv('MAT_C1' ,c1 ,is_available, lsubmodel, unitab)
146 CALL hm_get_floatv('MAT_C2' ,c2 ,is_available, lsubmodel, unitab)
147 CALL hm_get_floatv('MAT_C3' ,c3 ,is_available, lsubmodel, unitab)
148 CALL hm_get_floatv('MAT_C4' ,c4 ,is_available, lsubmodel, unitab)
149 CALL hm_get_floatv('MAT_C5' ,c5 ,is_available, lsubmodel, unitab)
150 CALL hm_get_floatv('MAT_E0' ,e0 ,is_available, lsubmodel, unitab)
151 CALL hm_get_floatv('MAT_PSH' ,psh ,is_available, lsubmodel, unitab)
152 ! allocate EoS data structure
153 matparam%IEOS = 1
154 matparam%EOS%NUPARAM = 7
155 matparam%EOS%NIPARAM = 0
156 matparam%EOS%NFUNC = 0
157 matparam%EOS%NTABLE = 0
158 CALL matparam%EOS%CONSTRUCT() !allocations
159 matparam%EOS%UPARAM(1) = c0-psh
160 matparam%EOS%UPARAM(2) = c1
161 matparam%EOS%UPARAM(3) = c2
162 matparam%EOS%UPARAM(4) = c3
163 matparam%EOS%UPARAM(5) = c4
164 matparam%EOS%UPARAM(6) = c5
165 matparam%EOS%UPARAM(7) = zero
166 matparam%EOS%PSH = psh
167 matparam%EOS%E0 = e0
168 bulk = c1
169 pm(23)= e0
170 pm(31)= c0-psh
171 pm(32) = c1 !legacy value
172 !!PM(32)= C1+C4*ABS(E0) ! bulk modulus
173 pm(88)= psh
174 ieos = 1 ! POLYNOMIAL EOS
175 ELSE
176 ieos = 18 ! linear EOS is used by default
177 ENDIF
178!
179 matparam%IEOS = ieos
180 ipm(4) = ieos ! keep this temporarily for output and fluid sections
181
182 IF (pmin==zero) pmin=-ep20
183 IF (rhor==zero) rhor=rho0
184 pm(1) = rhor
185 pm(89)= rho0
186 IF(cn==zero.OR.cn==one) cn = onep0001
187 IF(epsm==zero) epsm = ep20
188 IF(sigm==zero) sigm = ep20
189 g = young/(two*(one +anu))
190 IF(.NOT.eos_embedded)bulk = young/(three*(one - two*anu))
191 IF(.NOT.eos_embedded)e0 = zero
192
193 pm(20)=young
194 pm(21)=anu
195 pm(22)=g
196 pm(23)=e0
197 pm(32)=bulk
198 pm(37)=pmin
199 pm(38)=ca
200 pm(39)=cb
201 pm(40)=cn
202 pm(41)=epsm
203 pm(42)=sigm
204 pm(79)=three100
205 !Formulation for solid elements time step computation.
206 ipm(252)= 2
207 pm(105) = (one -two*anu)/(one - anu) ! TWO*G/(BULK + FOUR_OVER_3*G)
208
209 mtag%G_PLA = 1
210 mtag%L_PLA = 1
211 mtag%G_TEMP = 1
212 mtag%L_TEMP = 1
213
214 CALL init_mat_keyword(matparam,"ELASTO_PLASTIC")
215 CALL init_mat_keyword(matparam,"INCREMENTAL")
216 CALL init_mat_keyword(matparam,"LARGE_STRAIN")
217 ! Material compatibility with /EOS
218 CALL init_mat_keyword(matparam,"EOS")
219
220 ! EOS/Thermo keyword
221 CALL init_mat_keyword(matparam,"HYDRO_EOS")
222
223 ! Properties compatibility
224 CALL init_mat_keyword(matparam,"SOLID_ISOTROPIC")
225 CALL init_mat_keyword(matparam,"SPH")
226
227 WRITE(iout,2001) titr,mat_id,3
228 WRITE(iout,1000)
229
230 IF(is_encrypted)THEN
231 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
232 ELSE
233 WRITE(iout,2002)rho0,rhor
234 WRITE(iout,1300)young,anu,g,bulk
235 WRITE(iout,1400)ca,cb,cn,epsm,sigm
236 WRITE(iout,1500)pmin
237 ENDIF
238 !--------------------------------------------------!
239 ! OLD FORMAT 100 (OBSOLETE) !
240 !--------------------------------------------------!
241 IF(eos_embedded)THEN !
242 WRITE(iout,2000) !
243 IF(is_encrypted)THEN !
244 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA' !
245 ELSE !
246 WRITE(iout,2500)c0,c1,c2,c3,c4,c5,e0,pm(1),psh !
247 ENDIF !
248 ENDIF !
249 !--------------------------------------------------!
250C-----------
251 RETURN
252 1000 FORMAT(
253 & 5x,' ELASTIC PLASTIC HYDRODYNAMIC ',/,
254 & 5x,' ---------------------------- ',//)
255 1300 FORMAT(
256 & 5x,'YOUNG',1h','s modulus . . . . . . . . . . . .=',1PG20.13/,
257 & 5X,'poisson',1H','S RATIO . . . . . . . . . . . .=',1pg20.13/,
258 & 5x,'SHEAR MODULUS . . . . . . . . . . . . . .=',1pg20.13/,
259 & 5x,'BULK MODULUS . . . . . . . . . . . . . .=',1pg20.13//)
260 1400 FORMAT(
261 & 5x,'YIELD COEFFICIENT CA. . . . . . . . . .=',1pg20.13/,
262 & 5x,'YIELD COEFFICIENT CB. . . . . . . . . .=',1pg20.13/,
263 & 5x,'YIELD COEFFICIENT CN. . . . . . . . . .=',1pg20.13/,
264 & 5x,'EPS-MAX . . . . . . . . . . . . . . . .=',1pg20.13/,
265 & 5x,'SIG-MAX . . . . . . . . . . . . . . . .=',1pg20.13//)
266 1500 FORMAT(
267 & 5x,'PRESSURE CUTOFF . . . . . . . . . . . .=',1pg20.13//)
268 2000 FORMAT(
269 & 5x,' POLYNOMIAL EOS ',/,
270 & 5x,' -------------- ',/,
271 & 5x,' INPUT 100 IS OBSOLETE ',/,
272 & 5x,' IT CAN BE ERASED BY /EOS CARD ',/)
273 2500 FORMAT(
274 & 5x,'C0. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
275 & 5x,'C1. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
276 & 5x,'C2. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
277 & 5x,'C3. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
278 & 5x,'C4. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
279 & 5x,'C5. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
280 & 5x,'INITIAL INTERNAL ENERGY PER UNIT VOLUME .=',1pg20.13/,
281 & 5x,'REFERENCE DENSITY . . . . . . . . . . . .=',1pg20.13/,
282 & 5x,'PRESSURE SHIFT. . . . . . . . . . . . . .=',1pg20.13//)
283 2001 FORMAT(/
284 & 5x,a,/,
285 & 5x,' MATERIAL NUMBER . . . . . . . . . . . .=',i10/,
286 & 5x,' MATERIAL LAW. . . . . . . . . . . . . .=',i10/)
287 2002 FORMAT(
288 & 5x,' INITIAL DENSITY . . . . . . . . . . . .=',1pg20.13/,
289 & 5x,' REFERENCE DENSITY . . . . . . . . . . .=',1pg20.13/)
290
291 RETURN
292 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: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_mat03
subroutine hm_read_mat03(lsubmodel, mtag, unitab, ipm, pm, mat_id, titr, israte, matparam)
Definition hm_read_mat03.F:41
init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:153
elbuftag_mod
Definition elbuftag_mod.F:230
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
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
starter
program starter
Definition starter.F:39
elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:235
submodel_mod::submodel_data
Definition submodel_mod.F:728
unitab_mod::unit_type_
Definition unitab_mod.F:524