OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fvmesh0.F File Reference
#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "units_c.inc"
#include "scr17_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine fvmesh0 (t_monvol, xyzini, ixs, ixc, ixtg, pm, ipm, igrsurf, xyzref, nb_node)
subroutine applysort2fvm (t_monvol)

Function/Subroutine Documentation

◆ applysort2fvm()

subroutine applysort2fvm ( type(monvol_struct_), dimension(nvolu), intent(inout) t_monvol)

Definition at line 326 of file fvmesh0.F.

327C Description: Apply renumbering to FVMBAGS (ELTG et TBA)
328C after S*HEAD S*TAILS routines
329C-----------------------------------------------
330C M o d u l e s
331C-----------------------------------------------
332 USE fvbag_mod
333 USE fvmbag_meshcontrol_mod
334 USE reorder_mod
335 USE monvol_struct_mod
336C-----------------------------------------------
337C I m p l i c i t T y p e s
338C-----------------------------------------------
339#include "implicit_f.inc"
340C-----------------------------------------------
341C C o m m o n B l o c k s
342C-----------------------------------------------
343#include "com04_c.inc"
344#include "scr17_c.inc"
345C-----------------------------------------------
346C D u m m y A r g u m e n t s
347C-----------------------------------------------
348 TYPE(MONVOL_STRUCT_), DIMENSION(NVOLU), INTENT(INOUT) :: T_MONVOL
349C-----------------------------------------------
350C L o c a l V a r i a b l e s
351C-----------------------------------------------
352 INTEGER :: IFV, N, ITYP, NNS, NTG, NBA, NTGA, NNA, NNI, NTGI, NNT, IREF, ID, NTGT
353 INTEGER J, IEL
354 INTEGER, DIMENSION(:), ALLOCATABLE :: TMP
355C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
356 ifv=0
357 DO n = 1, nvolu
358 ityp = t_monvol(n)%TYPE
359 id = t_monvol(n)%ID
360 IF (ityp == 6 .OR. ityp == 8 .OR. ityp == 11) THEN
361 iref = t_monvol(n)%IVOLU(59)
362 ifv = ifv + 1
363 nns = t_monvol(n)%NNS
364 ntg = t_monvol(n)%NTG
365 nni = t_monvol(n)%NNI
366 ntgi = t_monvol(n)%NTGI
367 nba = t_monvol(n)%NBRIC
368 ntga = t_monvol(n)%NTGA
369 nna = t_monvol(n)%NNA
370 nnt = nns + nni
371 ntgt = ntg + ntgi
372
373C----------------------------------------------
374C Parameters for automatic volume splitting
375C----------------------------------------------
376C PERMUTATION ELTG
377 IF (ntg > 0) THEN
378 ALLOCATE(tmp(ntg))
379 DO j = 1, ntg
380 iel = t_monvol(n)%ELTG(j)
381 IF (iel <= numelc) THEN
382 tmp(j) = permutation%SHELL(iel + numelc)
383 ELSE IF (iel > numelc) THEN
384 tmp(j) = numelc + permutation%TRIANGLE(iel - numelc + numeltg)
385 ENDIF
386 ENDDO
387 t_monvol(n)%ELTG(1:ntg) = tmp(1:ntg)
388 DEALLOCATE(tmp)
389 ENDIF
390
391 IF (ntgi > 0) THEN
392 ALLOCATE(tmp(ntgi))
393 DO j = 1, ntgi
394 iel = t_monvol(n)%ELTG(ntg + j)
395 IF (iel <= numelc) THEN
396 tmp(j) = permutation%SHELL(iel + numelc)
397 ELSE IF (iel > numelc) THEN
398 tmp(j) = numelc + permutation%TRIANGLE(iel - numelc + numeltg)
399 ENDIF
400 ENDDO
401 t_monvol(n)%ELTG(ntg + 1: ntg + ntgi) = tmp(1:ntgi)
402 DEALLOCATE(tmp)
403 ENDIF
404
405C PERMUTATION TBA
406 IF (ntg > 0 .AND. (.NOT. tetramesher_used)) THEN
407 DO j = 1, nba
408 iel = t_monvol(n)%TBRIC(1, j)
409 t_monvol(n)%TBRIC(1, j) = permutation%SOLID(iel + numels)
410 ENDDO
411 ELSE IF(ntg > 0 .AND. tetramesher_used) THEN
412C Not possible
413 ENDIF
414 ENDIF
415 ENDDO
416C
417
418 RETURN
419C
id
initmumps id
Definition multiplerhs_example.m:4
fvbag_mod
Definition fvbag_mod.F:44
fvmbag_meshcontrol_mod
Definition fvmbag_meshcontrol_mod.F:46
fvmbag_meshcontrol_mod::tetramesher_used
logical tetramesher_used
Definition fvmbag_meshcontrol_mod.F:49
monvol_struct_mod
Definition monvol_struct_mod.F:66
reorder_mod
Definition reorder_mod.F:45
reorder_mod::permutation
type(reorder_struct_) permutation
Definition reorder_mod.F:54

◆ fvmesh0()

subroutine fvmesh0 ( type(monvol_struct_), dimension(nvolu), intent(inout) t_monvol,
xyzini,
integer, dimension(nixs,*) ixs,
integer, dimension(nixc,*) ixc,
integer, dimension(nixtg,*) ixtg,
pm,
integer, dimension(npropmi,*) ipm,
type (surf_), dimension(nsurf) igrsurf,
xyzref,
integer nb_node )

Definition at line 54 of file fvmesh0.F.

55C-----------------------------------------------
56C M o d u l e s
57C-----------------------------------------------
58 USE my_alloc_mod
59 USE fvbag_mod
60 USE fvmbag_meshcontrol_mod
61 USE reorder_mod
62 USE message_mod
63 USE groupdef_mod
64 USE monvol_struct_mod
65 USE names_and_titles_mod , ONLY : nchartitle
66 USE fvelinte_mod , ONLY : fvelinte
67 USE fvvolu_mod , ONLY : fvvolu
68 USE fvinjnormal_mod , ONLY : fvinjnormal
69 USE fvlength_mod , ONLY : fvlength
70 use element_mod , only : nixs,nixc,nixtg
71C-----------------------------------------------
72C I m p l i c i t T y p e s
73C-----------------------------------------------
74#include "implicit_f.inc"
75C-----------------------------------------------
76C C o m m o n B l o c k s
77C-----------------------------------------------
78#include "com01_c.inc"
79#include "com04_c.inc"
80#include "param_c.inc"
81#include "units_c.inc"
82#include "scr17_c.inc"
83C-----------------------------------------------
84C D u m m y A r g u m e n t s
85C-----------------------------------------------
86 TYPE(MONVOL_STRUCT_), DIMENSION(NVOLU), INTENT(INOUT) :: T_MONVOL
87 INTEGER IXS(NIXS,*), IXC(NIXC,*), IXTG(NIXTG,*)
88 INTEGER IPM(NPROPMI,*)
89 INTEGER NB_NODE
90 my_real xyzini(3,nb_node), pm(npropm,*), xyzref(3,nb_node)
91 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
92C-----------------------------------------------
93C L o c a l V a r i a b l e s
94C-----------------------------------------------
95 INTEGER IFV,
96 . N, ITYP, NNS, NTG, NBRIC, NBX, NBY, NNB,
97 . NBA, NTGA, NNA,
98 . PFLAG, NNI, NTGI, ILVOUT, NNFV, NNT, NSURFI,
99 . NSEG, IREF, NTRFV, NPOLH, ID,
100 . NTGT, INODE
101 INTEGER I
102 INTEGER IBID
103 CHARACTER(len=nchartitle) :: TITR
104 INTEGER, DIMENSION(:), ALLOCATABLE :: MINUS_SIGN_REVERSE
105 my_real, DIMENSION(:,:), ALLOCATABLE :: x
106 my_real dirx,diry,dirz,dir2x,dir2y,dir2z,origx,origy,origz,lx,ly,lz
107C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
108 CALL my_alloc(x,3,nb_node)
109 pflag = 0
110 ifv=0
111 DO n = 1, nvolu
112 ityp = t_monvol(n)%TYPE
113 id = t_monvol(n)%ID
114 titr = t_monvol(n)%TITLE
115 IF (ityp == 6 .OR. ityp == 8) THEN
116 IF (pflag == 0) THEN
117 WRITE(iout,1000)
118 pflag = 1
119 ENDIF
120C
121 iref = t_monvol(n)%IVOLU(59)
122 IF(iref==0) THEN
123 x=xyzini
124 ELSE
125 x=xyzref
126 ENDIF
127C
128 ifv=ifv+1
129 nns = t_monvol(n)%NNS
130 ntg = t_monvol(n)%NTG
131 nni = t_monvol(n)%NNI
132 ntgi = t_monvol(n)%NTGI
133 nba = t_monvol(n)%NBRIC
134 ntga = t_monvol(n)%NTGA
135 nna = t_monvol(n)%NNA
136 nnt = nns + nni
137 ntgt = ntg + ntgi
138
139 t_monvol(n)%KR5 = 1+nrvolu*nvolu+lrcbag+lrbagjet+lrbaghol+t_monvol(n)%IVOLU(34) + nnt*6
140C----------------------------------------------
141C Automatic volume cutting parameters
142C----------------------------------------------
143 dirx = t_monvol(n)%RVOLU(35)
144 diry = t_monvol(n)%RVOLU(36)
145 dirz = t_monvol(n)%RVOLU(37)
146 dir2x = t_monvol(n)%RVOLU(38)
147 dir2y = t_monvol(n)%RVOLU(39)
148 dir2z = t_monvol(n)%RVOLU(40)
149 origx = t_monvol(n)%RVOLU(41)
150 origy = t_monvol(n)%RVOLU(42)
151 origz = t_monvol(n)%RVOLU(43)
152 lx = t_monvol(n)%RVOLU(44)
153 ly = t_monvol(n)%RVOLU(45)
154 lz = t_monvol(n)%RVOLU(53)
155 fvdata(ifv)%L_TYPE = 0
156 fvdata(ifv)%ID_DT_OPTION = t_monvol(n)%IVOLU(27)
157 fvdata(ifv)%LAMBDA = zero
158 fvdata(ifv)%DTOLD = zero
159 fvdata(ifv)%CFL_COEF = t_monvol(n)%RVOLU(71)
160 fvdata(ifv)%DTMIN = t_monvol(n)%RVOLU(72)
161 fvdata(ifv)%PDISP = zero
162 fvdata(ifv)%PDISP_OLD = zero
163C
164 CALL fvverif(
165 . ntga, t_monvol(n)%ELEMA, x, id,
166 . dirx, diry, dirz, dir2x, dir2y,
167 . dir2z, origx, origy, origz,
168 . lx, ly, lz, t_monvol(n)%NODES, t_monvol(n)%IBUFA, t_monvol(n)%TAGELA,
169 . titr)
170C
171 t_monvol(n)%RVOLU(44) = lx
172 t_monvol(n)%RVOLU(45) = ly
173 t_monvol(n)%RVOLU(53) = lz
174C
175 nbx = t_monvol(n)%IVOLU(54)
176 nby = t_monvol(n)%IVOLU(55)
177 nbric = nbx * nby
178 nnb = (nbx + 1) * (nby + 1) * 2
179 ALLOCATE(fvdata(ifv)%BRIC(8,nbric),
180 . fvdata(ifv)%TBRIC(13,nbric),
181 . fvdata(ifv)%XB(3,nnb),
182 . fvdata(ifv)%SFAC(6,4,nbric))
183C
184 CALL fvbric(t_monvol(n)%IVOLU, t_monvol(n)%RVOLU, t_monvol(n)%NODES, x, nns)
185C
186 t_monvol(n)%IVOLU(50) = t_monvol(n)%IVOLU(46)
187 t_monvol(n)%IVOLU(51) = t_monvol(n)%IVOLU(47)
188 t_monvol(n)%IVOLU(52) = t_monvol(n)%IVOLU(48)
189 t_monvol(n)%IVOLU(53) = t_monvol(n)%IVOLU(49)
190C
191 CALL fvmesh1(
192 . t_monvol(n)%NODES, t_monvol(n)%ELEM, x, t_monvol(n)%IVOLU, fvdata(ifv)%BRIC,
193 . fvdata(ifv)%XB, t_monvol(n)%RVOLU, ntg, ntgi, nbric, fvdata(ifv)%TBRIC,
194 . fvdata(ifv)%SFAC, fvdata(ifv)%DLH, nba, ntga,
195 . t_monvol(n)%TBRIC, t_monvol(n)%TFAC, t_monvol(n)%TAGELS, t_monvol(n)%IBUFA,
196 . t_monvol(n)%ELEMA, t_monvol(n)%TAGELA, ixs, id ,titr, nb_node, ityp)
197
198 IF (kmesh(n) >= 2) THEN
199 t_monvol(n)%KRA5 = 1 + nrvolu * nvolu + lrcbag + lrbagjet + lrbaghol +
200 . t_monvol(n)%IVOLU(34) + 7*nnt+4*ntgt+6*nna
201 t_monvol(n)%VELOCITY(1:3, 1:nna) = zero
202 DO i = 1, nna
203 inode = t_monvol(n)%IBUFA(i)
204 t_monvol(n)%NODE_COORD(1, i) = node_coord(1, inode)
205 t_monvol(n)%NODE_COORD(2, i) = node_coord(2, inode)
206 t_monvol(n)%NODE_COORD(3, i) = node_coord(3, inode)
207 ENDDO
208 ENDIF
209C
210 IF (ntgi > 0) THEN
211 ilvout = t_monvol(n)%IVOLU(44)
212 nsurfi=t_monvol(n)%INT_SURFID
213 nseg=igrsurf(nsurfi)%NSEG
214C-----------------------------
215C Internal surface elements
216C-----------------------------
217 t_monvol(n)%POROSITY(1:ntgi) = zero
218 CALL fvelinte(t_monvol(n)%NODES, t_monvol(n)%ELEM(1, ntg + 1), ixc, ixtg,
219 . pm, ipm, ilvout, ifv, nnt, ntg, t_monvol(n)%POROSITY,
220 . nseg,igrsurf(nsurfi)%ELTYP, ntgi, t_monvol(n)%ELTG,
221 . nb_node,igrsurf(nsurfi)%ELEM)
222C--------------------------------
223C Injectors on internal surface
224C--------------------------------
225 ALLOCATE(minus_sign_reverse(ntgi))
226 minus_sign_reverse(:) = 0
227 CALL fvinjectint(
228 . t_monvol(n)%NODES, t_monvol(n)%ELEM(1, ntg + 1), t_monvol(n)%IBAGJET,
229 . t_monvol(n)%NJET , igrsurf ,
230 . t_monvol(n)%ITAGEL(ntg + 1), nns+nni , ntgi,nb_node,
231 . minus_sign_reverse)
232C--------------------------------
233C Porous surface on internal surface
234C--------------------------------
235 CALL fvventholeint(
236 . t_monvol(n)%NODES, t_monvol(n)%ELEM(1, ntg + 1), t_monvol(n)%IBAGHOL,
237 . t_monvol(n)%NVENT, igrsurf ,
238 . t_monvol(n)%ITAGEL(ntg + 1), nns+nni , ntgi, nb_node)
239C
240 CALL fvelsurf(
241 . t_monvol(n)%NODES, t_monvol(n)%ELEM(1, ntg + 1), ibid, ixc, ixtg, ntgi,
242 . t_monvol(n)%ELTG(ntg + 1), t_monvol(n)%MATTG(ntg + 1), nb_node, .false.)
243C-----------------------------------------------------------------------
244C REDEFINE INTERNAL TRIANGLE FOR NORMAL CONSISTENCY FOR INJECTORS
245C-----------------------------------------------------------------------
246 CALL fvinjnormal(
247 . t_monvol(n)%NODES, t_monvol(n)%ELEM(1, ntg + 1), ixc, ixtg,
248 . t_monvol(n)%ELTG(ntg + 1), ntgi, ilvout,
249 . minus_sign_reverse)
250 DEALLOCATE(minus_sign_reverse)
251 CALL fvelprint(
252 . ixc, ixtg, ntgi, t_monvol(n)%ITAGEL(ntg + 1), t_monvol(n)%ELTG(ntg + 1),
253 . t_monvol(n)%IBAGHOL, ilvout , 1 )
254C------------------------------------------------------
255C Tag internal element for /TH/SURF mass flow option
256C------------------------------------------------------
257 DO i=1,nsurf
258 nseg=igrsurf(i)%NSEG
259 CALL fvthsurf(nseg, ntgi,igrsurf(i)%ELTYP, t_monvol(n)%ELTG(ntg + 1),
260 . t_monvol(n)%THSURF_TAG(i, 1:ntgi + 1), igrsurf(i)%ELEM)
261 ENDDO
262 ENDIF
263C---------------------------------
264C Additional brick nodes
265C---------------------------------
266 IF (nna > 0) THEN
267 nnfv = t_monvol(n)%IVOLU(46)
268 CALL fvnodbr(t_monvol(n)%IBUFA, nna, nnfv, ifv, nb_node)
269 ENDIF
270C
271 nnfv = t_monvol(n)%IVOLU(46)
272 ntrfv = t_monvol(n)%IVOLU(47)
273 npolh = t_monvol(n)%IVOLU(49)
274C------------------------------------
275C REFERENCE METRICS : VOLUME CHECK
276C------------------------------------
277 IF (iref /= 0) THEN
278 CALL fvvolu( ityp, nnfv, ntrfv, npolh,
279 1 t_monvol(n)%NODES, t_monvol(n)%IBUFA, t_monvol(n)%ELEMA, t_monvol(n)%TAGELA,
280 2 xyzini, t_monvol(n)%IVOLU, t_monvol(n)%RVOLU,
281 3 fvdata(ifv)%IFVNOD, fvdata(ifv)%RFVNOD, fvdata(ifv)%IFVTRI,
282 4 fvdata(ifv)%IFVPOLY,fvdata(ifv)%IFVTADR,fvdata(ifv)%IFVPOLH,
283 5 fvdata(ifv)%IFVPADR,fvdata(ifv)%MPOLH,
284 6 fvdata(ifv)%EPOLH, fvdata(ifv)%VPOLH_INI )
285C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
286
287 ENDIF
288C----------------------------------------
289C COMPUTE MINIMUM LENGTH FOR TIME STEP
290C----------------------------------------
291 CALL fvlength(nnfv, ntrfv, npolh,
292 1 t_monvol(n)%NODES, t_monvol(n)%IBUFA, t_monvol(n)%ELEMA, t_monvol(n)%TAGELA,
293 2 xyzini, t_monvol(n)%IVOLU,
294 3 fvdata(ifv)%IFVNOD, fvdata(ifv)%RFVNOD, fvdata(ifv)%IFVTRI,
295 4 fvdata(ifv)%IFVPOLY,fvdata(ifv)%IFVTADR,fvdata(ifv)%IFVPOLH,
296 5 fvdata(ifv)%IFVPADR,fvdata(ifv)%IBPOLH, fvdata(ifv)%DLH )
297C----------------------------
298C COMPUTE INITIAL SURFACE
299C----------------------------
300 CALL fvelarea(t_monvol(n)%NODES, t_monvol(n)%ELEM, xyzref, ntgt,
301 1 t_monvol(n)%ELAREA)
302C
303 ENDIF
304 ENDDO
305C
306 DEALLOCATE(x)
307 RETURN
308C
3091000 FORMAT(
310 . //,' FVMBAG: FINITE VOLUME MESH '/
311 . ' -------------------------- ')
312C
my_real
#define my_real
Definition cppsort.cpp:32
fvelarea
subroutine fvelarea(ibuf, elem, x, nel, elarea)
Definition fvelarea.F:32
fvelprint
subroutine fvelprint(ixc, ixtg, nel, itagel, eltg, ibaghol, ilvout, iflag)
Definition fvelprint.F:32
fvelsurf
subroutine fvelsurf(ibuf, elem, elem_id, ixc, ixtg, nel, eltg, mattg, nb_node, flag)
Definition fvelsurf.F:33
fvinjectint
subroutine fvinjectint(ibuf, elem, ibagjet, njet, igrsurf, itagel, nn, nel, nb_node, minus_sign_reverse)
Definition fvinjectint.F:32
fvverif
subroutine fvverif(nela, elema, x, monvid, vx3, vy3, vz3, vx1, vy1, vz1, xb0, yb0, zb0, lx, ly, lz, ibuf, ibufa, tagela, titr)
Definition fvmbag1.F:278
fvnodbr
subroutine fvnodbr(ibufa, nna, nnfv, ifv, nb_node)
Definition fvmbag1.F:530
fvthsurf
subroutine fvthsurf(nseg, ntgi, surf_eltyp, eltg, itag, surf_elem)
Definition fvthsurf.F:29
fvventholeint
subroutine fvventholeint(ibuf, elem, ibaghol, nvent, igrsurf, itagel, nn, nel, nb_node)
Definition fvventholeint.F:32
fvbag_mod::fvdata
type(fvbag_data), dimension(:), allocatable fvdata
Definition fvbag_mod.F:128
fvelinte_mod
Definition fvelinte.F:28
fvelinte_mod::fvelinte
subroutine fvelinte(ibuf, elem, ixc, ixtg, pm, ipm, ilvout, ifv, nnt, ntg, porosity, nseg, surf_eltyp, ntgi, eltg, nb_node, surf_elem)
Definition fvelinte.F:40
fvinjnormal_mod
Definition fvinjnormal.F:28
fvinjnormal_mod::fvinjnormal
subroutine fvinjnormal(ibuf, elem, ixc, ixtg, eltg, nel, ilvout, minus_sign_reverse)
Definition fvinjnormal.F:38
fvlength_mod
Definition fvlength.F:28
fvlength_mod::fvlength
subroutine fvlength(nns, nntr, npolh, ibuf, ibufa, elema, tagela, x, ivolu, ifvnod, rfvnod, ifvtri, ifvpoly, ifvtadr, ifvpolh, ifvpadr, ibpolh, dlh)
Definition fvlength.F:43
fvmbag_meshcontrol_mod::kmesh
integer, dimension(:), allocatable kmesh
Definition fvmbag_meshcontrol_mod.F:48
fvvolu_mod
Definition fvvolu.F:28
fvvolu_mod::fvvolu
subroutine fvvolu(ityp, nns, nntr, npolh, ibuf, ibufa, elema, tagela, x, ivolu, rvolu, ifvnod, rfvnod, ifvtri, ifvpoly, ifvtadr, ifvpolh, ifvpadr, mpolh, epolh, vpolh_ini)
Definition fvvolu.F:44
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
fvbric
subroutine fvbric(ivolu, rvolu, ibuf, x, nn)
Definition fvbric.F:31
fvmesh1
subroutine fvmesh1(ibuf, elem, x, ivolu, bric, xb, rvolu, nel, neli, nbric, tbric, sfac, dxm, nba, nela, tba, tfaca, tagels, ibufa, elema, tagela, ixs, id, titr, nb_node, ityp)
Definition fvmesh.F:53