OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_eos.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_eos ../starter/source/materials/eos/hm_read_eos.f
25!||--- called by ------------------------------------------------------
26!|| read_material_models ../starter/source/materials/read_material_models.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| hm_get_string ../starter/source/devtools/hm_reader/hm_get_string.F
30!|| hm_option_count ../starter/source/devtools/hm_reader/hm_option_count.F
31!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
32!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
33!|| hm_read_eos_compaction ../starter/source/materials/eos/hm_read_eos_compaction.F90
34!|| hm_read_eos_compaction2 ../starter/source/materials/eos/hm_read_eos_compaction2.F90
35!|| hm_read_eos_compaction_tab ../starter/source/materials/eos/hm_read_eos_compaction_tab.F90
36!|| hm_read_eos_exponential ../starter/source/materials/eos/hm_read_eos_exponential.F90
37!|| hm_read_eos_gruneisen ../starter/source/materials/eos/hm_read_eos_gruneisen.F
38!|| hm_read_eos_ideal_gas ../starter/source/materials/eos/hm_read_eos_ideal_gas.F
39!|| hm_read_eos_ideal_gas_vt ../starter/source/materials/eos/hm_read_eos_ideal_gas_vt.F
40!|| hm_read_eos_linear ../starter/source/materials/eos/hm_read_eos_linear.F
41!|| hm_read_eos_lszk ../starter/source/materials/eos/hm_read_eos_lszk.F
42!|| hm_read_eos_murnaghan ../starter/source/materials/eos/hm_read_eos_murnaghan.F
43!|| hm_read_eos_nasg ../starter/source/materials/eos/hm_read_eos_nasg.F
44!|| hm_read_eos_noble_abel ../starter/source/materials/eos/hm_read_eos_noble_abel.F
45!|| hm_read_eos_osborne ../starter/source/materials/eos/hm_read_eos_osborne.F
46!|| hm_read_eos_polynomial ../starter/source/materials/eos/hm_read_eos_polynomial.F
47!|| hm_read_eos_powderburn ../starter/source/materials/eos/hm_read_eos_powderburn.f90
48!|| hm_read_eos_puff ../starter/source/materials/eos/hm_read_eos_puff.F
49!|| hm_read_eos_sesame ../starter/source/materials/eos/hm_read_eos_sesame.F
50!|| hm_read_eos_stiffened_gas ../starter/source/materials/eos/hm_read_eos_stiffened_gas.F
51!|| hm_read_eos_tabulated ../starter/source/materials/eos/hm_read_eos_tabulated.F
52!|| hm_read_eos_tillotson ../starter/source/materials/eos/hm_read_eos_tillotson.F
53!||--- uses -----------------------------------------------------
54!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
55!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
56!|| hm_read_eos_compaction2_mod ../starter/source/materials/eos/hm_read_eos_compaction2.F90
57!|| hm_read_eos_compaction_mod ../starter/source/materials/eos/hm_read_eos_compaction.F90
58!|| hm_read_eos_compaction_tab_mod ../starter/source/materials/eos/hm_read_eos_compaction_tab.F90
59!|| hm_read_eos_exponential_mod ../starter/source/materials/eos/hm_read_eos_exponential.F90
60!|| message_mod ../starter/share/message_module/message_mod.F
61!|| submodel_mod ../starter/share/modules1/submodel_mod.F
62!|| table_mod ../starter/share/modules1/table_mod.F
63!||====================================================================
64 SUBROUTINE hm_read_eos(
65 . MAT_PARAM ,IPM ,PM ,BUFMAT ,NUMMAT ,NFUNCT ,
66 . BUFLEN ,IADBUF ,EOS_TAG ,UNITAB ,LSUBMODEL,
67 . MLAW_TAG ,NPC ,TF ,SNPC ,NPTS ,SBUFMAT,
68 . NTABLE ,TABLE)
69C-----------------------------------------------
70C A n a l y s e M o d u l e
71C-----------------------------------------------
72 USE matparam_def_mod
73 USE unitab_mod
74 USE message_mod
75 USE elbuftag_mod
76 USE submodel_mod , ONLY : submodel_data, nsubmod
77 USE hm_option_read_mod
78 USE names_and_titles_mod , ONLY : nchartitle, ncharkey
79 USE hm_read_eos_exponential_mod , ONLY : hm_read_eos_exponential
80 USE hm_read_eos_compaction_mod , ONLY : hm_read_eos_compaction
81 USE hm_read_eos_compaction2_mod , ONLY : hm_read_eos_compaction2
82 USE hm_read_eos_compaction_tab_mod , ONLY : hm_read_eos_compaction_tab
83 USE table_mod , ONLY : ttable
84 !USE EOS_PARAM_MOD , ONLY : EOS_PARAM_
85C-----------------------------------------------
86C I m p l i c i t T y p e s
87C-----------------------------------------------
88#include "implicit_f.inc"
89C-----------------------------------------------
90C C o m m o n B l o c k s
91C-----------------------------------------------
92#include "units_c.inc"
93#include "param_c.inc"
94C-----------------------------------------------
95C D u m m y A r g u m e n t s
96C-----------------------------------------------
97 INTEGER,INTENT(IN) :: NUMMAT, NTABLE, NFUNCT
98 TYPE(MLAW_TAG_) , DIMENSION(NUMMAT),INTENT(INOUT) :: MLAW_TAG
99 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
100 INTEGER :: BUFLEN,IADBUF
101 INTEGER,INTENT(IN) :: SBUFMAT
102 INTEGER ,DIMENSION(NPROPMI,NUMMAT) ,INTENT(INOUT) :: IPM
103 INTEGER,INTENT(IN) :: SNPC, NPTS
104 INTEGER,INTENT(IN) :: NPC(SNPC)
105 my_real,INTENT(IN) :: tf(npts)
106 my_real ,DIMENSION(NPROPM ,NUMMAT) ,INTENT(INOUT) :: pm
107 my_real ,DIMENSION(SBUFMAT),INTENT(INOUT) :: bufmat
108 TYPE(matparam_struct_) ,DIMENSION(NUMMAT) ,TARGET, INTENT(INOUT) :: MAT_PARAM
109 TYPE(eos_tag_),DIMENSION(0:MAXEOS) ,INTENT(INOUT) :: EOS_TAG
110 TYPE(submodel_data) ,DIMENSION(NSUBMOD) ,INTENT(IN) :: LSUBMODEL
111 TYPE(ttable) ,DIMENSION(NTABLE) ,INTENT(IN) :: TABLE
112C-----------------------------------------------
113C L o c a l V a r i a b l e s
114C-----------------------------------------------
115 INTEGER :: I,IMAT,IEOS,HM_NEOS,IUNIT,FLAGUNIT,FLAGMAT,IMID,EOS_uid,MTN,UNIT_ID
116 CHARACTER(LEN=NCHARTITLE) :: TITR
117 CHARACTER KEY2*32, EOS_key*32
118 LOGICAL IS_AVAILABLE
119 INTEGER ISFLUID
120 TYPE(eos_param_),POINTER :: EOS_STRUCT
121c--------------------------------------------------
122c Routine reading EOS models :
123 !------------------------------------!
124 ! IEOS ! EOS !
125 !-----------!------------------------!
126 ! 0 => 18 (LINEAR) !
127 ! 1 ! POLYNOMIAL !
128 ! 2 ! GRUNEISEN !
129 ! 3 ! TILLOTSON !
130 ! 4 ! PUFF !
131 ! 5 ! SESAME !
132 ! 6 ! NOBLE-ABEL ! 2017.0
133 ! 7 ! IDEAL GAS ! 2018.0
134 ! 8 ! MUNAGHAN ! 2018.0
135 ! 9 ! OSBORNE ! 2018.0
136 ! 10 ! STIFFENED GAS ! 2018.0
137 ! 11 ! LSZK ! 2018.0
138 ! 12 ! POWDER-BURN ! 2019.1
139 ! 13 ! COMPACTION ! 2019.1
140 ! 14 ! NASG ! 2020.0
141 ! 15 ! JWL ! internal use : /INIMAP
142 ! 16 ! IDEALGAS_VT ! 2022.0
143 ! 17 ! TABULATED ! 2022.2
144 ! 18 ! LINEAR ! 2019.0
145 ! 19 ! EXPONENTIAL ! 2024.1
146 ! 20 ! COMPACTION2 ! 2025.1
147 ! 21 ! COMPACTION_TAB ! 2026.0
148 !------------------------------------!
149c======================================================================-
150c COUNT EOS MODELS USING CFG FILES
151c--------------------------------------------------
152c
153 CALL hm_option_count('/EOS',hm_neos)
154c
155c--------------------------------------------------
156c START BROWSING EOS MODELS
157c--------------------------------------------------
158c
159 CALL hm_option_start('/EOS')
160c--------------------------------------------------
161c
162 DO i = 1,hm_neos
163 CALL hm_option_read_key(lsubmodel, option_id=eos_uid, option_titr=titr, unit_id=unit_id)
164 CALL hm_get_string('EOS_Keyword',key2, 32,is_available)
165 eos_key = key2(1:len_trim(key2))
166 ! Check MAT_Id
167 flagmat = 0
168 DO imat=1,nummat-1
169 imid = ipm(1,imat)
170 IF (imid == eos_uid) THEN
171 flagmat = 1
172 EXIT
173 ENDIF
174 ENDDO
175 IF (eos_uid > 0 .AND. flagmat == 0) THEN
176 CALL ancmsg(msgid=1663, anmode=aninfo, msgtype=msgerror, i1=eos_uid, c1='EOS', c2=eos_key, c3=titr)
177 ENDIF
178
179 ! Check Unit_ID
180 flagunit = 0
181 DO iunit=1,unitab%NUNITS
182 IF (unitab%UNIT_ID(iunit) == unit_id) THEN
183 flagunit = 1
184 EXIT
185 ENDIF
186 ENDDO
187 IF (unit_id > 0 .AND. flagunit == 0) THEN
188 CALL ancmsg(msgid=659,anmode=aninfo,msgtype=msgerror,i1=eos_uid,i2=unit_id,c1='MATERIAL',c2='EOS MODEL',c3=titr)
189 ENDIF
190
191 !initialize %EOS member for current material
192 mat_param(imat)%EOS%TITLE = ''
193 mat_param(imat)%EOS%TITLE(1:len_trim(titr)) = titr(1:len_trim(titr))
194 mat_param(imat)%EOS%NUPARAM = 0
195 mat_param(imat)%EOS%NIPARAM = 0
196 mat_param(imat)%EOS%NFUNC = 0
197 mat_param(imat)%EOS%NTABLE = 0
198 mat_param(imat)%EOS%ISFLUID = 0
199 eos_struct => mat_param(imat)%EOS
200 isfluid = 0
201
202 WRITE(iout,1000) trim(titr),eos_uid
203 eos_struct%TITLE = eos_key
204 ! READ EOS
205 SELECT CASE (eos_key)
206c---
207 CASE ('POLYNOMIAL')
208 ieos = 1
209 CALL hm_read_eos_polynomial(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,isfluid,mat_param(imat))
210c---
211 CASE ('GRUNEISEN')
212 ieos = 2
213 CALL hm_read_eos_gruneisen(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
214c---
215 CASE ('TILLOTSON')
216 ieos = 3
217 CALL hm_read_eos_tillotson(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,eos_tag,ieos,mat_param(imat))
218c---
219 CASE ('PUFF')
220 ieos = 4
221 CALL hm_read_eos_puff(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
222c---
223 CASE ('SESAME')
224 ieos = 5
225 CALL hm_read_eos_sesame(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,imid,titr,bufmat,buflen,iadbuf,mat_param(imat))
226c---
227 CASE ('NOBLE-ABEL','NA')
228 ieos = 6
229 CALL hm_read_eos_noble_abel(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
230 isfluid = 1
231c---
232 CASE ('IDEAL-GAS','IDEAL-GAS-VE')
233 ieos = 7
234 CALL hm_read_eos_ideal_gas(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid, mat_param(imat))
235 isfluid = 1
236c---
237 CASE ('MURNAGHAN')
238 ieos = 8
239 CALL hm_read_eos_murnaghan(npropm, iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid, mat_param(imat))
240c---
241 CASE ('OSBORNE')
242 ieos = 9
243 CALL hm_read_eos_osborne(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
244c---
245 CASE ('STIFF-GAS','SG')
246 ieos = 10
247 CALL hm_read_eos_stiffened_gas(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
248 isfluid = 1
249c---
250 CASE ('LSZK')
251 ieos = 11
252 CALL hm_read_eos_lszk(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
253 isfluid = 1
254c---
255 CASE ('POWDER-BURN')
256 ieos = 12
257 CALL hm_read_eos_powderburn(iout,pm(1,imat),ipm(1,imat),unitab,lsubmodel,eos_uid, mat_param(imat)%EOS,
258 . npropm, npropmi, mlaw_tag(imat), eos_tag, ieos )
259c---
260 CASE ('COMPACTION')
261 ieos = 13
262 CALL hm_read_eos_compaction(iout,pm(1,imat),unitab,lsubmodel,eos_uid,eos_tag,ieos,npropm,maxeos,eos_struct)
263c---
264 CASE ('NASG')
265 ieos = 14
266 CALL hm_read_eos_nasg(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
267 isfluid = 1
268c---
269 CASE ('IDEAL-GAS-VT')
270 ieos = 16
271 CALL hm_read_eos_ideal_gas_vt(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
272 isfluid = 1
273c---
274 CASE ('TABULATED')
275 ieos = 17
276 CALL hm_read_eos_tabulated(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid, eos_struct)
277c---
278 CASE ('LINEAR')
279 ieos = 18
280 CALL hm_read_eos_linear(iout,pm(1,imat),unitab,iunit,lsubmodel,eos_uid,mat_param(imat))
281c---
282 CASE ('EXPONENTIAL')
283 ieos = 19
284 CALL hm_read_eos_exponential(iout,pm(1,imat),unitab,lsubmodel,npropm, eos_struct)
285c---
286 CASE ('COMPACTION2')
287 ieos = 20
288 CALL hm_read_eos_compaction2(iout,pm(1,imat),unitab,lsubmodel,eos_uid,eos_tag,ieos,npropm,maxeos,
289 . eos_struct, nfunct, npc, tf ,snpc , npts, mat_param(imat)%REZON)
290 CASE ('COMPACTION_TAB')
291 ieos = 21
292 CALL hm_read_eos_compaction_tab(iout,pm(1,imat),unitab,lsubmodel,eos_uid,eos_tag,ieos,npropm,maxeos,
293 . eos_struct, ntable, table, mat_param(imat)%REZON)
294c---
295 CASE DEFAULT
296 ieos = -1
297 CALL ancmsg(msgid=67,anmode=aninfo,msgtype=msgerror,i1= eos_uid,c1=key2,c2='EOS KEYWORD IS NOT VALID')
298c---
299 END SELECT
300
301 ipm(4,imat) = ieos
302 mat_param(imat)%IEOS = ieos
303 mat_param(imat)%EOS%EOSTYPE = ieos
304 mat_param(imat)%EOS%ISFLUID = isfluid
305c--------------------------------------------------
306 ENDDO ! HM_NEOS
307c--------------------------------------------------
308 DO i=1,nummat-1
309 mtn = ipm(2,i)
310 IF (mtn /= 42) pm(100,i) = pm(32,i)
311 ENDDO
312c-----------
313 RETURN
314c--------------------------------------------------
315 1000 FORMAT(//
316 & 5x,'EQUATION OF STATE ',/,
317 & 5x,a,/,
318 & 5x,'MATERIAL ID . . . . . . . . . . . .=',i10/)
319c--------------------------------------------------
320 END
my_real
#define my_real
Definition cppsort.cpp:32
hm_get_string
subroutine hm_get_string(name, sval, size, is_available)
Definition hm_get_string.F:78
hm_option_count
subroutine hm_option_count(entity_type, hm_option_number)
Definition hm_option_count.F:105
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:253
hm_read_eos
subroutine hm_read_eos(mat_param, ipm, pm, bufmat, nummat, nfunct, buflen, iadbuf, eos_tag, unitab, lsubmodel, mlaw_tag, npc, tf, snpc, npts, sbufmat, ntable, table)
Definition hm_read_eos.F:69
hm_read_eos_gruneisen
subroutine hm_read_eos_gruneisen(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_gruneisen.F:36
hm_read_eos_ideal_gas
subroutine hm_read_eos_ideal_gas(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_ideal_gas.F:36
hm_read_eos_ideal_gas_vt
subroutine hm_read_eos_ideal_gas_vt(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_ideal_gas_vt.F:36
hm_read_eos_linear
subroutine hm_read_eos_linear(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_linear.F:35
hm_read_eos_lszk
subroutine hm_read_eos_lszk(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_lszk.F:36
hm_read_eos_murnaghan
subroutine hm_read_eos_murnaghan(npropm, iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_murnaghan.F:36
hm_read_eos_nasg
subroutine hm_read_eos_nasg(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_nasg.F:36
hm_read_eos_noble_abel
subroutine hm_read_eos_noble_abel(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_noble_abel.F:35
hm_read_eos_osborne
subroutine hm_read_eos_osborne(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_osborne.F:36
hm_read_eos_polynomial
subroutine hm_read_eos_polynomial(iout, pm, unitab, iunit, lsubmodel, imideos, isfluid, mat_param)
Definition hm_read_eos_polynomial.F:35
hm_read_eos_puff
subroutine hm_read_eos_puff(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_puff.F:35
hm_read_eos_sesame
subroutine hm_read_eos_sesame(iout, pm, unitab, iunit, lsubmodel, imideos, imid, titr, bufmat, mfi, idf, mat_param)
Definition hm_read_eos_sesame.F:41
hm_read_eos_stiffened_gas
subroutine hm_read_eos_stiffened_gas(iout, pm, unitab, iunit, lsubmodel, imideos, mat_param)
Definition hm_read_eos_stiffened_gas.F:36
hm_read_eos_tabulated
subroutine hm_read_eos_tabulated(iout, pm, unitab, iunit, lsubmodel, imideos, eos_struct)
Definition hm_read_eos_tabulated.F:37
hm_read_eos_tillotson
subroutine hm_read_eos_tillotson(iout, pm, unitab, iunit, lsubmodel, imideos, eos_tag, ieos, mat_param)
Definition hm_read_eos_tillotson.F:36
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:414
elbuftag_mod
Definition elbuftag_mod.F:230
elbuftag_mod::maxeos
integer, save maxeos
Definition elbuftag_mod.F:372
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
names_and_titles_mod::ncharkey
integer, parameter ncharkey
Definition names_and_titles_mod.F:1005
submodel_mod
Definition submodel_mod.F:717
submodel_mod::nsubmod
integer nsubmod
Definition submodel_mod.F:725
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
starter
program starter
Definition starter.F:39
elbuftag_mod::eos_tag_
Definition elbuftag_mod.F:356
submodel_mod::submodel_data
Definition submodel_mod.F:728
table_mod::ttable
Definition table_mod.F:126