OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i24buce.F File Reference
#include "implicit_f.inc"
#include "comlock.inc"
#include "com01_c.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine i24buce (x, v, irect, nsv, stfn, nmn, nrtm, nsn, cand_e, cand_n, gap, noint, ii_stok, ncontact, bminma, marge, curv_max, pmax_gap, vmaxdt, nb_n_b, eshift, ild, nin, stf, gap_s, nsnr, ncont, gap_m, itask, bgapsmx, i_mem, pene_old, itab, nbinflg, mbinflg, ilev, msegtyp, edge_l2, iedge, iseadd, isedge, cand_t, flagremnode, kremnod, remnod, cand_a, renum, nsnrold, irtse, is2se, nsne, dgapload, intheat, idt_therm, nodadt_therm)

Function/Subroutine Documentation

◆ i24buce()

subroutine i24buce ( x,
v,
integer, dimension(4,*) irect,
integer, dimension(*) nsv,
stfn,
integer nmn,
integer nrtm,
integer nsn,
integer, dimension(*) cand_e,
integer, dimension(*) cand_n,
gap,
integer noint,
integer ii_stok,
integer ncontact,
bminma,
marge,
curv_max,
pmax_gap,
vmaxdt,
integer nb_n_b,
integer eshift,
integer ild,
integer nin,
stf,
gap_s,
integer nsnr,
integer ncont,
gap_m,
integer itask,
bgapsmx,
integer i_mem,
pene_old,
integer, dimension(*) itab,
integer, dimension(*) nbinflg,
integer, dimension(*) mbinflg,
integer ilev,
integer, dimension(*) msegtyp,
edge_l2,
integer iedge,
integer, dimension(*) iseadd,
integer, dimension(*) isedge,
integer, dimension(*) cand_t,
integer flagremnode,
integer, dimension(*) kremnod,
integer, dimension(*) remnod,
integer, dimension(*) cand_a,
integer, dimension(*) renum,
integer nsnrold,
integer, dimension(5,*) irtse,
integer, dimension(2,*) is2se,
integer nsne,
intent(in) dgapload,
integer, intent(in) intheat,
integer, intent(in) idt_therm,
integer, intent(in) nodadt_therm )

Definition at line 33 of file i24buce.F.

45C============================================================================
46C M o d u l e s
47C-----------------------------------------------
48 USE tri7box
49 USE message_mod
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54#include "comlock.inc"
55C-----------------------------------------------
56C C o m m o n B l o c k s
57C-----------------------------------------------
58#include "com01_c.inc"
59#include "task_c.inc"
60C-----------------------------------------------
61C D u m m y A r g u m e n t s
62C-----------------------------------------------
63 INTEGER NMN, NRTM, NSN, NOINT,IDT,NIN, NSNR
64 INTEGER IRECT(4,*),NSV(*),ITASK,ILEV,IEDGE,NSNE
65 INTEGER CAND_E(*),CAND_N(*),MSEGTYP(*)
66 INTEGER NCONTACT,ESHIFT,ILD,NB_N_B, NCONT,I_MEM,
67 . II_STOK,ITAB(*),NBINFLG(*),MBINFLG(*),ISEADD(*),
68 . ISEDGE(*),CAND_T(*),FLAGREMNODE,KREMNOD(*),REMNOD(*),
69 . CAND_A(*) ,RENUM(*),NSNROLD,IRTSE(5,*),IS2SE(2,*)
70C REAL
71 my_real
72 . gap,pmax_gap,vmaxdt,
73 . bminma(6),curv_max(nrtm),bgapsmx, marge
74 my_real , INTENT(IN) :: dgapload
75 my_real
76 . x(3,*), v(3,*), stfn(*),pene_old(5,nsn),
77 . stf(*), gap_s(*), gap_m(*),edge_l2(*)
78C-----------------------------------------------
79C L o c a l V a r i a b l e s
80C-----------------------------------------------
81 INTEGER I, J, IP0, IP1, MAXSIZ,
82 . LOC_PROC, N, ISZNSNR,
83 . NSNFIOLD(NSPMD)
84 INTEGER, INTENT(IN) :: INTHEAT
85 INTEGER, INTENT(IN) :: IDT_THERM
86 INTEGER, INTENT(IN) :: NODADT_THERM
87C REAL
88 my_real
89 . xyzm(6), aaa
90 INTEGER NBX,NBY,NBZ,NSNF,NSNL
91 INTEGER (KIND=8) :: NBX8,NBY8,NBZ8,RES8,LVOXEL8
92
93C-----------------------------------------------
94C S o u r c e L i n e s
95C-----------------------------------------------
96c save max penetration (pene-gap) at sorting time
97c in PENE_OLD(4,i) used to compute
98c DELTA_PMAX_GAP_L = MAX(DELTA_PMAX_GAP_L,PENE_OLD(3,I)-PENE_OLD(4,I))
99c in I24BUCE_CRIT at each cycle
100C-----------------------------------------------
101 nsnf = 1 + itask*nsn / nthread
102 nsnl = (itask+1)*nsn / nthread
103 pene_old(4,nsnf:nsnl) = pene_old(3,nsnf:nsnl)
104C-----------------------------------------------
105C
106C----- TRI DES ELEMENTS ET DES NOEUDS
107C
108C-----------------------------------------------
109C
110C----- BORNES DU DOMAINE DEJA CALCULEES
111C
112 xyzm(1) = bminma(4)
113 xyzm(2) = bminma(5)
114 xyzm(3) = bminma(6)
115 xyzm(4) = bminma(1)
116 xyzm(5) = bminma(2)
117 xyzm(6) = bminma(3)
118 i_mem = 0
119C
120 isznsnr = nsnr
121C=============================================================================
122c
123c Si MARGE - DELTA_PMAX_GAP - SOMME(Vrel*dt) < ZERO => RETRI
124c DELTA_PMAX_GAP = MAX((PENEmax(i)-GAP)-(PENEmax_Tri(i)-GAP)
125c
126C=============================================================================
127
128 IF( nmn /= 0) THEN
129 aaa = sqrt(nmn /
130 . ((bminma(1)-bminma(4))*(bminma(2)-bminma(5))
131 . +(bminma(2)-bminma(5))*(bminma(3)-bminma(6))
132 . +(bminma(3)-bminma(6))*(bminma(1)-bminma(4))))
133 ELSE
134 aaa = 0
135 ENDIF
136
137 aaa = 0.75*aaa
138
139 nbx = nint(aaa*(bminma(1)-bminma(4)))
140 nby = nint(aaa*(bminma(2)-bminma(5)))
141 nbz = nint(aaa*(bminma(3)-bminma(6)))
142 nbx = max(nbx,1)
143 nby = max(nby,1)
144 nbz = max(nbz,1)
145
146 nbx8=nbx
147 nby8=nby
148 nbz8=nbz
149 res8=(nbx8+2)*(nby8+2)*(nbz8+2)
150 lvoxel8 = lvoxel
151
152 IF(res8 > lvoxel8) THEN
153 aaa = lvoxel
154 aaa = aaa/((nbx8+2)*(nby8+2)*(nbz8+2))
155 aaa = aaa**(third)
156 nbx = int((nbx+2)*aaa)-2
157 nby = int((nby+2)*aaa)-2
158 nbz = int((nbz+2)*aaa)-2
159 nbx = max(nbx,1)
160 nby = max(nby,1)
161 nbz = max(nbz,1)
162 ENDIF
163
164 nbx8=nbx
165 nby8=nby
166 nbz8=nbz
167 res8=(nbx8+2)*(nby8+2)*(nbz8+2)
168
169 IF(res8 > lvoxel8) THEN
170 nbx = min(100,max(nbx8,1))
171 nby = min(100,max(nby8,1))
172 nbz = min(100,max(nbz8,1))
173 ENDIF
174c IF(MAX(NBX,NBY,NBZ) > 100) stop 678
175
176C initialisation complete de VOXEL
177C (en // SMP il y a possibilite de redondance de traitement mais no pb)
178 DO i=inivoxel,(nbx+2)*(nby+2)*(nbz+2)
179 voxel1(i)=0
180 ENDDO
181 inivoxel = max(inivoxel,(nbx+2)*(nby+2)*(nbz+2)+1)
182
183c write(iout,*)'i24buce 2'
184 CALL i24trivox(
185 1 nsn ,nsnr ,isznsnr ,i_mem ,vmaxdt ,
186 2 irect ,x ,stf ,stfn ,xyzm ,
187 3 nsv ,ii_stok ,cand_n ,eshift ,cand_e ,
188 4 ncontact,noint ,v ,bgapsmx ,
189 5 voxel1 ,nbx ,nby ,nbz ,pmax_gap ,
190 6 nrtm ,gap_s ,gap_m ,marge ,curv_max ,
191 7 nin ,itask ,pene_old,itab ,nbinflg ,
192 8 mbinflg ,ilev ,msegtyp ,edge_l2 ,iedge ,
193 9 iseadd ,isedge ,cand_t ,flagremnode,kremnod,
194 a remnod ,cand_a ,renum ,nsnrold ,irtse ,
195 b is2se ,nsne ,dgapload,intheat,idt_therm,nodadt_therm)
196 234 continue
197c write(iout,*)'i24buce 3'
198
199C
200C I_MEM = 2 ==> PAS ASSEZ DE MEMOIRE CANDIDATS
201
202c write(iout,*)'i24buce 4'
203C
204 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
i24trivox
subroutine i24trivox(nsn, nsnr, isznsnr, i_mem, vmaxdt, irect, x, stf, stfn, xyzm, nsv, ii_stok, cand_n, eshift, cand_e, mulnsn, noint, v, bgapsmx, voxel, nbx, nby, nbz, pmax_gap, nrtm, gap_s, gap_m, marge, curv_max, nin, itask, pene_old, itab, nbinflg, mbinflg, ilev, msegtyp, edge_l2, iedge, iseadd, isedge, cand_t, flagremnode, kremnod, remnod, cand_a, renum, nsnrold, irtse, is2se, nsne, dgapload, intheat, idt_therm, nodadt_therm)
Definition i24trivox.F:47
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
message_mod
Definition message_mod.F:1249
tri7box::voxel1
integer, dimension(lvoxel) voxel1
Definition tri7box.F:53
tri7box::inivoxel
integer inivoxel
Definition tri7box.F:53
tri7box::lvoxel
integer lvoxel
Definition tri7box.F:51