OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
gendynain.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!|| gendynain ../engine/source/output/dynain/gendynain.F
25!||--- called by ------------------------------------------------------
26!|| sortie_main ../engine/source/output/sortie_main.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| arret ../engine/source/system/arret.F
30!|| close_c ../common_source/tools/input_output/write_routines.c
31!|| cur_fil_c ../common_source/tools/input_output/write_routines.c
32!|| dynain_c_strag ../engine/source/output/dynain/dynain_c_strag.F
33!|| dynain_c_strsg ../engine/source/output/dynain/dynain_c_strsg.F
34!|| dynain_node ../engine/source/output/dynain/dynain_node.F
35!|| dynain_shel_mp ../engine/source/output/dynain/dynain_shel_mp.F
36!|| dynain_shel_spmd ../engine/source/output/dynain/dynain_shel_spmd.F
37!|| dynain_size_c ../engine/source/output/dynain/dynain_size.f
38!|| open_c ../common_source/tools/input_output/write_routines.c
39!|| spmd_outpitab ../engine/source/mpi/interfaces/spmd_outp.F
40!|| strs_txt50 ../engine/source/output/sta/sta_txt.F
41!||--- uses -----------------------------------------------------
42!|| drape_mod ../engine/share/modules/drape_mod.F
43!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
44!|| element_mod ../common_source/modules/elements/element_mod.F90
45!|| inoutfile_mod ../common_source/modules/inoutfile_mod.F
46!|| matparam_def_mod ../common_source/modules/mat_elem/matparam_def_mod.F90
47!|| message_mod ../engine/share/message_module/message_mod.F
48!|| stack_mod ../engine/share/modules/stack_mod.F
49!|| state_mod ../common_source/modules/state_mod.F
50!||====================================================================
51 SUBROUTINE gendynain(X ,ELBUF_TAB, BUFEL ,IXC ,IXTG ,
52 2 IPARG ,IPM , IGEO ,ITAB ,IPART ,
53 3 PM ,GEO , IPARTC ,IPARTTG ,LENG ,
54 4 LENGC ,LENGTG , WEIGHT ,NODGLOB ,THKE ,
55 5 NPBY ,LPBY , STACK ,DRAPE_SH4N ,DRAPE_SH3N ,
56 6 DYNAIN_DATA,DRAPEG ,MAT_PARAM)
57C-----------------------------------------------
58C M o d u l e s
59C-----------------------------------------------
60 USE elbufdef_mod
61 USE matparam_def_mod
62 USE message_mod
63 USE inoutfile_mod
64 USE stack_mod
65 USE drape_mod
66 USE state_mod
67 use element_mod , only : nixc,nixtg
68C-----------------------------------------------
69C I m p l i c i t T y p e s
70C-----------------------------------------------
71#include "implicit_f.inc"
72C-----------------------------------------------
73C C o m m o n B l o c k s
74C-----------------------------------------------
75#include "com01_c.inc"
76#include "com04_c.inc"
77#include "param_c.inc"
78#include "units_c.inc"
79#include "scr17_c.inc"
80#include "chara_c.inc"
81#include "task_c.inc"
82C-----------------------------------------------
83C D u m m y A r g u m e n t s
84C-----------------------------------------------
85 INTEGER IPARG(*),
86 . IXC(NIXC,*),IXTG(NIXTG,*),IPM(*),IGEO(*),
87 . ITAB(*) ,IPART(LIPART1,*) ,IPARTC(*) ,IPARTTG(*),
88 . WEIGHT(*), NODGLOB(*), NPBY(NNPBY,*), LPBY(*)
89 INTEGER LENG,LENGC,LENGTG
90 my_real
91 . x(*), bufel(*),
92 . pm(npropm,*), geo(npropg,*) ,thke(*)
93 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
94 TYPE (STACK_PLY) :: STACK
95 TYPE (DRAPE_) :: DRAPE_SH4N(NUMELC_DRAPE),DRAPE_SH3N(NUMELTG_DRAPE)
96 TYPE (DRAPEG_) :: DRAPEG
97 TYPE (DYNAIN_DATABASE), INTENT(INOUT) :: DYNAIN_DATA
98 TYPE (MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(IN) :: MAT_PARAM
99C-----------------------------------------------
100C L o c a l V a r i a b l e s
101C-----------------------------------------------
102 CHARACTER CHSTAT*4, FILNAM*100, T10*10, MES*40
103 INTEGER FILEN, I, IERR, J, N
104 INTEGER SIZLOC, SIZP0
105 INTEGER , DIMENSION(:),ALLOCATABLE :: ITABG, NODTAG ,DYNAIN_INDXC ,
106 . DYNAIN_INDXTG
107 INTEGER CTEXT(2149)
108 double precision
109 . , DIMENSION(:),ALLOCATABLE :: wa,wap0
110
111 INTEGER :: LEN_TMP_NAME
112 CHARACTER(len=2148) :: TMP_NAME
113
114 CHARACTER*100 LINE
115C-----------------------------------------------
116C S o u r c e L i n e s
117C-----------------------------------------------
118C===============================================|
119C OPEN FILE
120C-----------------------------------------------
121 IF(dynain_data%IDYNAINF>=10000)dynain_data%IDYNAINF=1
122 WRITE(chstat,'(I4.4)')dynain_data%IDYNAINF
123 IF(dynain_data%ZIPDYNAIN==0) THEN ! Not zipped file
124 filnam=rootnam(1:rootlen)//'_'//chstat//'.dynain'
125 filen = rootlen + 12
126 len_tmp_name = outfile_name_len + filen
127 tmp_name=outfile_name(1:outfile_name_len)//filnam(1:filen)
128 IF(ispmd == 0) THEN
129 OPEN(unit=iudynain,file=tmp_name(1:len_tmp_name),access='SEQUENTIAL',form='FORMATTED',status='UNKNOWN')
130 WRITE(iudynain,'(2A)')'$RADIOSS DYNAIN FILE ',filnam(1:filen)
131 END IF
132 ELSE ! zipped file
133 filnam=rootnam(1:rootlen)//'_'//chstat//'.dynain'
134 filen = rootlen + 12
135 len_tmp_name = outfile_name_len + filen
136 tmp_name=outfile_name(1:outfile_name_len)//filnam(1:filen)
137 DO i=1,len_tmp_name
138 ctext(i)=ichar(tmp_name(i:i))
139 CALL cur_fil_c(0)
140 ENDDO
141 IF(ispmd == 0) THEN
142 CALL open_c(ctext,len_tmp_name,6)
143 WRITE(line,'(2A)') '$RADIOSS DYNAIN FILE ',filnam(1:filen)
144 CALL strs_txt50(line,100)
145 ENDIF
146 ENDIF
147
148c
149
150C-----------------------
151C Allocation Tabs
152C-----------------------
153 ALLOCATE(nodtag(numnod),stat=ierr)
154 ALLOCATE(itabg(leng),stat=ierr)
155 ALLOCATE(dynain_indxc(2*lengc),stat=ierr)
156 ALLOCATE(dynain_indxtg(2*lengtg),stat=ierr)
157C-----------------------------------------------
158C CONNECTIVITIES + NODAL COORDINATES
159C-----------------------------------------------
160c
161 IF (nspmd > 1)CALL spmd_outpitab(itab,weight,nodglob,itabg)
162
163 nodtag=0
164
165 dynain_data%DYNAIN_NUMELC =0
166 dynain_data%DYNAIN_NUMELTG =0
167
168 IF(nspmd == 1)THEN
169 ! - shells -
170 CALL dynain_shel_mp(itab ,itabg ,leng ,igeo ,ixc ,
171 . ixtg ,ipartc ,iparttg ,dynain_data ,
172 . nodtag ,dynain_indxc,dynain_indxtg,iparg ,
173 . elbuf_tab,thke ,ipart )
174 dynain_data%DYNAIN_NUMELC_G =dynain_data%DYNAIN_NUMELC
175 dynain_data%DYNAIN_NUMELTG_G =dynain_data%DYNAIN_NUMELTG
176 ELSE
177 ! - shells -
178 CALL dynain_shel_spmd(itab ,itabg ,leng ,igeo ,ixc ,
179 . ixtg ,ipartc ,iparttg ,dynain_data ,
180 . nodtag ,dynain_indxc,dynain_indxtg,iparg ,
181 . elbuf_tab,thke ,lengc ,lengtg ,ipart )
182 END IF
183
184
185C-----------------------------------------------
186C RIGID BODY'S PRIMARY NODE IS OUTPUTTED IF ONE OF ITS SECONDARY NODES ARE
187C-----------------------------------------------
188 DO i=1,nrbody
189 DO j=1,npby(2,i)
190 n=lpby(npby(11,i)+j)
191 IF (nodtag(n)/=0) THEN
192 nodtag(npby(1,i)) = 1
193 EXIT
194 END IF
195 ENDDO
196 ENDDO
197C-----------------------------------------------
198 CALL dynain_node(x,numnod,itab,itabg,leng,nodglob,weight,nodtag,dynain_data)
199
200C-----------------------------------------------
201 CALL dynain_size_c(iparg ,elbuf_tab, sizp0 ,sizloc ,dynain_data )
202
203C-----------------------------------------------
204C ALLOCATION OF TABLES
205C-----------------------------------------------
206 ierr = 0
207 IF(sizloc >= 1) THEN
208 ALLOCATE(wa(sizloc),stat=ierr)
209 ELSE
210 ALLOCATE(wa(1))
211 ENDIF
212 IF(ierr/=0)THEN
213 CALL ancmsg(msgid=252,anmode=aninfo,
214 . i1=ierr)
215 CALL arret(2)
216 END IF
217
218 ierr = 0
219 sizp0 = max(1,sizp0)
220 ALLOCATE(wap0(sizp0),stat=ierr)
221 IF(ierr/=0)THEN
222 CALL ancmsg(msgid=252,anmode=aninfo,
223 . i1=ierr)
224 CALL arret(2)
225 END IF
226C-----------------------------------------------
227C SHELL SCALAR
228C-----------------------------------------------
229
230
231 IF(dynain_data%DYNAIN_C(4)==1) THEN
232 CALL dynain_c_strsg(
233 1 elbuf_tab ,iparg ,igeo ,ixc ,
234 2 ixtg ,wa ,wap0 ,ipartc,iparttg,
235 3 dynain_data,dynain_indxc,dynain_indxtg,sizp0 ,
236 4 geo ,stack ,drape_sh4n ,drape_sh3n,x ,
237 5 thke , drapeg ,nummat ,mat_param )
238 ENDIF
239C------------------------------------------
240
241 IF(dynain_data%DYNAIN_C(5)==1) THEN
242 CALL dynain_c_strag(
243 1 elbuf_tab ,iparg ,ipm ,igeo ,ixc ,
244 2 ixtg ,wa ,wap0 ,ipartc,iparttg,
245 3 dynain_data,dynain_indxc,dynain_indxtg,sizp0 ,
246 4 geo ,stack ,drape_sh4n ,drape_sh3n,x ,
247 5 thke ,drapeg)
248 ENDIF
249
250C-----------------------------------------------
251
252 IF(sizloc >= 1) DEALLOCATE(wa)
253 IF(sizp0 >= 1) DEALLOCATE(wap0)
254C-----------------------
255C DEAllocation Tabs
256C-----------------------
257 DEALLOCATE(nodtag,itabg,dynain_indxc,dynain_indxtg)
258C-----------------------------------------------
259C END
260C-----------------------------------------------
261 IF(ispmd==0) THEN
262 IF(dynain_data%ZIPDYNAIN==0) THEN
263 WRITE(iudynain,'(A)')'*END '
264 CLOSE(unit=iudynain)
265 ELSE
266 CALL strs_txt50('*END ',7)
267 CALL close_c()
268 ENDIF
269
270 WRITE (iout,1000) filnam(1:filen)
271 WRITE (istdo,1000) filnam(1:filen)
272 ENDIF
273
274
275 1000 FORMAT (4x,' DYNAIN FILE:',1x,a,' WRITTEN')
276C
277 RETURN
278 END
my_real
#define my_real
Definition cppsort.cpp:32
dynain_c_strag
subroutine dynain_c_strag(elbuf_tab, iparg, ipm, igeo, ixc, ixtg, wa, wap0, ipartc, iparttg, dynain_data, dynain_indxc, dynain_indxtg, sizp0, geo, stack, drape_sh4n, drape_sh3n, x, thke, drapeg)
Definition dynain_c_strag.F:48
dynain_c_strsg
subroutine dynain_c_strsg(elbuf_tab, iparg, igeo, ixc, ixtg, wa, wap0, ipartc, iparttg, dynain_data, dynain_indxc, dynain_indxtg, sizp0, geo, stack, drape_sh4n, drape_sh3n, x, thke, drapeg, nummat, mat_param)
Definition dynain_c_strsg.F:49
dynain_node
subroutine dynain_node(x, numnod, itab, itabg, leng, nodglob, weight, nodtag, dynain_data)
Definition dynain_node.F:35
dynain_shel_mp
subroutine dynain_shel_mp(itab, itabg, leng, igeo, ixc, ixtg, ipartc, iparttg, dynain_data, nodtag, dynain_indxc, dynain_indxtg, iparg, elbuf_tab, thke, ipart)
Definition dynain_shel_mp.F:39
dynain_shel_spmd
subroutine dynain_shel_spmd(itab, itabg, leng, igeo, ixc, ixtg, ipartc, iparttg, dynain_data, nodtag, dynain_indxc, dynain_indxtg, iparg, elbuf_tab, thke, lengc, lengtg, ipart)
Definition dynain_shel_spmd.F:41
dynain_size_c
subroutine dynain_size_c(iparg, elbuf_tab, p0ars, wasz, dynain_data)
Definition dynain_size.F:34
gendynain
subroutine gendynain(x, elbuf_tab, bufel, ixc, ixtg, iparg, ipm, igeo, itab, ipart, pm, geo, ipartc, iparttg, leng, lengc, lengtg, weight, nodglob, thke, npby, lpby, stack, drape_sh4n, drape_sh3n, dynain_data, drapeg, mat_param)
Definition gendynain.F:57
max
#define max(a, b)
Definition macros.h:21
drape_mod
Definition drape_mod.F:54
inoutfile_mod
Definition inoutfile_mod.F:96
inoutfile_mod::outfile_name
character(len=outfile_char_len) outfile_name
Definition inoutfile_mod.F:111
inoutfile_mod::outfile_name_len
integer outfile_name_len
Definition inoutfile_mod.F:109
message_mod
Definition message_mod.F:1257
stack_mod
Definition stack_mod.F:75
state_mod
Definition state_mod.F:54
spmd_outpitab
subroutine spmd_outpitab(v, weight, nodglob, vglob)
Definition spmd_outp.F:1081
strs_txt50
subroutine strs_txt50(text, length)
Definition sta_txt.F:87
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895
arret
subroutine arret(nn)
Definition arret.F:86
close_c
void close_c()
Definition write_routines.c:298
cur_fil_c
void cur_fil_c(int *nf)
Definition write_routines.c:477
open_c
void open_c(int *ifil, int *len, int *mod)
Definition write_routines.c:180