OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i24buce.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!|| i24buce ../engine/source/interfaces/intsort/i24buce.F
25!||--- called by ------------------------------------------------------
26!|| i24main_tri ../engine/source/interfaces/intsort/i24main_tri.F
27!||--- calls -----------------------------------------------------
28!|| i24trivox ../engine/source/interfaces/intsort/i24trivox.f
29!||--- uses -----------------------------------------------------
30!|| message_mod ../engine/share/message_module/message_mod.F
31!|| tri7box ../engine/share/modules/tri7box.F
32!||====================================================================
33 SUBROUTINE i24buce(
34 1 X ,V ,IRECT ,NSV ,STFN ,
35 2 NMN ,NRTM ,NSN ,CAND_E ,CAND_N ,
36 3 GAP ,NOINT ,II_STOK ,NCONTACT ,BMINMA ,
37 4 MARGE ,CURV_MAX,PMAX_GAP,VMAXDT ,NB_N_B ,
38 5 ESHIFT ,ILD ,NIN ,STF ,GAP_S ,
39 6 NSNR ,NCONT ,GAP_M ,ITASK ,BGAPSMX ,
40 7 I_MEM ,PENE_OLD,ITAB ,NBINFLG ,MBINFLG ,
41 8 ILEV ,MSEGTYP ,EDGE_L2 ,IEDGE ,ISEADD ,
42 9 ISEDGE ,CAND_T ,FLAGREMNODE,KREMNOD,REMNOD ,
43 A CAND_A ,RENUM ,NSNROLD ,IRTSE ,IS2SE ,
44 B NSNE ,DGAPLOAD,INTHEAT,IDT_THERM,NODADT_THERM)
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
205 END
i24buce
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)
Definition i24buce.F:45
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