OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_radiation.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_radiation ../starter/source/loads/thermic/hm_read_radiation.F
25!||--- called by ------------------------------------------------------
26!|| lectur ../starter/source/starter/lectur.F
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_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
33!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
34!|| usr2sys ../starter/source/system/sysfus.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_radiation(IB, FAC, ITAB, IXS ,IGRSURF, UNITAB ,LSUBMODEL,
41 . NRADIA,NUMRADIA,NIRADIA,LFACTHER)
42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE unitab_mod
46 USE message_mod
47 USE groupdef_mod
48 USE submodel_mod
49 USE hm_option_read_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 "param_c.inc"
59#include "units_c.inc"
60#include "com04_c.inc"
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 INTEGER ,INTENT(IN) :: NRADIA
65 INTEGER ,INTENT(IN) :: NUMRADIA
66 INTEGER ,INTENT(IN) :: NIRADIA
67 INTEGER ,INTENT(IN) :: LFACTHER
68 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
69 INTEGER IB(NIRADIA,*), ITAB(*), IXS(NIXS,*)
70 my_real fac(lfacther,*)
71
72 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
73 TYPE(submodel_data),INTENT(IN)::LSUBMODEL(*)
74C-----------------------------------------------
75C L o c a l V a r i a b l e s
76C-----------------------------------------------
77 INTEGER K, M, I1, I2, I3, I4, IFU, I, ISENS,NPR0,NN,ISU,IS,
78 . id,j,uid,iflagunit,ity
79 my_real
80 . fcx, fcy, fac_l, fac_t, fac_m, emi, sigma, tstart, tstop,
81 . emiss(numradia),fcx_dim,fcy_dim,tstop_dim
82 CHARACTER MESS*40
83 CHARACTER(LEN=NCHARTITLE) :: TITR
84 LOGICAL IS_AVAILABLE
85C-----------------------------------------------
86C E x t e r n a l F u n c t i o n s
87C-----------------------------------------------
88 INTEGER USR2SYS
89 DATA mess/'RADIATIVE FLUX DEFINITION '/
90C=======================================================================
91 is_available = .false.
92 k = 0
93C--------------------------------------------------
94C START BROWSING MODEL RADIATION
95C--------------------------------------------------
96 CALL hm_option_start('/RADIATION')
97C--------------------------------------------------
98C BROWSING /RADIATION OPTIONS 1->NRADIA
99C--------------------------------------------------
100 DO i=1,nradia
101 titr = ''
102 CALL hm_option_read_key(lsubmodel,
103 . unit_id = uid,
104 . option_id = id,
105 . option_titr = titr)
106 iflagunit = 0
107 DO j=1,unitab%NUNITS
108 IF (unitab%UNIT_ID(j) == uid) THEN
109 iflagunit = 1
110 EXIT
111 ENDIF
112 ENDDO
113 IF (uid /= 0.AND.iflagunit == 0) THEN
114 CALL ancmsg(msgid=659,anmode=aninfo,msgtype=msgerror,
115 . i2=uid,i1=id,c1='CONVECTION HEAT',
116 . c2='CONVECTION HEAT',
117 . c3=titr)
118 ENDIF
119C--------------------------------------------------
120C EXTRACT DATAS (INTEGER VALUES)
121C--------------------------------------------------
122 CALL hm_get_intv('entityid',isu,is_available,lsubmodel)
123 CALL hm_get_intv('curveid',ifu,is_available,lsubmodel)
124 CALL hm_get_intv('rad_sensor_id',isens,is_available,lsubmodel)
125C--------------------------------------------------
126C EXTRACT DATAS (REAL VALUES)
127C--------------------------------------------------
128 CALL hm_get_floatv('xscale',fcx,is_available,lsubmodel,unitab)
129 CALL hm_get_floatv_dim('xscale',fcx_dim,is_available,lsubmodel,unitab)
130 CALL hm_get_floatv('magnitude',fcy,is_available,lsubmodel,unitab)
131 CALL hm_get_floatv_dim('magnitude',FCY_DIM,IS_AVAILABLE,LSUBMODEL,UNITAB)
132 CALL HM_GET_FLOATV('rad_tstart',TSTART,IS_AVAILABLE,LSUBMODEL,UNITAB)
133 CALL HM_GET_FLOATV('rad_tstop',TSTOP,IS_AVAILABLE,LSUBMODEL,UNITAB)
134 CALL HM_GET_FLOATV_DIM('rad_tstop',TSTOP_DIM,IS_AVAILABLE,LSUBMODEL,UNITAB)
135 CALL HM_GET_FLOATV('flux',EMI,IS_AVAILABLE,LSUBMODEL,UNITAB)
136C--------------------------------------------------
137 IF (FCX == ZERO) FCX = FCX_DIM
138 IF (FCY == ZERO) FCY = FCY_DIM
139 IF(TSTOP == ZERO) TSTOP= EP30 * TSTOP_DIM
140C
141 FAC_M = UNITAB%FAC_M_WORK
142 FAC_T = UNITAB%FAC_T_WORK
143 SIGMA=STEFBOLTZ*(FAC_T*FAC_T*FAC_T)/FAC_M
144C
145 IS=0
146 DO J=1,NSURF
147 IF (ISU == IGRSURF(J)%ID) IS=J
148 ENDDO
149 IF(IS > 0)THEN
150 NN=IGRSURF(IS)%NSEG
151 DO J=1,NN
152 K=K+1
153 IB(1,K)=IGRSURF(IS)%NODES(J,1)
154 IB(2,K)=IGRSURF(IS)%NODES(J,2)
155 IB(3,K)=IGRSURF(IS)%NODES(J,3)
156 ITY =IGRSURF(IS)%ELTYP(J)
157 IF(ITY==7)THEN
158C true triangles (not segments built from 3 nodes).
159 IB(4,K)=0
160 ELSE
161 IB(4,K)=IGRSURF(IS)%NODES(J,4)
162 ENDIF
163 IB(5,K) = IFU
164 IB(6,K) = ISENS
165 IB(7,K) = IGRSURF(IS)%ELTYP(J)
166 IB(8,K) = IGRSURF(IS)%ELEM(J)
167 IF(ITY == 1) THEN
168 IB(9,K) = IXS(11,IGRSURF(IS)%ELEM(J))
169 ELSE
170 IB(9,K) = 0
171 ENDIF
172C
173 FAC(1,K) = FCY
174 FAC(2,K) = ONE/FCX
175 FAC(3,K) = EMI*SIGMA
176 FAC(4,K) = TSTART
177 FAC(5,K) = TSTOP
178 FAC(6,K) = ONE
179C
180C temporary storage for print out.
181 EMISS(K)=EMI
182 ENDDO
183 ENDIF
184 ENDDO
185C
186 I1=1
187 I2=MIN0(50,NUMRADIA)
188C
189 90 WRITE (IOUT,2000)
190 WRITE (IOUT,2001)
191 DO I=I1,I2
192 WRITE (IOUT,'(5(1x,i10),1x,1f10.3,2(1x,i10),1x,4g20.13)') I,
193 . ITAB(IB(1,I)),ITAB(IB(2,I)),ITAB(IB(3,I)),ITAB(IB(4,I)),
194 . EMISS(I),IB(5,I),IB(6,I),FAC(4,I),FAC(5,I),ONE/FAC(2,I),
195 . FAC(1,I)
196 ENDDO
197 IF(I2 == NUMRADIA)GOTO 200
198 I1=I1+50
199 I2=MIN0(I2+50,NUMRADIA)
200 GOTO 90
201 200 RETURN
202 300 CALL ANCMSG(MSGID=157,
203 . MSGTYPE=MSGERROR,
204 . ANMODE=ANINFO,
205 . I1=K)
206C---
207 2000 FORMAT(//
208 .' radiation heat '/
209 .' ---------------- ')
210 2001 FORMAT(/
211 .' segment node1 node2 node3 node4 emissivity',
212 .' curve sensor t-start t-stop', 8X,
213 .' scale-x scale-y')
214C-----------
215 RETURN
216 END
217
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_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:159
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:251
hm_read_radiation
subroutine hm_read_radiation(ib, fac, itab, ixs, igrsurf, unitab, lsubmodel, nradia, numradia, niradia, lfacther)
Definition hm_read_radiation.F:42
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:410
id
initmumps id
Definition multiplerhs_example.m:4
groupdef_mod
Definition groupdef_mod.F:662
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
radiation
subroutine radiation(ibcr, fradia, npc, tf, x, temp, nsensor, sensor_tab, fthe, iad, fthesky, python, glob_therm)
Definition radiation.F:38
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
submodel_mod::submodel_data
Definition submodel_mod.F:723