OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_inter_type21.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_inter_type21 ../starter/source/interfaces/int21/hm_read_inter_type21.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_inter_struct ../starter/source/interfaces/reader/hm_read_inter_struct.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!|| inter_dcod_friction ../starter/source/interfaces/reader/inter_dcod_friction.F
33!|| inter_dcod_function ../starter/source/interfaces/reader/inter_dcod_function.F
34!|| inter_dcod_sensor ../starter/source/interfaces/reader/inter_dcod_sensor.F
35!|| ngr2usr ../starter/source/system/nintrr.F
36!||--- uses -----------------------------------------------------
37!|| message_mod ../starter/share/message_module/message_mod.F
38!|| submodel_mod ../starter/share/modules1/submodel_mod.F
39!||====================================================================
40 SUBROUTINE hm_read_inter_type21(
41 1 IPARI ,STFAC ,FRIGAP ,NOINT ,NI ,
42 2 IGRSURF ,XFILTR ,FRIC_P ,NPC1 ,SENSORS ,
43 3 NOM_OPT ,UNITAB ,LSUBMODEL ,TITR ,NPC ,
44 4 TF ,NPARI ,NPARIR ,SNPC ,SNPC1 ,
45 5 LNOPT1 ,ITHERM_FE ,INTHEAT ,NOM_OPTFRIC ,INTBUF_FRIC_TAB )
46C============================================================================
47C
48C-----------------------------------------------
49C M o d u l e s
50C-----------------------------------------------
51 USE message_mod
52 USE intbuf_fric_mod
53 USE groupdef_mod
54 USE submodel_mod , ONLY : nsubmod, submodel_data
55 USE unitab_mod
56 USE sensor_mod
57 USE intbufdef_mod
58 USE names_and_titles_mod , ONLY : ncharkey, nchartitle
59C-----------------------------------------------
60C I m p l i c i t T y p e s
61C-----------------------------------------------
62#include "implicit_f.inc"
63C-----------------------------------------------
64C D u m m y A r g u m e n t s
65C-----------------------------------------------
66 INTEGER,INTENT(IN) :: NPARI,NPARIR,SNPC,SNPC1,LNOPT1
67 INTEGER,INTENT(IN) :: ITHERM_FE
68 INTEGER,INTENT(INOUT) :: INTHEAT
69 INTEGER,INTENT(IN) :: NOM_OPTFRIC(LNOPT1,NINTERFRIC)
70 INTEGER NOM_OPT(LNOPT1,*)
71 INTEGER ISU1, ISU2, NI
72 INTEGER IPARI(NPARI),NPC1(SNPC1),NPC(SNPC)
73 my_real stfac,xfiltr
74 my_real frigap(nparir),fric_p(10),tf(*)
75 CHARACTER(LEN=NCHARTITLE),INTENT(IN) :: TITR
76C-----------------------------------------------
77 TYPE (surf_) ,TARGET , DIMENSION(NSURF) :: igrsurf
78 TYPE(submodel_data),INTENT(IN) :: LSUBMODEL(NSUBMOD)
79 TYPE (UNIT_TYPE_),INTENT(IN) :: UNITAB
80 TYPE (SENSORS_) ,INTENT(IN) :: SENSORS
81 TYPE(INTBUF_FRIC_STRUCT_),INTENT(INOUT) :: INTBUF_FRIC_TAB(NINTERFRIC)
82C-----------------------------------------------
83C C o m m o n B l o c k s
84C-----------------------------------------------
85#include "scr06_c.inc"
86#include "com01_c.inc"
87#include "com04_c.inc"
88#include "intstamp_c.inc"
89#include "units_c.inc"
90#include "scr12_c.inc"
91#include "remesh_c.inc"
92C-----------------------------------------------
93C L o c a l V a r i a b l e s
94C-----------------------------------------------
95
96 INTEGER I,IBC1, IBC2, IBC3, NOINT, NTYP,INACTI, IBC1M,
97 . IBC2M, IBC3M, IGSTI,IS1, IS2,ILEV, IGAP,MULTIMP,
98 . MFROT,IFQ,MODFR,IADM,NRADM,INTTH,IFORM,IKTHE,
99 . IDEL21,IDSENS,IDELKEEP,INTKG,INVN,
100 . IFRICTH,IFTLIM,FCOND,INTFRIC
101 my_real
102 . fric,gap,startt,bumult,stopt,fheat,tint,padm,
103 . angladm, cadm, depth,c1,c2,c3,c4,c5,c6, alpha,gapscale,gapmax,
104 . stmin,stmax, pmax, kthe, xthe, frad, drad,
105 . visc,xfricth,pskid,dcond,xthe_unit,kthe_unit
106 INTEGER, DIMENSION(:), POINTER :: INGR2USR
107 LOGICAL IS_AVAILABLE
108C-----------------------------------------------
109C E x t e r n a l F u n c t i o n s
110C-----------------------------------------------
111 INTEGER NGR2USR
112CC=======================================================================
113C READING PENALTY INTERFACE /INTER/TYPE21
114C=======================================================================
115
116C Initializations
117 is1=0
118 is2=0
119 ibc1=0
120 ibc2=0
121 ibc3=0
122 ibc1m=0
123 ibc2m=0
124 ibc3m=0
125 mfrot=0
126 ifq=0
127 igsti = 0
128 invn = 0
129 iftlim=1
130 iadm =0
131 nradm=1
132 inacti = 0
133 ilev = 0
134 idsens = 0
135 idel21= 0
136 idelkeep=0
137 intth = 0
138 ifricth = 0
139 ikthe= 0
140 iform = 0
141 intkg = 0
142 fcond = 0
143 igap = 0
144 multimp = 0
145 intfric = 0
146C
147 padm =one
148 angladm=zero
149 stopt=ep30
150 fric = zero
151 gap = zero
152 gapscale = zero
153 startt = zero
154 visc = zero
155 xfiltr = zero
156
157 DO i = 1, 10
158 fric_p(i) = zero
159 ENDDO
160
161 kthe = zero
162 xthe = zero
163 tint = zero
164 frad = zero
165 drad = zero
166 fheat= zero
167 c1=zero
168 c2=zero
169 c3=zero
170 c4=zero
171 c5=zero
172 c6=zero
173 gapmax=ep30
174 gapscale = zero
175 dcond = zero
176 depth = zero
177C
178 nintstamp=nintstamp+1
179 ipari(43)=nintstamp
180C
181 ntyp=21
182 ipari(15)=noint
183 ipari(7)=ntyp
184C
185 is_available = .false.
186C--------------------------------------------------
187C EXTRACT DATAS (INTEGER VALUES)
188C--------------------------------------------------
189 CALL hm_get_intv('secondaryentityids',isu1,is_available,lsubmodel)
190 CALL hm_get_intv('mainentityids',isu2,is_available,lsubmodel)
191 CALL hm_get_intv('Istf',igsti,is_available,lsubmodel)
192 CALL hm_get_intv('Ithe',intth,is_available,lsubmodel)
193 CALL hm_get_intv('Igap',igap,is_available,lsubmodel)
194 CALL hm_get_intv('Multimp',multimp,is_available,lsubmodel)
195 CALL hm_get_intv('TYPE21_Idel',idel21,is_available,lsubmodel)
196 CALL hm_get_intv('TYPE21_Invn',invn,is_available,lsubmodel)
197 CALL hm_get_intv('Iadm',iadm,is_available,lsubmodel)
198C
199 IF(igap==1.OR.igap==2) CALL hm_get_intv('TYPE21_ITim',iftlim,is_available,lsubmodel)
200C
201 CALL hm_get_intv('Deactivate_X_BC',ibc1,is_available,lsubmodel)
202 CALL hm_get_intv('Deactivate_Y_BC',ibc2,is_available,lsubmodel)
203 CALL hm_get_intv('Deactivate_Z_BC',ibc3,is_available,lsubmodel)
204 CALL hm_get_intv('INACTIV',inacti,is_available,lsubmodel)
205C
206 CALL hm_get_intv('Ifric',mfrot,is_available,lsubmodel)
207 CALL hm_get_intv('Ifiltr',ifq,is_available,lsubmodel)
208 CALL hm_get_intv('ISENSOR',idsens,is_available,lsubmodel)
209 CALL hm_get_intv('Crx_Fun',ifricth,is_available,lsubmodel)
210 CALL hm_get_intv('Ifric',mfrot,is_available,lsubmodel)
211 CALL hm_get_intv('Fric_ID',intfric,is_available,lsubmodel)
212C
213 IF(iadm==2) CALL hm_get_intv('NRadm',nradm,is_available,lsubmodel)
214C
215 IF(intth > 0) THEN
216 CALL hm_get_intv('Ithe_form',iform,is_available,lsubmodel)
217 CALL hm_get_intv('fct_ID_k',ikthe,is_available,lsubmodel)
218 CALL hm_get_intv('F_COND',fcond,is_available,lsubmodel)
219 ENDIF
220C
221C--------------------------------------------------
222C EXTRACT DATAS (REAL VALUES)
223C--------------------------------------------------
224 IF(igap==1.OR.igap==2)THEN
225 CALL hm_get_floatv('GAPSCALE',gapscale,is_available,lsubmodel,unitab)
226 CALL hm_get_floatv('GAPMAX',gapmax,is_available,lsubmodel,unitab)
227 CALL hm_get_floatv('DIST',depth,is_available,lsubmodel,unitab)
228 CALL hm_get_floatv('PMAX',pmax,is_available,lsubmodel,unitab)
229 ENDIF
230C
231 CALL hm_get_floatv('STMIN',stmin,is_available,lsubmodel,unitab)
232 CALL hm_get_floatv('STMAX',stmax,is_available,lsubmodel,unitab)
233 CALL hm_get_floatv('Pskid',pskid,is_available,lsubmodel,unitab)
234C
235 CALL hm_get_floatv('STFAC',stfac,is_available,lsubmodel,unitab)
236 CALL hm_get_floatv('FRIC',fric,is_available,lsubmodel,unitab)
237 CALL hm_get_floatv('MINI',gap,is_available,lsubmodel,unitab)
238 CALL hm_get_floatv('TSTART',startt,is_available,lsubmodel,unitab)
239 CALL hm_get_floatv('TSTOP',stopt,is_available,lsubmodel,unitab)
240C
241 CALL hm_get_floatv('VISC',visc,is_available,lsubmodel,unitab)
242 CALL hm_get_floatv('SORT_FACT',bumult,is_available,lsubmodel,unitab)
243C
244 CALL hm_get_floatv('Xfreq',alpha,is_available,lsubmodel,unitab)
245 CALL hm_get_floatv('scale1',xfricth,is_available,lsubmodel,unitab)
246C
247 IF (mfrot>0) THEN
248 CALL hm_get_floatv('C1',c1,is_available,lsubmodel,unitab)
249 CALL hm_get_floatv('C2',c2,is_available,lsubmodel,unitab)
250 CALL hm_get_floatv('C3',c3,is_available,lsubmodel,unitab)
251 CALL hm_get_floatv('C4',c4,is_available,lsubmodel,unitab)
252 CALL hm_get_floatv('C5',c5,is_available,lsubmodel,unitab)
253 ENDIF
254 IF (mfrot>1) THEN
255 CALL hm_get_floatv('C6',c6,is_available,lsubmodel,unitab)
256 ENDIF
257C
258 IF(iadm==2)THEN
259 CALL hm_get_floatv('Padm',padm,is_available,lsubmodel,unitab)
260 CALL hm_get_floatv('Angladm',angladm,is_available,lsubmodel,unitab)
261 ENDIF
262C
263C
264 IF(intth > 0) THEN
265 CALL hm_get_floatv('Kthe',kthe,is_available,lsubmodel,unitab)
266 CALL hm_get_floatv('Tint',tint,is_available,lsubmodel,unitab)
267 CALL hm_get_floatv('A_scale_k',xthe,is_available,lsubmodel,unitab)
268 CALL hm_get_floatv('f_rad',FRAD,IS_AVAILABLE,LSUBMODEL,UNITAB)
269 CALL HM_GET_FLOATV('d_rad',DRAD,IS_AVAILABLE,LSUBMODEL,UNITAB)
270 CALL HM_GET_FLOATV('heat_al',FHEAT,IS_AVAILABLE,LSUBMODEL,UNITAB)
271 CALL HM_GET_FLOATV('d_cond',DCOND,IS_AVAILABLE,LSUBMODEL,UNITAB)
272 ENDIF
273
274C------------------------------------------------------------
275C Card1 :flags
276C------------------------------------------------------------
277
278C
279C.....* CHECKS *........
280C
281 IF (IDEL21 < 0) THEN
282 IDELKEEP=1
283 IDEL21=ABS(IDEL21)
284 END IF
285C
286.OR. IF (IDEL21>2N2D==1) IDEL21 = 0
287C
288 IF(ILEV==0)THEN
289 ILEV=1
290 ELSEIF(ILEV==-1)THEN
291C emul past
292 ILEV=0
293 END IF
294C
295 IS1=1
296 INGR2USR => IGRSURF(1:NSURF)%ID
297 ISU1=NGR2USR(ISU1,INGR2USR,NSURF)
298 IS2=1
299 ISU2=NGR2USR(ISU2,INGR2USR,NSURF)
300C
301.AND. IF (IADM/=0NADMESH==0) THEN
302 CALL ANCMSG(MSGID=647,
303 . MSGTYPE=MSGWARNING,
304 . ANMODE=ANINFO_BLIND_2,
305 . I1=NOINT,
306 . C1=TITR)
307 IADM=0
308 ENDIF
309
310C.......* Storage IPARI FRIGAP *........
311 IPARI(21) = IGAP
312 IPARI(61) = IDELKEEP
313 IPARI(17) = IDEL21
314 IPARI(20) = ILEV
315 IPARI(34) = IGSTI
316 IPARI(51) = INVN ! flag to detect if normals are inverted (default =0)
317 IPARI(45)=ISU1
318 IPARI(46)=ISU2
319 IPARI(13)=IS1*10+IS2
320 IPARI(44)=IADM
321
322 IF(MULTIMP==0)THEN
323 MULTIMP=12
324 END IF
325 IPARI(23)=MULTIMP
326
327C------------------------------------------------------------
328C Card2
329C------------------------------------------------------------
330.OR. IF(IGAP==1IGAP==2)THEN
331C
332C.....* CHECKS *........
333C
334 IF(GAPSCALE==ZERO)GAPSCALE=ONE
335
336cE FLAG says if you have to calculate the nodal thickness in the engine
337 IF(IGAP==2) INTER_ITHKNOD=1 !defined in interface module (common_source directory)
338C
339 IF(IGAP==2)THEN
340 IF(PMAX==ZERO) PMAX=EP30
341 ELSE
342C Igap=1, Pmax is not used
343 PMAX=EP30
344 END IF
345C
346
347 ELSE
348C Igap=0, Pmax is not used
349 PMAX=EP30
350 ENDIF
351
352C.......* Storage IPARI FRIGAP *........
353 FRIGAP(19) = GAPSCALE
354 FRIGAP(16) = GAPMAX
355 FRIGAP(23) = DEPTH
356 FRIGAP(15) = PMAX
357 IPARI(52) = IFTLIM ! flag to deactivate tangential force limitation (default =1)
358
359C------------------------------------------------------------
360C Card3
361C------------------------------------------------------------
362
363C
364C.....* CHECKS *........
365C
366 IF(IGSTI==0)THEN
367 I7STIFS=1
368 IF(STMAX==ZERO)STMAX=EP30
369 FRIGAP(17) = STMIN
370 FRIGAP(18) = STMAX
371 ELSE
372 STMIN = ZERO
373 STMAX = EP30
374 END IF
375C
376 IF(PSKID==ZERO) PSKID=EP30
377 FRIGAP(35) = PSKID
378C------------------------------------------------------------
379C Card3
380C------------------------------------------------------------
381
382C
383C.....* CHECKS *........
384C
385 IF(STFAC==ZERO) THEN
386 STFAC=ONE
387 ENDIF
388 IF(IGSTI==1)STFAC=-STFAC
389C
390 IF (STOPT == ZERO) STOPT = EP30
391
392C.....* Storage IPARI FRIGAP *.......
393 IPARI(34)=IGSTI
394 FRIGAP(1)=FRIC
395 FRIGAP(2)=GAP
396 FRIGAP(3)=STARTT
397 FRIGAP(11)=STOPT
398
399C------------------------------------------------------------
400C Card4
401C------------------------------------------------------------
402 IF(VISC==ZERO) VISC=ONE
403
404C.....* Storage IPARI FRIGAP *.........
405
406 IPARI(11)=4*IBC1+2*IBC2+IBC3 + 8 *(4*IBC1M+2*IBC2M+IBC3M)
407 IPARI(22)=INACTI
408 FRIGAP(14)=VISC
409
410C BUMULT is increased for big models
411 IF(BUMULT==ZERO) THEN
412 BUMULT = BMUL0
413 IF(NUMNOD > 2500000) THEN
414 BUMULT = BMUL0*TWO
415 ELSEIF(NUMNOD > 1500000) THEN
416 BUMULT = BMUL0*THREE/TWO
417 END IF
418 END IF
419 FRIGAP(4)=BUMULT
420C------------------------------------------------------------
421C Card5 : FRICTION data
422C------------------------------------------------------------
423
424C
425C.....* CHECKS *........
426C
427 IF(FRIC == ZERO)THEN
428 IF(IFRICTH /= 0)THEN
429 FRIC = ONE
430 ENDIF
431 ENDIF
432
433 IF(XFRICTH == ZERO) XFRICTH=ONE
434C
435 IF (ALPHA==ZERO) IFQ = 0
436c
437 MODFR=2
438 IFQ = IFQ + 10
439c
440 IF (IFQ>0) THEN
441 IF (IFQ==10) XFILTR = ONE
442 IF (MOD(IFQ,10)==1) XFILTR = ALPHA
443 IF (MOD(IFQ,10)==2) XFILTR=FOUR*ATAN2(ONE,ZERO) / ALPHA
444 IF (MOD(IFQ,10)==3) XFILTR=FOUR*ATAN2(ONE,ZERO) * ALPHA
445 IF (XFILTR<ZERO) THEN
446 CALL ANCMSG(MSGID=554,
447 . MSGTYPE=MSGERROR,
448 . ANMODE=ANINFO_BLIND_1,
449 . I1=NOINT,
450 . C1=TITR,
451 . R1=ALPHA)
452.AND. ELSEIF (XFILTR>1MOD(IFQ,10)<=2) THEN
453 CALL ANCMSG(MSGID=554,
454 . MSGTYPE=MSGERROR,
455 . ANMODE=ANINFO_BLIND_1,
456 . I1=NOINT,
457 . C1=TITR,
458 . R1=ALPHA)
459 ENDIF
460 ELSE
461 XFILTR = ZERO
462 ENDIF
463C------------------------------------------------------------
464C Optional Card5 Card6 : C1...C5.C6 friction data
465C------------------------------------------------------------
466
467C.....* Storage IPARI FRIGAP *.........
468 FRIC_P(1) = C1
469 FRIC_P(2) = C2
470 FRIC_P(3) = C3
471 FRIC_P(4) = C4
472 FRIC_P(5) = C5
473 FRIC_P(6) = C6
474 IPARI(50) = IFRICTH
475 FRIGAP(34) = XFRICTH
476 IPARI(30) = MFROT
477 IPARI(31) = IFQ
478 IPARI(72) = INTFRIC
479
480C------------------------------------------------------------
481C Option Card8 :Adm = 2 adaptative mesh
482C------------------------------------------------------------
483C
484 IF(IADM==2)THEN
485 IF(NRADM==0) NRADM =3
486 IF(PADM==ZERO) PADM =ONE
487 ELSE
488 NRADM =1
489 PADM =ONE
490 ANGLADM=ZERO
491 END IF
492
493C.....* Storage IPARI FRIGAP *.........
494 CADM =COS(ANGLADM*PI/HUNDRED80)
495 KCONTACT =MAX(KCONTACT,IADM)
496 IPARI(49) =NRADM
497 FRIGAP(24)=PADM
498 FRIGAP(25)=CADM
499C------------------------------------------------------------
500C Option Card9 :Thermal input
501C------------------------------------------------------------
502 IF(INTTH > 0 ) THEN
503 INTHEAT = 1
504 ENDIF
505C
506C....* CHECKS *.............
507C
508 IF(KTHE == ZERO)THEN
509 IF(IKTHE /= 0)THEN
510 KTHE = ONE
511 CALL HM_GET_FLOATV_DIM('kthe' ,KTHE_UNIT ,IS_AVAILABLE, LSUBMODEL, UNITAB)
512 KTHE = ONE * KTHE_UNIT
513 ENDIF
514 ENDIF
515
516 IF (XTHE == ZERO) THEN
517 !units
518 CALL HM_GET_FLOATV_DIM('a_scale_k' ,XTHE_UNIT ,IS_AVAILABLE, LSUBMODEL, UNITAB)
519 XTHE = ONE * XTHE_UNIT
520 ENDIF
521
522 IF(IFORM /= 0)THEN
523 INTTH = 2
524 IPARI(47) = 2
525 ENDIF
526
527C.....* Storage IPARI FRIGAP *.........
528 IPARI(47) = INTTH
529 IPARI(48) = IFORM
530 IPARI(42) = IKTHE
531 FRIGAP(21 ) = TINT
532 FRIGAP(30) = XTHE
533C------------------------------------------------------------
534C Card5 : THERMAL MODELLING card2
535C------------------------------------------------------------
536
537C
538C....* CHECKS *.............
539C
540.AND. IF(ITHERM_FE == 0 INTTH > 0 ) THEN
541 INTHEAT = 0
542 IPARI(47) = 0
543 CALL ANCMSG(MSGID=702,
544 . MSGTYPE=MSGWARNING,
545 . ANMODE=ANINFO,
546 . I1=NOINT,
547 . C1=TITR,
548 . R1=DCOND,
549 . R2=DRAD)
550 ENDIF
551
552 IF(INTTH == 2 ) FTEMPVAR21 = 1
553
554 IF(FCOND ==0) DCOND = ZERO
555
556 IF(FRAD==ZERO ) DRAD = ZERO
557
558.AND..AND. IF(FCOND /= 0DCOND/=ZERODRAD==ZERO) THEN
559 DRAD = DCOND
560 CALL ANCMSG(MSGID=1810,
561 . MSGTYPE=MSGWARNING,
562 . ANMODE=ANINFO,
563 . I1=NOINT,
564 . C1=TITR,
565 . R1=DCOND,
566 . R2=DRAD)
567 ENDIF
568
569 IF(DCOND > DRAD) THEN
570 DCOND = DRAD
571 CALL ANCMSG(MSGID=1809,
572 . MSGTYPE=MSGWARNING,
573 . ANMODE=ANINFO,
574 . I1=NOINT,
575 . C1=TITR,
576 . R1=DCOND,
577 . R2=DRAD)
578 ENDIF
579
580C
581C.....* Storage IPARI FRIGAP *.........
582 FRIGAP(20) = KTHE
583 FRIGAP(31) = FRAD
584 FRIGAP(32) = DRAD
585 FRIGAP(33) = FHEAT
586 FRIGAP(36) = DCOND ! max conduction distance
587 IPARI(53) = FCOND ! function of variation of heat exchange as funct of distance
588C
589 IPARI(65) = INTKG
590C FRIGAP(10) is initialized but used only in engine for storing number of couples candidates
591 FRIGAP(10)=FLOAT(0)
592C
593C
594C------------------------------------------------------------
595C RENUMBERING OF FUNCTIONS AND SENSOR - USER TO INTERNAL ID
596C------------------------------------------------------------
597C
598 CALL INTER_DCOD_FUNCTION(NTYP,NI,IPARI,NPC1,NOM_OPT,NPC,TF)
599 CALL INTER_DCOD_SENSOR (NTYP,NI,IPARI,NOM_OPT,SENSORS)
600 CALL INTER_DCOD_FRICTION(NTYP,NI,IPARI,NOM_OPT,NOM_OPTFRIC,
601 . INTBUF_FRIC_TAB)
602
603C
604C------------------------------------------------------------
605C PRINTOUT
606C------------------------------------------------------------
607C
608 IF(IDSENS/=0) THEN
609 WRITE(IOUT,2101)IBC1,IBC2,IBC3,IBC1M,IBC2M,IBC3M,
610 . ILEV,IGSTI,STFAC,STMIN,STMAX,
611 . FRIC,IGAP,GAP,GAPMAX,GAPSCALE,IDSENS,
612 . BUMULT,INACTI,VISC,PMAX,MULTIMP,INVN,IFTLIM,
613 . PSKID
614 ELSE
615 WRITE(IOUT,2104)IBC1,IBC2,IBC3,IBC1M,IBC2M,IBC3M,
616 . ILEV,IGSTI,STFAC,STMIN,STMAX,
617 . FRIC,IGAP,GAP,GAPMAX,GAPSCALE,STARTT,STOPT,
618 . BUMULT,INACTI,VISC,PMAX,MULTIMP,INVN,IFTLIM,
619 . PSKID
620 ENDIF
621
622 IF(INTFRIC > 0 ) THEN
623 WRITE(IOUT,1527) INTFRIC
624 ELSE
625C
626 WRITE(IOUT,1520)MOD(IFQ,10), XFILTR
627C
628 IF(MFROT==0)THEN
629 IF(IFRICTH ==0 ) THEN
630 WRITE(IOUT,1524) FRIC
631 ELSE
632 WRITE(IOUT,1525) IFRICTH,XFRICTH,FRIC
633 ENDIF
634 ELSEIF(MFROT==1)THEN
635 WRITE(IOUT,1515)FRIC_P(1),FRIC_P(2),FRIC_P(3),
636 . FRIC_P(4),FRIC_P(5)
637 ELSEIF(MFROT==2)THEN
638 WRITE(IOUT,1522)FRIC,FRIC_P(1),FRIC_P(2),FRIC_P(3),
639 . FRIC_P(4),FRIC_P(5),FRIC_P(6)
640 ELSEIF(MFROT==3)THEN
641 WRITE(IOUT,1523)FRIC_P(1),FRIC_P(2),FRIC_P(3),
642 . FRIC_P(4),FRIC_P(5),FRIC_P(6)
643 ELSEIF(MFROT==4)THEN
644 WRITE(IOUT,1526) FRIC,FRIC_P(1),FRIC_P(2)
645 ENDIF
646 ENDIF
647C
648 IF(IDEL21/=0) THEN
649 WRITE(IOUT,'(a,i5/)')
650 . ' deletion flag on failure(1:yes) : ',IDEL21
651 IF(IDELKEEP == 1)THEN
652 WRITE(IOUT,'(a)')
653 . ' idel: DO not remove non-connected nodes from secondary surface'
654 ENDIF
655 ENDIF
656 IF(IADM/=0) THEN
657 WRITE(IOUT,*)' mesh refinement CASE of contact',
658 . ' (0:no/1:due to curvature/2:due to curvature or penetration)',
659 . ' set to ',IADM
660 IF(IADM==2)THEN
661 WRITE(IOUT,1557) NRADM,PADM,ANGLADM
662 END IF
663 ENDIF
664C
665 IF(INTTH > 0 )THEN
666 IF(IKTHE==0)THEN
667 WRITE(IOUT,2102) KTHE,TINT,FRAD,DRAD,FHEAT,IFORM,FCOND,DCOND
668 ELSE
669 WRITE(IOUT,2103) IKTHE,XTHE,KTHE,TINT,FRAD,DRAD,FHEAT,IFORM,FCOND,DCOND
670 END IF
671 END IF
672C--------------------------------------------------------------
673 IF(IS1==0)THEN
674 WRITE(IOUT,'(6x,a)')'no secondary surface input'
675 ELSEIF(IS1==1)THEN
676 WRITE(IOUT,'(6x,a)')'secondary surface input by segments'
677 ELSEIF(IS1==2)THEN
678 WRITE(IOUT,'(6x,a)')'secondary surface input by nodes'
679 ELSEIF(IS1==3)THEN
680 WRITE(IOUT,'(6x,a)')'secondary surface input by segments'
681 ELSEIF(IS1==4 )THEN
682 WRITE(IOUT,'(6x,a)')'secondary side input by bricks'
683 ELSEIF(IS1==5 )THEN
684 WRITE(IOUT,'(6x,a)')'secondary side input by solid elements'
685 ENDIF
686 IF(IS2==0)THEN
687 WRITE(IOUT,'(6x,a)')'no main surface input'
688 ELSEIF(IS2==1)THEN
689 WRITE(IOUT,'(6x,a)')'main surface input by segments'
690 ELSEIF(IS2==2)THEN
691 WRITE(IOUT,'(6x,a)')'main surface input by nodes'
692 ELSEIF(IS2==3)THEN
693 WRITE(IOUT,'(6x,a)')'main surface input by segments'
694 ELSEIF(IS2==4)THEN
695 WRITE(IOUT,'(6x,a)')'main surface refers ',
696 . 'to hyper-ellipsoidal surface'
697 ENDIF
698C
699C--------------------------------------------------------------
700C------------
701 RETURN
702 2502 FORMAT(' rupture parameters '
703 . /10X,'scal_f . . . . . . . . . . . . . . ',1PG20.13
704 . /10X,'scal_disp . . . . . . . . . . . . . ',1PG20.13
705 . /10X,'scal_sr . . . . . . . . . . . . . . ',1PG20.13
706 . /10X,'filtering coeff . . . . . . . . . . ',1PG20.13
707 . /10X,'default secondary area. . . . . . . . . ',1PG20.13
708 . /10X,'dn_max . . . . . . . . . . . . . . ',1PG20.13
709 . /10X,'dt_max . . . . . . . . . . . . . . ',1PG20.13
710 . /10X,'ifunn . . . . . . . . . . . . . . ',I10
711 . /10X,'ifunt . . . . . . . . . . . . . . ',I10
712 . /10X,'ifuns . . . . . . . . . . . . . . ',I10
713 . /10X,'imod . . . . . . . . . . . . . . ',I10
714 . /10X,'isym . . . . . . . . . . . . . . ',I10
715 . /10X,'ifiltr . . . . . . . . . . . . . . ',I10//)
716 1515 FORMAT(//
717 . ' friction model 1 (viscous polynomial)'/,
718 . ' mu = muo + c1 p + c2 v + c3 pv + c4 p^2 + c5 v^2'/,
719 . ' c1 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
720 . ' c2 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
721 . ' c3 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
722 . ' c4 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
723 . ' c5 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
724 . ' tangential pressure limit. . .. . . . . .',1PG20.13/)
725 1522 FORMAT(/
726 . ' friction model 2 (darmstad law) :'/,
727 . ' mu = muo+c1*exp(c2*v)*p^2+c3*exp(c4*v)*p+c5*exp(c6*v)'/,
728 . ' muo. . . . . . . . . . . . . . . . . . . ',1PG20.13/,
729 . ' c1 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
730 . ' c2 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
731 . ' c3 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
732 . ' c4 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
733 . ' c5 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
734 . ' c6 . . . . . . . . . . . . . . . . . . . ',1PG20.13/)
735 1523 FORMAT(/
736 . ' friction model 3 (renard law) :'/,
737 . ' c1 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
738 . ' c2 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
739 . ' c3 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
740 . ' c4 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
741 . ' c5 . . . . . . . . . . . . . . . . . . . ',1PG20.13/,
742 . ' c6 . . . . . . . . . . . . . . . . . . . ',1PG20.13/)
743 1524 FORMAT(/
744 . ' friction model 0 (coulomb law) :'/,
745 . ' friction coefficient . . . . . . . . . ',1PG20.13/)
746 1525 FORMAT(//
747 . ' friction model 0 (coulomb law) :'/,
748 . ' Function for friction coefficient wrt temperature',I10/,
749 . ' abscissa scale factor on ifuntcf. . . . . ',1PG20.13/,
750 . ' ordinate scale factor on ifuntcf . . . . ',1PG20.13/)
751 1526 FORMAT(/
752 . ' exponential decay friction law '/
753 . ' mu = c1+(MUo-c1)*exp(-c2*v)'/
754 . ' static coefficient muo . . . . . . . . . ',1PG20.13/,
755 . ' dynamic coefficient c1 . . . . . . . . . ',1PG20.13/,
756 . ' exponential decay coefficient c2 . . . . ',1PG20.13/)
757C
758 1527 FORMAT(/
759 . ' interface friction model. . . . . . . . . ',I10)
760
761
762 1520 FORMAT(
763 . ' friction filtering flag. . . . . . . . . ',I10/,
764 . ' filtering factor . . . . . . . . . . . . ',1PG20.13)
765 1557 FORMAT(
766 .' number of elements within a 90 degrees fillet ',I5/,
767 .' --------------------------------------------- '/,
768 .' criteria for refinement due to penetration : '/,
769 .' ------------------------------------------ '/,
770 .' minimum percentage of penetration ',
771 . 1PG20.13/,
772 .' maximum angle on main side at contact location ',
773 . 1PG20.13//)
774 2101 FORMAT(//
775 . ' type==21 parallel/stamping ' //,
776 . ' bound. cond. deleted after impact in x dir ',I1/,
777 . ' secondary node (1:YES 0:NO) y dir ',I1/,
778 . ' z dir ',I1/,
779 . ' bound. cond. deleted after impact in x dir ',I1/,
780 . ' main node (1:YES 0:NO) y dir ',I1/,
781 . ' z dir ',I1/,
782 . ' formulation level. . . . . . . . . . . . . .',I1/,
783 . ' (0:OLD, 1:OPTIMIZED). . . . . . . .',/,
784 . ' stiffness formulation. . . . . . . . . . ',I1/,
785 .' 0 : stiffness is computed from stiffness on secondary side'/,
786 .' 1 : stfac is a stiffness value '/,
787 . ' stiffness factor or stiffness value . . . . ',1PG20.13/,
788 . ' minimum stiffness. . . . . . . . . . . . . ',1PG20.13/,
789 . ' maximum stiffness. . . . . . . . . . . . . ',1PG20.13/,
790 . ' friction factor . . . . . . . . . . . . . . ',1PG20.13/,
791 . ' variable gap flag . . . . . . . . . . . . . ',I5/,
792 . ' minimum gap . . . . . . . . . . . . . . . . ',1PG20.13/,
793 . ' maximum gap (= 0. <=> NO MAXIMUM GAP) . . . ',1PG20.13/,
794 . ' gap scale factor. . . . . . . . . . . . . . ',1PG20.13/,
795 . ' start time/stop time activated by sensor id ',I10/,
796 . ' bucket factor . . . . . . . . . . . . . . . ',1PG20.13/,
797 . ' de-activation of initial penetrations . . . ',I10/,
798 . ' critical damping factor . . . . . . . . . . ',1PG20.13/,
799 . ' maximum pressure due to thickness variation '/,
800 . ' case of igap=2. . . . . . . ',1PG20.13/,
801 . ' mean possible number of impact/node . . . . ',I5/,
802 . ' inverted normals detection flag . . . . . . ',I5/,
803 . ' tangential force limitation flag . . . . . . ',I5/,
804 . ' maximum pressure for skid line output ',1PG20.13)
805 2104 FORMAT(//
806 . ' type==21 parallel/stamping ' //,
807 . ' bound. cond. deleted after impact in x dir ',I1/,
808 . ' secondary node (1:YES 0:NO) y dir ',I1/,
809 . ' z dir ',I1/,
810 . ' bound. cond. deleted after impact in x dir ',I1/,
811 . ' main node (1:YES 0:NO) y dir ',I1/,
812 . ' z dir ',I1/,
813 . ' formulation level. . . . . . . . . . . . . .',I1/,
814 . ' (0:OLD, 1:OPTIMIZED). . . . . . . .',/,
815 . ' stiffness formulation. . . . . . . . . . ',I1/,
816 .' 0 : stiffness is computed from stiffness on secondary side'/,
817 .' 1 : stfac is a stiffness value '/,
818 . ' stiffness factor or stiffness value . . . ',1PG20.13/,
819 . ' minimum stiffness. . . . . . . . . . . . ',1PG20.13/,
820 . ' maximum stiffness. . . . . . . . . . . . ',1PG20.13/,
821 . ' friction factor . . . . . . . . . . . . . ',1PG20.13/,
822 . ' variable gap flag . . . . . . . . . . . . ',I5/,
823 . ' minimum gap . . . . . . . . . . . . . . . ',1PG20.13/,
824 . ' maximum gap (= 0. <=> NO MAXIMUM GAP) . . ',1PG20.13/,
825 . ' gap scale factor. . . . . . . . . . . . . ',1PG20.13/,
826 . ' start time. . . . . . . . . . . . . . . . ',1PG20.13/,
827 . ' stop time . . . . . . . . . . . . . . . . ',1PG20.13/,
828 . ' bucket factor . . . . . . . . . . . . . . ',1PG20.13/,
829 . ' de-activation of initial penetrations . . ',I10/,
830 . ' critical damping factor . . . . . . . . . ',1PG20.13/,
831 . ' maximum pressure due to thickness variation'/,
832 . ' case of igap=2. . . . . . ',1PG20.13/,
833 . ' mean possible number of impact/node . . . ',I5/,
834 . ' inverted normals detection flag . . . . . . ',I5/,
835 . ' tangential force limitation flag . . . . . . ',I5/,
836 . ' maximum pressure for skid line output ',1PG20.13)
837 2102 FORMAT(//
838 . ' thermal interface ' //,
839 . ' thermal heat exchange coefficient . . . . .',1PG20.13/,
840 . ' interface temperature . . . . . . . . . . ',1PG20.13/,
841 . ' radiation factor . . . . . . . . . . . . . ',1PG20.13/,
842 . ' maximum distance for radiation computation.',1PG20.13/,
843 . ' frictional heat transfer. . . . . . . . . .',1PG20.13/,
844 . ' formulation choice : . . . . . . . . . . . ',I10,/,
845 . ' 0 : heat transfer between secondary side',/,
846 . ' and constant temperature in interface',/,
847 . ' 1 : heat exchange between pieces in contact'/
848 . ' Function for thermal heat exchange coefficient wrt distance',i10/,
849 . ' Maximum distance for conductive heat exchange',1pg20.13)
850 2103 FORMAT(//
851 . ' Thermal interface ' //,
852 . ' Function for thermal heat exchange coefficient wrt contact pressure',i10/,
853 . ' Abscissa scale factor on IFUNTCK. . . . . ',1pg20.13/,
854 . ' Ordinate scale factor on IFUNTCK . . . . ',1pg20.13/,
855 . ' Interface temperature . . . . . . . . . . ',1pg20.13/,
856 . ' Radiation factor . . . . . . . . . . . . . ',1pg20.13/,
857 . ' Maximum distance for radiation computation.',1pg20.13/,
858 . ' Frictional heat generation . . . . . . . . ',1pg20.13/,
859 . ' FORMULATION CHOICE : . . . . . . . . . . . ',i10,/,
860 . ' 0 : heat transfer between SECONDARY side',/,
861 . ' and constant temperature in interface',/,
862 . ' 1 : heat Exchange between pieces in contact'/
863 . ' Function for thermal heat exchange coefficient wrt distance',i10/,
864 . ' Maximum distance for conductive heat exchange',1pg20.13)
865C
866 END
my_real
#define my_real
Definition cppsort.cpp:32
damping
subroutine damping(nodft, nodlt, v, vr, a, ar, damp, ms, in, igrnod, dim, itask, weight, tagslv_rby, wfext)
Definition damping.F:882
alpha
#define alpha
Definition eval.h:35
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528
hm_read_inter_type21
subroutine hm_read_inter_type21(ipari, stfac, frigap, noint, ni, igrsurf, xfiltr, fric_p, npc1, sensors, nom_opt, unitab, lsubmodel, titr, npc, tf, npari, nparir, snpc, snpc1, lnopt1, itherm_fe, intheat, nom_optfric, intbuf_fric_tab)
Definition hm_read_inter_type21.F:46
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
id
initmumps id
Definition multiplerhs_example.m:4
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
groupdef_mod
Definition groupdef_mod.F:668
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
unitab_mod
Definition unitab_mod.F:517
radiation
subroutine radiation(ibcr, fradia, npc, tf, x, temp, nsensor, sensor_tab, fthe, iad, fthesky, python, glob_therm)
Definition radiation.F:38
main
int main(int argc, char *argv[])
Definition CMakeCCompilerId.c:841
static
subroutine static(v, vr, a, ar, ms, in, igrnod, weight_md, wfext)
Definition static.F:33
groupdef_mod::surf_
Definition groupdef_mod.F:751
submodel_mod::submodel_data
Definition submodel_mod.F:728