OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
anioffc_ply.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!|| anioffc_ply ../engine/source/output/anim/generate/anioffc_ply.F
25!||--- called by ------------------------------------------------------
26!|| genani ../engine/source/output/anim/generate/genani.F
27!||--- calls -----------------------------------------------------
28!|| spmd_iget_partn ../engine/source/mpi/anim/spmd_iget_partn.F
29!|| write_c_c ../common_source/tools/input_output/write_routtines.c
30!||--- uses -----------------------------------------------------
31!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
32!|| my_alloc_mod ../common_source/tools/memory/my_alloc.F90
33!|| plyxfem_mod ../engine/share/modules/plyxfem_mod.F
34!|| stack_mod ../engine/share/modules/stack_mod.F
35!||====================================================================
36 SUBROUTINE anioffc_ply( IPLY, NEL_PLY, ELBUF_TAB, IPARG,
37 . IOFF, EL2FA, NBF, IADD,
38 . NBF_L, NBPART, IADG, NODGLOB,
39 . IPART, IPARTC, IPARTTG, NBF_PXFEMG,
40 . IPM ,IGEO ,IXC, STACK )
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE elbufdef_mod
45 USE plyxfem_mod
46 USE stack_mod
47 USE my_alloc_mod
48C-----------------------------------------------
49C I m p l i c i t T y p e s
50C-----------------------------------------------
51#include "implicit_f.inc"
52C-----------------------------------------------
53C C o m m o n B l o c k s
54C-----------------------------------------------
55#include "mvsiz_p.inc"
56#include "com01_c.inc"
57#include "com04_c.inc"
58#include "param_c.inc"
59#include "remesh_c.inc"
60#include "scr17_c.inc"
61#include "task_c.inc"
62C-----------------------------------------------
63C D u m m y A r g u m e n t s
64C-----------------------------------------------
65 INTEGER IPARG(NPARG,*),EL2FA(*),NBF,IOFF(*),
66 . IADD(*),NBF_L,NBPART, IADG(NSPMD,*),NODGLOB(*),
67 . IPART(LIPART1,*), IPARTC(*), IPARTTG(*)
68 INTEGER IPLY,NEL_PLY,PLYS,NBF_PXFEMG,IXC(NIXC,*),IPM(NPROPMI,*),
69 . igeo(npropgi,*)
70 TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
71 TYPE (STACK_PLY) :: STACK
72C-----------------------------------------------
73C L o c a l V a r i a b l e s
74C-----------------------------------------------
75 INTEGER I, NG, NEL, NFT, IAD, ITY, LFT, NPT,
76 . N, J, LLT, MLW, IP, NN, K1, K2,MT,JALE, IMID,
77 . N1,N2,N3,N4, NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9,NN10,
78 . N_FA, IHBE, SH_IH, ISTRAIN, IEXPAN,ISHPLYXFEM,NELXFE,IE,
79 . plyelems(numelc),ipt,elc,mpt,ngl(mvsiz),isubstack,ipmat_iply,
80 . ipmat, ippid
81
82 INTEGER RBUF,NEL3,NEL5,NEL8,NPTM,NN12,NN13,NN14,
83 . NN15,NN16,NN16A,NPG,PID(MVSIZ),MATLY,NUVARR,NUVARD,NBM_S,
84 . IFAILURE,NPTR,NPTS,IS,IR,OFFIPLY
85 INTEGER,DIMENSION(:),ALLOCATABLE::IOFFBUF
86 REAL R4
87 TYPE(G_BUFEL_) ,POINTER :: GBUF
88 TYPE(BUF_LAY_) ,POINTER :: BUFLY
89 TYPE(l_bufel_) ,POINTER :: LBUF
90C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
91 CALL my_alloc(ioffbuf,nbf_l)
92 npg = 4
93 ioffbuf = -1
94 DO plys = 1,nplypart
95 iply = indx_ply(plys)
96
97C-----------------------------------------------
98C PLY ELEMS INIT
99C-----------------------------------------------
100 plyelems=0
101 DO i=1,plyshell(iply)%PLYNUMSHELL
102 ipt = plyshell(iply)%SHELLIPT(i)
103 elc = plyshell(iply)%SHID(i)
104 plyelems(elc)=ipt
105 ENDDO
106
107C-----------------------------------------------
108C
109cc NN1 = EPLYXFE*(IPLY - 1)
110
111C-----------------------------------------------
112C
113 nelxfe = 0
114 ie = 0
115 DO ng=1,ngroup
116 mlw =iparg(1,ng)
117 nel =iparg(2,ng)
118 ity =iparg(5,ng)
119 nft =iparg(3,ng)
120 iad =iparg(4,ng)
121 mpt = iparg(6,ng)
122 ishplyxfem = iparg(50,ng)
123 isubstack = iparg(71,ng)
124 lft=1
125 llt=nel
126C-----------------------------------------------
127C COQUES 4 N
128C-----------------------------------------------
129 IF(ishplyxfem > 0 ) THEN
130C for batoz shell
131 nptr = elbuf_tab(ng)%NPTR
132 npts = elbuf_tab(ng)%NPTS
133C
134 gbuf => elbuf_tab(ng)%GBUF
135 bufly => elbuf_tab(ng)%BUFLY(ipt) ! Ipt is no right
136 IF (ity == 3) THEN
137 ihbe = iparg(23,ng)
138 npt =iparg(6,ng)
139 istrain=iparg(44,ng)
140 iexpan=iparg(49,ng)
141 ifailure = iparg(43,ng)
142 n_fa =nel_ply
143 nelxfe = nelxfe + nel
144
145 DO i=1,nel
146 pid(i) = ixc(6,nft + i)
147 ngl(i) = ixc(7,nft + i)
148 ENDDO
149C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
150 IF (ihbe >= 11.AND.ihbe <= 19)THEN
151 nel3 = nel*3
152 nel5 = nel*5
153 nel8 = nel*8
154 nptm = max(1,mpt)
155 IF (mpt >= one) THEN
156 ENDIF
157 nbm_s = 6*nel*mpt+nel
158 nuvarr = zero
159 nuvard = zero
160C
161 ippid = 2
162 ipmat = ippid + npt ! layer material address ( NLAY = NPT )
163 ipmat_iply = ipmat + npt
164 DO j= 1,mpt-1
165 DO i=1,nel
166 matly = stack%IGEO(ipmat + j,isubstack)
167 nuvarr = max(nuvarr,ipm(221,matly))
168 matly = stack%IGEO(ipmat_iply + j,isubstack)
169 nuvard = max(nuvard, ipm(221,matly))
170 ENDDO
171 ENDDO
172 DO i=1,nel
173 matly = stack%IGEO(ipmat + j,isubstack)
174 nuvarr = max(nuvarr,ipm(221,matly))
175 ENDDO
176 ENDIF
177 nbm_s = nbm_s + nel*mpt*nuvarr
178C-----------------------------------------------
179C COQUES 3 N
180C-----------------------------------------------
181 ELSEIF(ity==7)THEN
182c it's not available for T3
183 ENDIF
184C-----------------------------------------------
185 IF(ity == 3)THEN
186 IF(mlw == 0 .OR. mlw == 13)THEN
187C-----------------------------------------------
188C DUMMY ELEMENTS
189C-----------------------------------------------
190 DO i=lft,llt
191 n = i + nft
192 IF(plyelems(n) > 0) THEN
193 ie = ie + 1
194 ioff(el2fa(n_fa+ie)) = 1
195 ENDIF
196 ENDDO
197 ELSE
198C-----------------------------------------------
199C OFF Is Only for batoz for now. Because Plyxfem is available only with Bato z shell
200C-----------------------------------------------
201 IF(nadmesh==0.OR.(ity/=3.AND.ity/=7))THEN
202 DO i=lft,llt
203 n = i + nft
204 ipt = plyelems(n)
205 IF(plyelems(n) > 0) THEN
206 ie = ie +1
207 ioff(el2fa(n_fa+ie)) = 1
208 DO is = 1,npts
209 DO ir = 1,nptr
210 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(ir,is,1)
211 offiply = nint(min(gbuf%OFF(i),lbuf%OFF(i),one))
212 ioff(el2fa(n_fa+ie)) = min(ioff(el2fa(n_fa+ie)),offiply)
213 ENDDO
214 ENDDO
215 ENDIF
216 ENDDO
217 ELSE
218 DO i=lft,llt
219 n = i + nft
220 ipt = plyelems(n)
221 IF(plyelems(n) > 0) THEN
222 ip=ipartc(nft+i)
223 ie = ie + 1
224 ioff(el2fa(n_fa+ie)) = 1
225 IF(ipart(10,ip)>0)THEN
226 DO is = 1,npts
227 DO ir = 1,nptr
228 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(ir,is,1)
229 offiply = nint(min(gbuf%OFF(i),lbuf%OFF(i),one))
230 ioff(el2fa(n_fa+ie))=min(ioff(el2fa(n_fa+ie)),offiply)
231 ENDDO
232 ENDDO
233 ELSE
234 DO is = 1,npts
235 DO ir = 1,nptr
236 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(ir,is,1)
237 offiply = nint(min(gbuf%OFF(i),lbuf%OFF(i),one))
238 ioff(el2fa(n_fa+ie)) =min(ioff(el2fa(n_fa+ie)),offiply)
239 ENDDO
240 ENDDO
241 END IF
242 ENDIF
243 ENDDO
244 ENDIF
245 ENDIF
246 ENDIF
247C
248 ENDIF
249C-----------------------------------------------
250 ENDDO
251C-----------------------------------------------
252 iadd(iply) = nel_ply + ie
253 nel_ply = nel_ply + plyshell(iply)%PLYNUMSHELL
254 ENDDO
255C
256 IF (nspmd==1)THEN
257 CALL write_c_c(ioff,nbf)
258 ELSE
259
260 DO i = 1, nbf_l
261 ioffbuf(i) = ioff(i)
262 ENDDO
263
264 IF (ispmd==0) THEN
265 rbuf = nbf_pxfemg
266 CALL spmd_iget_partn(1,nbf_l,ioffbuf,nplypart,iadg,rbuf,2)
267 ELSE
268 rbuf = 1
269 CALL spmd_iget_partn(1,nbf_l,ioffbuf,nplypart,iadg,rbuf,2)
270 END IF
271 ENDIF
272 DEALLOCATE(ioffbuf)
273 RETURN
274 END
anioffc_ply
subroutine anioffc_ply(iply, nel_ply, elbuf_tab, iparg, ioff, el2fa, nbf, iadd, nbf_l, nbpart, iadg, nodglob, ipart, ipartc, iparttg, nbf_pxfemg, ipm, igeo, ixc, stack)
Definition anioffc_ply.F:41
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
plyxfem_mod
Definition plyxfem_mod.F:31
plyxfem_mod::indx_ply
integer, dimension(:), allocatable indx_ply
Definition plyxfem_mod.F:60
plyxfem_mod::nplypart
integer nplypart
Definition plyxfem_mod.F:59
plyxfem_mod::plyshell
type(plyshells), dimension(:), allocatable plyshell
Definition plyxfem_mod.F:56
stack_mod
Definition stack_mod.F:72
spmd_iget_partn
subroutine spmd_iget_partn(size, nbf_l, np, nbpart, iadg, srbuf, iflag)
Definition spmd_iget_partn.F:50
write_c_c
void write_c_c(int *w, int *len)
Definition write_routtines.c:591