OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat11.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_mat11 ../starter/source/materials/mat/mat011/hm_read_mat11.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_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!|| init_mat_keyword ../starter/source/materials/mat/init_mat_keyword.F
34!||--- uses -----------------------------------------------------
35!|| elbuftag_mod ../starter/share/modules1/elbuftag_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_mat11(LSUBMODEL, MTAG, UNITAB, IPM, PM,
40 . MAT_ID, TITR ,MATPARAM)
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE unitab_mod
45 USE elbuftag_mod
46 USE message_mod
47 USE submodel_mod
48 USE ale_mod
49 USE matparam_def_mod
50 USE names_and_titles_mod , ONLY : nchartitle
51C-----------------------------------------------
52C I m p l i c i t T y p e s
53C-----------------------------------------------
54#include "implicit_f.inc"
55C-----------------------------------------------
56C C o m m o n B l o c k s
57C-----------------------------------------------
58#include "units_c.inc"
59#include "param_c.inc"
60C-----------------------------------------------
61C D u m m y A r g u m e n t s
62C-----------------------------------------------
63 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
64 INTEGER, INTENT(IN) :: MAT_ID
65 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
66 TYPE(submodel_data), DIMENSION(*),INTENT(IN) :: LSUBMODEL
67 TYPE(mlaw_tag_), INTENT(INOUT) :: MTAG
68 INTEGER, DIMENSION(NPROPMI) ,INTENT(INOUT) :: IPM
69 my_real, DIMENSION(NPROPM) ,INTENT(INOUT) :: pm
70 TYPE(matparam_struct_),INTENT(INOUT) :: MATPARAM
71C-----------------------------------------------
72C D e s c r i p t i o n
73C-----------------------------------------------
74C
75C INITIALISATION MATERIAU 11 LOI SPECIALE CONDITIONS AUX LIMITES
76C ISENTHALPIQUES FLUIDES PARFAIT
77C
78C ITYP = 0 - (GAS) INLET USING STAGNATION POINT STATE
79C ITYP = 1 - (LIQUID)INLET USING STAGNATION POINT STATE
80C ITYP = 2 - GENERAL INLET/OUTLET APPLYING USER FUNCTIONS
81C ITYP = 3 - NON RFLECTING FRONTIER (NRF)
82C
83C-----------------------------------------------
84C L o c a l V a r i a b l e s
85C-----------------------------------------------
86 INTEGER IFR, IFP, IFE, IFT, IFQ, ITYP, INOD, ICOMP
87 my_real
88 . gam, vcrt2, gamrp, gam1, p0, c1, ssp2, ef,
89 . rho0, psh, carl, gam2, e0, vcrt, rcrt, pcrt, ssp,
90 . t0, dc, alp0, rhoc,
91 . tscal,rhor,
92 . fac_length_, fac_time_
93 LOGICAL :: IS_ENCRYPTED, IS_AVAILABLE
94C-----------------------------------------------
95C S o u r c e L i n e s
96C-----------------------------------------------
97
98 is_encrypted = .false.
99 is_available = .false.
100 fac_length_ = one
101 fac_time_ = one
102
103C Check input encryption
104 CALL hm_option_is_encrypted(is_encrypted)
105C Initial and reference density
106 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
107 CALL hm_get_floatv('Refer_Rho' ,rhor ,is_available, lsubmodel, unitab)
108 IF (rhor == zero) THEN
109 rhor = rho0
110 ENDIF
111 pm(1) = rhor
112 pm(89) = rho0
113 ale%GLOBAL%IS_BOUNDARY_MATERIAL = .true.
114C
115 carl = zero
116 t0 = zero
117 gam = zero
118 vcrt2 = zero
119 gamrp = zero
120 gam1 = zero
121 p0 = zero
122 c1 = zero
123 ssp = zero
124 ssp2 = zero
125 ef = zero
126 e0 = zero
127 dc = zero
128 rhoc = zero
129 alp0 = infinity
130 ifr = 0
131 ifp = 0
132 ife = 0
133 ift = 0
134 ifq = 0
135 ityp = 0
136 inod = 0
137 icomp = 0
138
139C Formulation type
140 CALL hm_get_intv('Itype', ityp, is_available, lsubmodel)
141C
142 CALL hm_get_floatv('MAT_PSH', psh, is_available, lsubmodel, unitab)
143 CALL hm_get_floatv('SCALE', tscal, is_available, lsubmodel, unitab)
144 CALL hm_get_floatv_dim('SCALE', fac_time_, is_available, lsubmodel, unitab)
145 CALL hm_get_floatv_dim('h', fac_length_, is_available, lsubmodel, unitab)
146
147C Switch according to input type
148 SELECT CASE(ityp)
149 CASE (0)
150 CALL hm_get_intv('NODE1', inod, is_available, lsubmodel)
151 CALL hm_get_floatv('GAMMA', gam, is_available, lsubmodel, unitab)
152 CALL hm_get_floatv('K_cdi', dc, is_available, lsubmodel, unitab)
153
154 CALL hm_get_intv('FUN_A1', ifr, is_available, lsubmodel)
155
156 CALL hm_get_intv('FUN_A2', ifp, is_available, lsubmodel)
157 CALL hm_get_floatv('MAT_PScale', p0, is_available, lsubmodel, unitab)
158 CASE (1)
159 CALL hm_get_intv('NODE1', inod, is_available, lsubmodel)
160 CALL hm_get_floatv('MAT_C1', c1, is_available, lsubmodel, unitab)
161 CALL hm_get_floatv('K_cdi', dc, is_available, lsubmodel, unitab)
162
163 CALL hm_get_intv('FUN_A1', ifr, is_available, lsubmodel)
164
165 CALL hm_get_intv('FUN_A2', ifp, is_available, lsubmodel)
166 CALL hm_get_floatv('MAT_PScale', p0, is_available, lsubmodel, unitab)
167
168 CALL hm_get_intv('FUN_A6', ife, is_available, lsubmodel)
169 CALL hm_get_floatv('MAT_E0', e0, is_available, lsubmodel, unitab)
170 CASE (2)
171 CALL hm_get_intv('FUN_A1', ifr, is_available, lsubmodel)
172
173 CALL hm_get_intv('FUN_A2', ifp, is_available, lsubmodel)
174 CALL hm_get_floatv('MAT_PScale', p0, is_available, lsubmodel, unitab)
175
176 CALL hm_get_intv('FUN_A6', ife, is_available, lsubmodel)
177 CALL hm_get_floatv('MAT_E0', e0, is_available, lsubmodel, unitab)
178 CASE (3)
179 CALL hm_get_floatv('MAT_C0', ssp, is_available, lsubmodel, unitab)
180 CALL hm_get_floatv('h', carl, is_available, lsubmodel, unitab)
181 IF(carl*ssp==zero)THEN
182 CALL ancmsg(msgid=304,
183 . msgtype=msgerror,
184 . anmode=aninfo,
185 . i2=mat_id,
186 . c1=titr)
187 ENDIF
188 CASE DEFAULT
189 CALL ancmsg(msgid=1665, msgtype=msgerror, anmode=aninfo,
190 . i1=mat_id,
191 . i2 = ityp,
192 . c1=titr)
193 END SELECT
194
195 CALL hm_get_intv('Xt_fun', ift, is_available, lsubmodel)
196 CALL hm_get_intv('Yt_fun', ifq, is_available, lsubmodel)
197
198C----------------------
199C DEFAULT VALUES
200C----------------------
201 IF(carl==zero)carl = em01 * fac_length_
202 IF(t0==zero) t0 = three100
203 IF(tscal == zero) tscal = one * fac_time_
204 SELECT CASE(ityp)
205 CASE (0)
206 gam1 = one/(gam - one)
207 gam2 = half/gam
208 gamrp= gam2*(gam-one)
209 e0 = gam1*(p0+psh)
210 ef = gam1*(psh)
211 ssp2 = gam*p0/rho0
212 vcrt2= two*gam*(p0+psh)/rho0/(gam + one)
213 vcrt = sqrt(vcrt2)
214 rcrt = (one-gamrp*rho0*vcrt2/(p0+psh))**gam1
215 pcrt = (p0+psh)*rcrt**gam
216 rcrt = rho0*rcrt
217 CASE (1)
218 ssp2 = c1/rho0
219 gam1 = one
220 CASE (2)
221 gam1 = one
222 ssp2 = ssp**2
223 CASE (3)
224 IF(carl>zero)alp0=half*ssp/carl
225 rhoc=rho0*ssp
226 CASE DEFAULT
227C Already exited at input reading
228 END SELECT
229
230C----------------------
231C PRINTOUT
232C----------------------
233 WRITE(iout, 800)trim(titr),mat_id,11
234
235 IF(is_encrypted)THEN
236 WRITE(iout,900)
237 WRITE(iout,'(5x,a,//)')'confidential data'
238 ELSE
239 WRITE(IOUT,1000)ITYP,PSH,TSCAL
240 WRITE(IOUT, 850) RHO0,RHOR
241 SELECT CASE(ITYP)
242 CASE(0)
243 WRITE(IOUT,1100)GAM,P0,E0,VCRT,RCRT,PCRT,INOD
244 WRITE(IOUT,1350)EF,DC
245 WRITE(IOUT,1300)IFR,IFP
246 CASE(1)
247 WRITE(IOUT,1200)C1,P0,E0,INOD
248 WRITE(IOUT,1349)DC
249 WRITE(IOUT,1301)IFR,IFP,IFE
250 CASE(2)
251 WRITE(IOUT,1400)P0,E0
252 WRITE(IOUT,1302)IFR,IFP,IFE
253 CASE(3)
254 WRITE(IOUT,1353)SSP,CARL,INOD
255 CASE DEFAULT
256C Already exited at input reading
257 END SELECT
258 IF(ITYP/=3)WRITE(IOUT,2001) IFT,IFQ
259 ENDIF
260
261C----------------------
262C STORAGE
263C----------------------
264 IPM(11) = IFR
265 IPM(12) = IFP
266 IPM(13) = IFE
267 IPM(14) = IFT
268 IPM(17) = IFQ
269
270 PM(23) = E0
271 PM(25) = GAM
272 PM(27) = VCRT2
273 PM(28) = GAMRP
274 PM(29) = GAM1
275 PM(31) = P0
276 PM(32) = C1
277 PM(34) = -HALF*SSP2/CARL**2
278 PM(37) = EF
279 PM(40) = ONE/TSCAL
280 PM(50) = ITYP
281 PM(51) = INOD ! USR2SYS APRES FSDCOD, IF(IFORM8==2) INOD=USR2SYS(INOD,ITABM1,MESS)
282 PM(79) = T0
283 PM(80) = EP30
284 PM(88) = PSH
285 PM(97) = RHOC
286 PM(98) = ALP0
287 PM(99) = DC+ONE
288
289C----------------------
290C MATERIAL BUFFER
291C----------------------
292 MTAG%G_RK = 1
293 MTAG%G_RE = 1
294 MTAG%G_TEMP = 1
295 MTAG%L_RK = 1
296 MTAG%L_TEMP = 1
297 MTAG%L_RE = 1
298 MTAG%L_VK = 1
299 MTAG%L_EINS = 1
300
301 ! EOS/Thermo keyword for pressure treatment in elements
302 CALL INIT_MAT_KEYWORD(MATPARAM,"HYDRO_EOS")
303
304 CALL INIT_MAT_KEYWORD(MATPARAM,"INCOMPRESSIBLE")
305
306 ! Properties compatibility
307 CALL INIT_MAT_KEYWORD(MATPARAM,"SOLID_ISOTROPIC")
308C--------------------------------
309 RETURN
310
311C----------------------
312C FORMAT
313C----------------------
314 800 FORMAT(/
315 & 5X,A,/,
316 & 5X,'material number. . . . . . . . . . . . . . .=',I10/,
317 & 5X,'material law . . . . . . . . . . . . . . . .=',I10/)
318 850 FORMAT(
319 & 5X,'initial density . . . . . . . . . . . .=',1PG20.13/,
320 & 5X,'reference density . . . . . . . . . . .=',1PG20.13/)
321 900 FORMAT(
322 & 5X,40H LAW FOR FLUID BOUNDARY ELEMENTS ,/,
323 & 5X,40H ------------------------------- ,/)
324 1000 FORMAT(
325 & 5X,40H LAW FOR FLUID BOUNDARY ELEMENTS ,/,
326 & 5X,40H ------------------------------- ,/,
327 & 5X,40HITYP. . . . . . . . . . . . . . . . . .=,I10/,
328 & 5X,'ityp = 0 : gas inlet',/,
329 & 5X,'ityp = 1 : liquid inlet',/,
330 & 5X,'ityp = 2 : imposed inlet/outlet - time dependent',/,
331 & 5X,'ityp = 3 : non-reflecting boundary',/,
332 & 5X,'psh: pressure shift . . . . . . . . . .=',1PG20.13/,
333 & 5X,'time scale factor . . . . . . . . . . .=',1PG20.13/)
334 1100 FORMAT(
335 & 5X,40HGAMMA CONSTANT. . . . . . . . . . . . .=,E12.4/,
336 & 5X,40HSTAGNATION PRESSURE . . . . . . . . . .=,E12.4/,
337 & 5X,40HSTAGNATION ENERGY . . . . . . . . . . .=,E12.4/,
338 & 5X,40HCRITICAL VELOCITY . . . . . . . . . . .=,E12.4/,
339 & 5X,40HCRITICAL DENSITY. . . . . . . . . . . .=,E12.4/,
340 & 5X,40HCRITICAL PRESSURE . . . . . . . . . . .=,E12.4/,
341 & 5X,40HREFERENCE NODE (VELOCITY) . . . . . . .=,I10/)
342 1200 FORMAT(
343 & 5X,40HBULK MODULUS. . . . . . . . . . . . . .=,E12.4/,
344 & 5X,40HSTAGNATION PRESSURE . . . . . . . . . .=,E12.4/,
345 & 5X,40HSTAGNATION ENERGY . . . . . . . . . . .=,E12.4/,
346 & 5X,40HREFERENCE NODE (VELOCITY) . . . . . . .=,I10/)
347 1300 FORMAT(
348 & 5X,40HSTAGNATION DENSITY LOAD CURVE. . . . .=,I10/,
349 & 5X,40HSTAGNATION PRESSURE LOAD CURVE. . . . .=,I10/)
350 1301 FORMAT(
351 & 5X,40HSTAGNATION DENSITY LOAD CURVE. . . . .=,I10/,
352 & 5X,40HSTAGNATION PRESSURE LOAD CURVE. . . . .=,I10/,
353 & 5X,40HSTAGNATION ENERGY LOAD CURVE. . . . .=,I10/)
354 1302 FORMAT(
355 & 5X,40HDENSITY LOAD CURVE . . . . . . . . .=,I10/,
356 & 5X,40HPRESSURE LOAD CURVE . . . . . . . . .=,I10/,
357 & 5X,40HENERGY LOAD CURVE . . . . . . . . .=,I10/)
358 1349 FORMAT(
359 & 5X,40HDISCHARGE COEFFICIENT (ENTRY LOSS). . .=,E12.4/)
360 1350 FORMAT(
361 & 5X,40HFINAL STAGNATION ENERGY . . . . . . . .=,E12.4/,
362 & 5X,40HDISCHARGE COEFFICIENT (ENTRY LOSS). . .=,E12.4/)
363 1353 FORMAT(
364 & 5X,40HCHARACTERISTIC SOUND SPEED. . . . . . .=,E12.4/,
365 & 5X,40HCHARACTERISTIC LENGTH . . . . . . . . .=,E12.4/,
366 & 5X,40HREFERENCE NODE (VELOCITY) . . . . . . .=,I10/)
367 1400 FORMAT(
368 & 5X,40HINITIAL PRESSURE(P-PSH) . . . . . . . .=,E12.4/,
369 & 5X,40HINITIAL ENERGY. . . . . . . . . . . . .=,E12.4/)
370 2001 FORMAT(
371 & 5X,' thermal boundary ',/,
372 & 5X,' ---------------- ',/,
373 & 5X,'temperature load curve. . . . . . . . .=',I10/,
374 & 5X,'flux load curve. . . . . . . . . . . .=',I10/)
375C-----------
376 RETURN
377 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_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_mat11
subroutine hm_read_mat11(lsubmodel, mtag, unitab, ipm, pm, mat_id, titr, matparam)
Definition hm_read_mat11.F:41
ale_mod
Definition ale_mod.F:150
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:253
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
elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:235
submodel_mod::submodel_data
Definition submodel_mod.F:728