OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_interfaces.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_interfaces ../starter/source/interfaces/reader/hm_read_interfaces.F
25!||--- called by ------------------------------------------------------
26!|| lectur ../starter/source/starter/lectur.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
30!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
31!|| hm_read_inter_fsi ../starter/source/interfaces/reader/hm_read_inter_fsi.F
32!|| hm_read_inter_lagdt_type07 ../starter/source/interfaces/int07/hm_read_inter_lagdt_type07.F
33!|| hm_read_inter_lagmul ../starter/source/interfaces/reader/hm_read_inter_lagmul.F
34!|| hm_read_inter_struct ../starter/source/interfaces/reader/hm_read_inter_struct.F
35!|| hm_sz_r2r ../starter/source/coupling/rad2rad/routines_r2r.F
36!|| ifrontplus ../starter/source/spmd/node/frontplus.F
37!|| usr2sys ../starter/source/system/sysfus.F
38!||--- uses -----------------------------------------------------
39!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
40!|| message_mod ../starter/share/message_module/message_mod.F
41!|| r2r_mod ../starter/share/modules1/r2r_mod.F
42!|| submodel_mod ../starter/share/modules1/submodel_mod.F
43!||====================================================================
44 SUBROUTINE hm_read_interfaces(
45 1 IPARI ,FRIGAP ,ITAB ,ITABM1 ,
46 2 IGRNOD ,IGRSURF ,IGRSLIN ,IGRBRIC ,IGRSH3N ,
47 3 IGRTRUSS ,NPC ,ISKN ,XFILTR ,STFAC ,
48 4 FRIC_P ,I2RUPT ,AREASL ,UNITAB ,NOM_OPT ,
49 5 DEF_INTER ,NPC1 ,SENSORS ,MULTI_FVM ,NOM_OPTFRIC,
50 6 INTBUF_FRIC_TAB ,LSUBMODEL,TF ,NPTS ,NPARI ,
51 7 KLOADPINTER ,DGAPINT ,INTERFACES ,SITAB ,NPARIR ,
52 8 SITABM1 ,SISKWN ,LISKN ,SNPC ,SNPC1 ,
53 9 ITHERM_FE ,INTHEAT )
54C-----------------------------------------------
55C D e s c r i p t i o n
56C-----------------------------------------------
57C General Subroutine which manages interface Readers (/INTER/TYPE..)
58C-----------------------------------------------
59C M o d u l e s
60C-----------------------------------------------
61 USE unitab_mod , ONLY : unit_type_
62 USE r2r_mod
63 USE message_mod
64 USE multi_fvm_mod , ONLY : multi_fvm_struct
65 USE intbuf_fric_mod , ONLY : intbuf_fric_struct_
66 USE groupdef_mod , ONLY : group_, surf_
67 USE submodel_mod , ONLY : nsubmod, submodel_data
68 USE hm_option_read_mod
69 USE sensor_mod , ONLY : sensors_
70 USE interfaces_mod , ONLY : interfaces_
71 USE names_and_titles_mod , ONLY : nchartitle, ncharkey
72C-----------------------------------------------
73C I m p l i c i t T y p e s
74C-----------------------------------------------
75#include "implicit_f.inc"
76C-----------------------------------------------
77C C o m m o n B l o c k s
78C-----------------------------------------------
79#include "com01_c.inc"
80#include "com04_c.inc"
81#include "units_c.inc"
82#include "scr17_c.inc"
83#include "r2r_c.inc"
84#include "inter22.inc"
85C-----------------------------------------------
86C D u m m y A r g u m e n t s
87C-----------------------------------------------
88 INTEGER,INTENT(IN) :: NPARI, NPARIR !< array sizes (IPARI and FRIGAP)
89 INTEGER,INTENT(IN) :: NPTS !< array size TF
90 INTEGER,INTENT(IN) :: SNPC, SNPC1 !< array size (NPC and NPC1)
91 INTEGER,INTENT(IN) :: SITAB, SITABM1 !< array size (ITAB, ITABM1)
92 INTEGER,INTENT(IN) :: LISKN, SISKWN !< array size (ISKN)
93 INTEGER,INTENT(IN) :: ITHERM_FE
94 INTEGER,INTENT(INOUT) :: INTHEAT
95 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
96 INTEGER NOM_OPT(LNOPT1,*)
97 INTEGER IPARI(NPARI,NINTER), ITAB(NUMNOD), ITABM1(SITABM1),NPC(SNPC),ISKN(SISKWN),
98 . def_inter(100),npc1(snpc1),nom_optfric(lnopt1,*),kloadpinter(ninter+1)
99 my_real xfiltr(ninter),stfac(ninter),fric_p(10,ninter),i2rupt(6,ninter),frigap(nparir,ninter),areasl(ninter),tf(npts)
100 my_real,INTENT(IN) :: dgapint(ninter)
101 TYPE(multi_fvm_struct), INTENT(IN) :: MULTI_FVM
102 TYPE(intbuf_fric_struct_) INTBUF_FRIC_TAB(NINTERFRIC)
103 TYPE (SUBMODEL_DATA), DIMENSION(NSUBMOD),INTENT(IN) :: LSUBMODEL
104 TYPE (SENSORS_) ,INTENT(IN) :: SENSORS
105 TYPE (GROUP_) , DIMENSION(NGRNOD) :: IGRNOD
106 TYPE (GROUP_) , DIMENSION(NGRBRIC) :: IGRBRIC
107 TYPE (GROUP_) , DIMENSION(NGRSH3N) :: IGRSH3N
108 TYPE (GROUP_) , DIMENSION(NGRTRUS) :: IGRTRUSS
109 TYPE (SURF_) , DIMENSION(NSURF) ,TARGET :: IGRSURF
110 TYPE (SURF_) , DIMENSION(NSLIN) ,TARGET :: IGRSLIN
111 TYPE (INTERFACES_) ,INTENT(INOUT):: INTERFACES
112C-----------------------------------------------
113C L o c a l V a r i a b l e s
114C-----------------------------------------------
115 INTEGER I,J,P,NI,NIN,ILAGM,NOINT,NB,UID,SUB_ID,IFLAGUNIT,K
116 CHARACTER(LEN=NCHARKEY) :: KEY, KEY2, STR_TMP
117 CHARACTER(LEN=NCHARTITLE) :: TITR
118 CHARACTER(LEN=40) :: MESS
119 DATA mess/'INTERFACE INPUT '/
120C-----------------------------------------------
121C E x t e r n a l F u n c t i o n s
122C-----------------------------------------------
123 INTEGER,EXTERNAL :: USR2SYS
124C-----------------------------------------------
125C S o u r c e L i n e s
126C-----------------------------------------------
127C READING INTERFACES
128
129 ni = 0
130 nb = 0
131 intheat = 0
132 ftempvar21 = 0
133 nitsche = 0
134 interadhesion = 0
135 i22len = 0
136 int22 = 0
137
138 ! START BROWSING MODEL INTERFACES
139 CALL hm_option_start('/INTER')
140 ni = 0
141 ! BROWSING MODEL INTERFACES 1->HM_NINTER
142 DO nin=1,hm_ninter
143 nb = nb+1
144 !Multidomains --> ignoring untaged interface and sub-intergaces
145 IF (nsubdom > 0) THEN
146 IF(tagint(nb)==0) CALL hm_sz_r2r(tagint,nb,lsubmodel)
147 ENDIF
148 ilagm = 0
149 key = ''
150 ! EXTRACT DATAS OF /INTER/... LINE
151 CALL hm_option_read_key(lsubmodel,option_id=noint,unit_id=uid,submodel_id=sub_id,
152 . option_titr=titr,keyword2=key,keyword3=key2)
153 ! CHECK IF READ OPTION IS /INTER/SUB
154 IF(key(1:len_trim(key)) == 'SUB') cycle
155 ni=ni+1
156 DO i = 1,npari
157 ipari(i,ni) = 0
158 ENDDO
159 ipari(23,ni) = 1
160 ! CHECK IF UID IS EXISTING
161 iflagunit = 0
162 DO j=1,unitab%NUNITS
163 IF (unitab%UNIT_ID(j) == uid) THEN
164 iflagunit = 1
165 EXIT
166 ENDIF
167 ENDDO
168 IF (uid /= 0.AND.iflagunit == 0) THEN
169 CALL ancmsg(msgid=659,anmode=aninfo,msgtype=msgerror,i2=uid,i1=noint,c1='INTERFACE',c2='INTERFACE',c3=titr)
170 ENDIF
171 ! WRITE TITLE IN OUT FILE
172 WRITE(iout,1300)
173 str_tmp = key(1:len_trim(key))
174 ! CHECK INTERFACE AVAILABILITY WITH 2D ANALYSIS
175 IF(n2d /= 0 .AND.
176 . str_tmp /= 'TYPE1' .AND. str_tmp /= 'type2.AND.'
177 . STR_TMP /= 'type3.AND.' STR_TMP /= 'type5') THEN
178 CALL ANCMSG(MSGID=2097,ANMODE=ANSTOP,MSGTYPE=MSGERROR,I1=NOINT,C1=TITR,C2=STR_TMP)
179 ENDIF
180C--------------------------------------------------
181 SELECT CASE(STR_TMP)
182 CASE ('kine', 'type1', 'type9', 'type12', 'type18', 'type22')
183 ! SPECIFIC CASE OF FSI INTERFACES READING ROUTINES
184 CALL HM_READ_INTER_FSI(
185 1 IPARI(1,NI) ,STFAC(NI) ,FRIGAP(1,NI) ,XFILTR(NI) ,FRIC_P(1,NI),
186 2 IGRNOD ,IGRSURF ,ILAGM ,UNITAB ,NI ,
187 3 NOM_OPT ,TITR ,IGRBRIC ,IGRSH3N ,IGRTRUSS ,
188 4 NOINT ,KEY ,LSUBMODEL ,ITAB ,
189 5 ITABM1 ,ISKN ,MULTI_FVM ,NPARI ,NPARIR ,
190 6 NGRNOD ,NGRBRIC ,NGRSH3N ,NGRTRUS ,NSURF ,
191 7 LNOPT1 ,SITABM1 ,SISKWN ,LISKN ,
192 8 LTITR ,SITAB)
193
194 CASE ('lagmul')
195 ! /INTER/LAGMUL
196 CALL HM_READ_INTER_LAGMUL(
197 1 IPARI(1,NI) ,STFAC(NI) ,FRIGAP(1,NI) ,IGRNOD ,IGRSURF ,
198 2 ILAGM ,NI ,NOM_OPT ,TITR ,DEF_INTER ,
199 3 IGRBRIC ,UNITAB ,LSUBMODEL ,NOINT ,KEY2 ,
200 4 NPARI ,NPARIR)
201
202 CASE ('lagdt')
203 ! /INTER/LAGDT
204 CALL HM_READ_INTER_LAGDT_TYPE07(
205 1 IPARI(1,NI) ,STFAC(NI) ,FRIGAP(1,NI) ,XFILTR(NI),FRIC_P(1,NI),
206 2 IGRNOD ,IGRSURF ,ILAGM ,UNITAB ,NI ,
207 4 NOM_OPT ,TITR ,DEF_INTER ,NOINT ,LSUBMODEL ,
208 5 NPARI ,NPARIR)
209
210 CASE DEFAULT
211 ! ALL OTHER INTERFACES : STRUCTURE + TIED TYPE2
212 CALL HM_READ_INTER_STRUCT(
213 1 IPARI(1,NI) ,STFAC(NI) ,FRIGAP(1,NI) ,XFILTR(NI) ,FRIC_P(1,NI),
214 2 IGRNOD ,IGRSURF ,IGRSLIN ,NPC ,ILAGM ,
215 3 UNITAB ,I2RUPT(1,NI) ,
216 4 AREASL(NI) ,NI ,NOM_OPT ,TITR ,DEF_INTER ,
217 5 NPC1 ,SENSORS ,
218 6 NOM_OPTFRIC ,INTBUF_FRIC_TAB,IGRBRIC ,NOINT ,KEY ,
219 7 LSUBMODEL ,TF ,INTERFACES ,SNPC ,SNPC1 ,
220 8 NPARI ,NPARIR ,NPTS ,
221 9 LNOPT1 ,LTITR ,NINTER25 ,
222 . NGRNOD ,NGRBRIC ,NSURF ,NSLIN ,
223 . ITHERM_FE ,INTHEAT )
224
225 END SELECT
226
227 NOINT = IPARI(15,NI)
228 IF (IPARI(71,NI) <= 0) THEN
229 DO K=1,NI-1
230 IF (NOINT == IPARI(15,K)) THEN
231 CALL ANCMSG(MSGID=117,MSGTYPE=MSGERROR,ANMODE=ANINFO_BLIND_1,I1=NOINT,C1=TITR,I2=NOINT)
232 ENDIF
233 ENDDO
234 ENDIF
235
236 IF(IPARI(39,NI) == 1)THEN
237 IPARI(40,NI)=USR2SYS(IPARI(40,NI),ITABM1,MESS,IPARI(15,NI))
238 DO P = 1, NSPMD
239 CALL IFRONTPLUS(IPARI(40,NI),P)
240 END DO
241 ELSEIF(IPARI(39,NI) == 2)THEN
242 IPARI(40,NI)=USR2SYS(IPARI(40,NI),ITABM1,MESS,IPARI(15,NI))
243 IPARI(41,NI)=USR2SYS(IPARI(41,NI),ITABM1,MESS,IPARI(15,NI))
244 DO P = 1, NSPMD
245 CALL IFRONTPLUS(IPARI(40,NI),P)
246 CALL IFRONTPLUS(IPARI(41,NI),P)
247 END DO
248 ENDIF
249
250 IF(NINTLOADP > 0)IPARI(95,NI) = KLOADPINTER(NI+1) - KLOADPINTER(NI)
251 IF(NINTLOADP > 0)FRIGAP(46,NI) = DGAPINT(NI)
252
253 ENDDO ! NI
254
255C=======================================================================
256 1300 FORMAT( /1X,' interfaces ' /
257 . 1X,' -------------- '// )
258 RETURN
259 END
260
my_real
#define my_real
Definition cppsort.cpp:32
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:251
hm_read_interfaces
subroutine hm_read_interfaces(ipari, frigap, itab, itabm1, igrnod, igrsurf, igrslin, igrbric, igrsh3n, igrtruss, npc, iskn, xfiltr, stfac, fric_p, i2rupt, areasl, unitab, nom_opt, def_inter, npc1, sensors, multi_fvm, nom_optfric, intbuf_fric_tab, lsubmodel, tf, npts, npari, kloadpinter, dgapint, interfaces, sitab, nparir, sitabm1, siskwn, liskn, snpc, snpc1, itherm_fe, intheat)
Definition hm_read_interfaces.F:54
i2rupt
subroutine i2rupt(x, v, a, ms, in, stifn, fsav, weight, irect, nsv, msr, irtl, irupt, crst, mmass, miner, smass, siner, area, uvar, xsm0, dsm, fsm, prop, ipari, nsn, nmn, nuvar, igtyp, pid, npf, tf, itab, fncont, pdama2, isym, inorm, h3d_data, fncontp, ftcontp)
Definition int2rupt.F:122
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:410
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
names_and_titles_mod::ncharkey
integer, parameter ncharkey
Definition names_and_titles_mod.F:1000
r2r_mod
Definition r2r_mod.F:125
r2r_mod::tagint
integer, dimension(:), allocatable tagint
Definition r2r_mod.F:132
submodel_mod
Definition submodel_mod.F:712
submodel_mod::nsubmod
integer nsubmod
Definition submodel_mod.F:720
unitab_mod
Definition unitab_mod.F:513
hm_sz_r2r
subroutine hm_sz_r2r(tag, val, lsubmodel)
Definition routines_r2r.F:300
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
groupdef_mod::group_
Definition groupdef_mod.F:700
groupdef_mod::surf_
Definition groupdef_mod.F:745
submodel_mod::submodel_data
Definition submodel_mod.F:723
unitab_mod::unit_type_
Definition unitab_mod.F:520