OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_therm.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!||====================================================================
25!|| hm_read_therm ../starter/source/materials/therm/hm_read_therm.F
26!||--- called by ------------------------------------------------------
27!|| read_material_models ../starter/source/materials/read_material_models.F
28!||--- calls -----------------------------------------------------
29!|| ancmsg ../starter/source/output/message/message.F
30!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
31!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
32!|| hm_option_count ../starter/source/devtools/hm_reader/hm_option_count.F
33!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.f
34!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
35!||--- uses -----------------------------------------------------
36!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
37!|| message_mod ../starter/share/message_module/message_mod.F
38!|| submodel_mod ../starter/share/modules1/submodel_mod.F
39!||====================================================================
40 SUBROUTINE hm_read_therm(MAT_PARAM ,BUFMAT ,BUFLEN ,IPM ,PM ,
41 . UNITAB ,LSUBMODEL )
42C-----------------------------------------------
43C D e s c r i p t i o n
44C-----------------------------------------------
45C Material Parameters (hm_read_therm.F)
46C PM(71) = JTHE is a material flag. 1: thermal parameters are input 0: no thermal parameter
47C
48C GLOBAL FLAGS (hm_read_part.F)
49C ITHERM_FE : 1 if at least one PART has thermal parameter and lagrangian framework
50C ITHERM : 1 if at least one PART has thermal parameter and ale/euler framework
51C
52C GROUP FLAG (./source/element/....tails.F) IPARG(13) will be flag for thermics.
53C IN ALL CASES
54C IPARG(13) : 0 if current group do not require temperature calculation
55C
56C SOLIDS
57C IPARG(13) : -1 if current group requires temperature calculation at nodes (FEM with Lagrange)
58C IPARG(13) : +1 if current group requires temperature calculation at centroids (FVM with ALE or Euler)
59C
60C SHELLS
61C IPARG(13) : +1 if current group requires temperature calculation at nodes (FEM with Lagrange)
62C-----------------------------------------------
63C M o d u l e s
64C-----------------------------------------------
65 USE unitab_mod
66 USE message_mod
67 USE submodel_mod
68 USE hm_option_read_mod
69 USE names_and_titles_mod , ONLY : nchartitle
70 USE matparam_def_mod
71C============================================================================
72C I m p l i c i t T y p e s
73C-----------------------------------------------
74#include "implicit_f.inc"
75C-----------------------------------------------
76C C o m m o n B l o c k s
77C-----------------------------------------------
78#include "scr17_c.inc"
79#include "com01_c.inc"
80#include "com04_c.inc"
81#include "param_c.inc"
82#include "units_c.inc"
83C-----------------------------------------------
84C D u m m y A r g u m e n t s
85C-----------------------------------------------
86 INTEGER,INTENT(IN) :: BUFLEN
87 INTEGER, DIMENSION(NPROPMI,NUMMAT), INTENT(INOUT) :: IPM
88 my_real, DIMENSION(NPROPM ,NUMMAT), INTENT(INOUT) :: pm
89 my_real, DIMENSION(BUFLEN), INTENT(INOUT) :: bufmat
90 TYPE (UNIT_TYPE_) ,INTENT(IN) :: UNITAB
91 TYPE(submodel_data) ,DIMENSION(NSUBMOD) ,INTENT(IN) :: LSUBMODEL
92 TYPE(matparam_struct_) ,DIMENSION(NUMMAT) ,INTENT(INOUT) :: MAT_PARAM
93C-----------------------------------------------
94C L o c a l V a r i a b l e s
95C-----------------------------------------------
96 INTEGER :: I,ITH,FLAGMAT,FLAGUNIT,IUNIT,UID,MAT_ID,IMAT,ILAW,ALE,
97 . lag,eul,jale,jtur,iad_therm,ntherm,iform,law2_form
98 my_real :: tini,tmelt,rho_cp,as,bs,al,bl,rho_cpm1,efrac
99 CHARACTER(LEN=NCHARTITLE) :: TITR
100 CHARACTER :: KEY*80
101 LOGICAL IS_AVAILABLE,IS_AVAILABLE_IFORM
102C=======================================================================
103c COUNT /HEAT MODELS USING CFG FILES
104c--------------------------------------------------
105c
106 CALL hm_option_count('/HEAT',ntherm)
107c
108c--------------------------------------------------
109c START BROWSING /HEAT MODELS
110c--------------------------------------------------
111c
112 CALL hm_option_start('/HEAT')
113c
114c--------------------------------------------------
115 DO ith = 1,ntherm
116c
117 CALL hm_option_read_key(lsubmodel, option_id=mat_id , option_titr=titr , unit_id=uid , keyword2=key )
118c--------------------------------------------------
119c Check MAT_Id
120c--------------------------------------------------
121 flagmat = 0
122 ilaw = 0
123 imat = 0
124 DO i =1,nummat-1
125 IF (mat_id == mat_param(i)%MAT_ID) THEN
126 flagmat = 1
127 imat = i
128 EXIT
129 ENDIF
130 ENDDO
131 IF (mat_id > 0 .AND. flagmat == 0) THEN
132 CALL ancmsg(msgid=1663,anmode=aninfo,msgtype=msgerror,
133 & i1= mat_id, c1='HEAT/MAT', c2='HEAT/MAT', c3='')
134 ENDIF
135 ilaw = mat_param(imat)%ILAW
136c--------------------------------------------------
137c Check Unit_ID
138c--------------------------------------------------
139 flagunit = 0
140 DO iunit=1,unitab%NUNITS
141 IF (unitab%UNIT_ID(iunit) == uid) THEN
142 flagunit = 1
143 EXIT
144 ENDIF
145 ENDDO
146 IF (uid > 0 .AND. flagunit == 0) THEN
147 CALL ancmsg(msgid=659, anmode=aninfo, msgtype=msgerror,
148 & i1= mat_id,
149 & i2= uid,
150 & c1='HEAT/MAT',
151 & c2='HEAT/MAT',
152 & c3= '')
153 ENDIF
154c--------------------------------------------------
155c Reading user parameters
156c--------------------------------------------------
157 CALL hm_get_floatv('HEAT_T0' ,tini ,is_available, lsubmodel, unitab)
158 CALL hm_get_floatv('HEAT_RHocp' ,rho_cp ,is_available, lsubmodel, unitab)
159 CALL hm_get_floatv('HEAT_AS' ,as ,is_available, lsubmodel, unitab)
160 CALL hm_get_floatv('HEAT_BS' ,bs ,is_available, lsubmodel, unitab)
161 CALL hm_get_intv ('HEAT_Iform' ,iform ,is_available_iform, lsubmodel)
162c
163 CALL hm_get_floatv('HEAT_T1' ,tmelt ,is_available, lsubmodel, unitab)
164 CALL hm_get_floatv('HEAT_AL' ,al ,is_available, lsubmodel, unitab)
165 CALL hm_get_floatv('HEAT_BL' ,bl ,is_available, lsubmodel, unitab)
166 CALL hm_get_floatv('HEAT_EFRAC' ,efrac ,is_available, lsubmodel, unitab)
167c--------------------------------------------------
168c Specific Case : law2
169c--------------------------------------------------
170 IF (ilaw == 2) THEN
171 law2_form = nint(pm(50,imat))
172 IF (tmelt /= pm(54,imat) .AND. law2_form == 1) THEN ! zerilli
173 CALL ancmsg(msgid=764, msgtype=msgwarning, anmode=aninfo_blind_1,
174 & i1=mat_id, c1=titr, i2=mat_id, i3=mat_id)
175 END IF
176 IF (rho_cp /= pm(69,imat)) THEN
177 CALL ancmsg(msgid=765, msgtype=msgwarning, anmode=aninfo_blind_1,
178 & i1=mat_id, c3=titr, i2=mat_id, i3=mat_id)
179 ENDIF
180c--------------------------------------------------
181c Specific Case : law73
182c--------------------------------------------------
183 ELSEIF (ilaw == 73)THEN
184 iad_therm = ipm(7,imat)-1
185 IF (bufmat(iad_therm+20) /= pm(79,imat) ) THEN
186 CALL ancmsg(msgid=764, msgtype=msgwarning, anmode=aninfo_blind_1,
187 & i1=mat_id, c1=titr, i2=mat_id, i3=mat_id)
188 ENDIF
189 IF (rho_cp == zero) THEN
190 rho_cpm1 = zero
191 ELSE
192 rho_cpm1 = one/rho_cp
193 ENDIF
194 IF (bufmat(iad_therm+21) /= rho_cpm1) THEN
195 CALL ancmsg(msgid=765, msgtype=msgwarning, anmode=aninfo_blind_1,
196 & i1=mat_id, c1=titr, i2=mat_id, i3=mat_id)
197 ENDIF
198 ENDIF
199c-------
200c Check at least if the corresponding material is indeed lagrangian !!
201 lag = 0
202 eul = 0
203 ale = 0
204 jale = nint(pm(72,imat))
205 IF (jale == 0 .AND. ilaw/=18 .AND. ilaw/=11) THEN
206 ilag= 1
207 lag = 1
208 ELSEIF(jale == 1)THEN
209 iale= 1
210 ale = 1
211 ELSEIF(jale == 2)THEN
212 ieuler= 1
213 eul = 1
214 ELSEIF(jale == 3) THEN
215 ilag= 1
216 lag = 1
217 ENDIF
218c
219 IF(is_available_iform)THEN
220 IF(iform /=0 .AND. lag == 0)THEN
221 iform = 0
222 CALL ancmsg(msgid=1609, msgtype=msgwarning, anmode=aninfo_blind_1,
223 & i1=mat_id, c1="warning", C2=TITR,
224 & C3="iform = 1 is incompatible with /ale and /euler materials, it will be ignored")
225 ENDIF
226 ENDIF
227c
228.AND. IF (RHO_CP == ZERO LAG == 0) THEN
229 CALL ANCMSG(MSGID=1609, MSGTYPE=MSGERROR, ANMODE=ANINFO_BLIND_1,
230 & I1=MAT_ID, C1="error", C2=TITR,
231 & C3="rho_0 cp PARAMETER must be greater than zero")
232 ENDIF
233c--------------------------------------------------
234c Default Values
235c--------------------------------------------------
236 IF (TINI == ZERO) TINI = PM(23,IMAT) / MAX(EM20,RHO_CP) ! E0 / rho_cp
237 IF (TINI == ZERO) TINI = THREE100
238 IF (TMELT == ZERO) TMELT = EP20
239 IF (EFRAC < ZERO) EFRAC = ZERO
240 IF (EFRAC > ONE ) EFRAC = ONE
241 IF (EFRAC == ZERO) EFRAC = ONE
242
243 PM(71,IMAT) = ONEP1 ! jthe
244 PM(69,IMAT) = RHO_CP
245 PM(75,IMAT) = AS
246 PM(76,IMAT) = BS
247 PM(77,IMAT) = AL
248 PM(78,IMAT) = BL
249 PM(79,IMAT) = TINI
250 PM(80,IMAT) = TMELT
251 PM(90,IMAT) = EFRAC
252
253 MAT_PARAM(IMAT)%ITHERM = 1 ! setting jthe=1
254 MAT_PARAM(IMAT)%THERM%RHOCP = RHO_CP
255 MAT_PARAM(IMAT)%THERM%TINI = TINI
256 MAT_PARAM(IMAT)%THERM%TMELT = TMELT
257 MAT_PARAM(IMAT)%THERM%AS = AS
258 MAT_PARAM(IMAT)%THERM%BS = BS
259 MAT_PARAM(IMAT)%THERM%AL = AL
260 MAT_PARAM(IMAT)%THERM%BL = BL
261 MAT_PARAM(IMAT)%THERM%EFRAC = EFRAC
262c--------------------------------------------------
263c Output
264c--------------------------------------------------
265 WRITE(IOUT,2000) MAT_ID,TINI,TMELT,RHO_CP,AS,BS,AL,BL,EFRAC
266C
267 ENDDO ! I = 1,NTHERM
268
269c-----------------------------------------
270 2000 FORMAT(/
271 & 5X,' THERMAL PARAMETERS ',/,
272 & 5X,' ------------------ ',/,
273 & 5X,'MATERIAL NUMBER . . . . . . . . . . . . . . . =',I10/,
274 & 5X,'T0 (INITIAL TEMPERATURE). . . . . . . . . . . =',1PG20.13/,
275 & 5X,'TMELT (MELTING TEMPERATURE) . . . . . . . . . =',1PG20.13/,
276 & 5X,'SPECIFIC HEAT . . . . . . . . . . . . . . . . =',1PG20.13/,
277 & 5X,'AS (SOLID PHASE). . . . . . . . . . . . . . . =',1PG20.13/,
278 & 5X,'BS (SOLID PHASE). . . . . . . . . . . . . . . =',1PG20.13/,
279 & 5X,'AL (LIQUID PHASE) . . . . . . . . . . . . . . =',1PG20.13/,
280 & 5X,'BL (LIQUID PHASE) . . . . . . . . . . . . . . =',1PG20.13/,
281 & 5X,'FRACTION OF STRAIN ENERGY CONVERTED INTO HEAT =',1PG20.13/)
282c-----------------------------------------
283 RETURN
284 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_count
subroutine hm_option_count(entity_type, hm_option_number)
Definition hm_option_count.F:105
hm_option_read_key
subroutine hm_option_read_key(lsubmodel, option_id, unit_id, submodel_index, submodel_id, option_titr, keyword1, keyword2, keyword3, keyword4, opt_pos)
Definition hm_option_read_key.F:249
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:251
hm_read_therm
subroutine hm_read_therm(mat_param, bufmat, buflen, ipm, pm, unitab, lsubmodel)
Definition hm_read_therm.F:42
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:410
hm_option_read_mod
Definition hm_option_read_mod.F:408
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
submodel_mod::submodel_data
Definition submodel_mod.F:723