OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
dfuncs.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!|| dfuncs ../starter/source/output/anim/dfuncs.F
25!||--- called by ------------------------------------------------------
26!|| genani1 ../starter/source/output/anim/genani1.F
27!||--- calls -----------------------------------------------------
28!|| initbuf ../starter/source/output/anim/initbuf.F
29!||--- uses -----------------------------------------------------
30!||====================================================================
31 SUBROUTINE dfuncs(ELBUF_TAB,FUNC ,IFUNC ,IPARG ,
32 2 IXS ,PM ,EL2FA ,NBF ,ISPH3D )
33C-----------------------------------------------
34C M o d u l e s
35C-----------------------------------------------
36 USE elbufdef_mod
37C-----------------------------------------------
38C I m p l i c i t T y p e s
39C-----------------------------------------------
40#include "implicit_f.inc"
41C-----------------------------------------------
42C C o m m o n B l o c k s
43C-----------------------------------------------
44#include "vect01_c.inc"
45#include "mvsiz_p.inc"
46#include "com01_c.inc"
47#include "com04_c.inc"
48#include "param_c.inc"
49C-----------------------------------------------
50C D u m m y A r g u m e n t s
51C-----------------------------------------------
52C REAL
53 my_real
54 . func(*), pm(npropm,*)
55 INTEGER IPARG(NPARG,*),EL2FA(*),
56 . ixs(nixs,*),ifunc,nbf,isph3d
57 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
58 REAL WAL(NBF)
59C-----------------------------------------------
60C L o c a l V a r i a b l e s
61C-----------------------------------------------
62C REAL
63 my_real
64 . evar(mvsiz),
65 . off, p, vonm2, vonm, s1, s2, s12, s3, VALUE
66 INTEGER I, NG, NEL,N, J, MLW,NN, JTURB,MT, IALEL,
67 . nn1,nn2,nn3,offset,ii(6),inod, isolnod,
68 . jhbe, jivf, jclose, jplasol, irep, igtyp,
69 . icsen, isorthg, ifailure, iint
70 TYPE(g_bufel_) ,POINTER :: GBUF
71 REAL R4
72C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
73C La routine ne fonctionne que pour les IFUNC 3,6,7,14-19 (stress)
74C-----------------------------------------------
75 nn1 = 1
76 nn2 = 1
77 nn3 = nn2 + numels
78C
79 DO 900 ng=1,ngroup
80 CALL initbuf (iparg ,ng ,
81 2 mlw ,nel ,nft ,iad ,ity ,
82 3 npt ,jale ,ismstr ,jeul ,jtur ,
83 4 jthe ,jlag ,jmult ,jhbe ,jivf ,
84 5 jpor ,jcvt ,jclose ,jplasol ,
85 6 irep ,iint ,igtyp ,israt ,isrot ,
86 7 icsen ,isorth ,isorthg ,ifailure)
87 DO offset = 0,nel-1,nvsiz
88 nft =iparg(3,ng) + offset
89 isolnod = iparg(28,ng)
90 lft=1
91 llt=min(nvsiz,nel-offset)
92!
93 DO i=1,6
94 ii(i) = (i-1)*llt
95 ENDDO
96!
97C-----------------------------------------------
98C SOLID 8N
99C-----------------------------------------------
100 IF (ity == 1) THEN
101C-----------
102 gbuf => elbuf_tab(ng)%GBUF
103c
104 IF (mlw == 0 .OR. mlw == 13 . or. igtyp == 0) THEN
105 DO i=lft,llt
106 evar(i) = zero
107 ENDDO
108 ELSE
109 jturb=iparg(12,ng)*(iparg(7,ng)+iparg(11,ng))
110C
111 IF (ifunc == 1) THEN
112 DO i=lft,llt
113 IF (gbuf%G_PLA > 0) THEN
114 evar(i) = gbuf%PLA(i)
115 ENDIF
116 ENDDO
117 ELSEIF(ifunc == 2)THEN
118 DO i=lft,llt
119 evar(i) = gbuf%RHO(i)
120 ENDDO
121 ELSEIF(ifunc == 3)THEN
122 DO i=lft,llt
123 n = i + nft
124 ialel=iparg(7,ng)+iparg(11,ng)
125 IF (ialel == 0) THEN
126 mt=ixs(1,n)
127 evar(i) = gbuf%EINT(i)/max(em30,pm(1,mt))
128 ELSE
129 evar(i) = gbuf%EINT(i)/max(em30,gbuf%RHO(i))
130 ENDIF
131 ENDDO
132 ELSEIF (ifunc == 4) THEN
133 DO i=lft,llt
134 IF (gbuf%G_TEMP > 0) THEN
135 evar(i) = gbuf%TEMP(i)
136 ENDIF
137 ENDDO
138 ELSEIF(ifunc == 6 .OR. ifunc == 7)THEN
139 DO i=lft,llt
140 n = i + nft
141 p = - (gbuf%SIG(ii(1) + i)
142 . + gbuf%SIG(ii(2) + i)
143 . + gbuf%SIG(ii(3) + i)) * third
144 VALUE = p
145 IF (ifunc == 7) THEN
146 s1=gbuf%SIG(ii(1) + i)+p
147 s2=gbuf%SIG(ii(2) + i)+p
148 s3=gbuf%SIG(ii(3) + i)+p
149 vonm2= three*(gbuf%SIG(ii(4) + i)**2 +
150 . gbuf%SIG(ii(5) + i)**2 +
151 . gbuf%SIG(ii(6) + i)**2 +
152 . half*(s1*s1+s2*s2+s3*s3) )
153 vonm= sqrt(vonm2)
154 VALUE = vonm
155 ENDIF
156 evar(i) = VALUE
157 ENDDO
158c-----------
159 ELSEIF(ifunc >= 14 .AND. ifunc <= 19)THEN
160 DO i=lft,llt
161 evar(i) = gbuf%SIG(ii(ifunc-13) + i)
162 ENDDO
163 ENDIF
164c-----------
165 IF (isolnod == 16) THEN
166 DO i=lft,llt
167 n = nn2 + i + nft
168 IF(el2fa(n)/=0)THEN
169 func(el2fa(n)) = evar(i)
170 func(el2fa(n)+1) = evar(i)
171 func(el2fa(n)+2) = evar(i)
172 func(el2fa(n)+3) = evar(i)
173 ENDIF
174 ENDDO
175 ELSE
176 DO i=lft,llt
177 n = nn2 + i + nft
178 IF(el2fa(n)/=0)THEN
179 func(el2fa(n)) = evar(i)
180 ENDIF
181 ENDDO
182 ENDIF
183 ENDIF
184C
185C-----------------------------------------------
186 ELSEIF (isph3d == 1.AND.ity == 51) THEN
187C TETRAS SPH.
188C-----------------------------------------------
189 gbuf => elbuf_tab(ng)%GBUF
190 IF (ifunc >= 14 .AND. ifunc <= 19) THEN
191 DO i=lft,llt
192 n = i + nft
193 IF (el2fa(nn3+n)/=0) THEN
194 func(el2fa(nn3+n)) = gbuf%SIG(ii(ifunc-13) + i)
195 ENDIF
196 ENDDO
197 ELSE
198 DO i=lft,llt
199 n = i + nft
200 IF(el2fa(nn3+n)/=0)THEN
201 func(el2fa(nn3+n)) = zero
202 ENDIF
203 ENDDO
204 ENDIF
205C
206 ENDIF
207C-----------------------------------------------
208C FIN DE BOUCLE SUR LES OFFSET
209C-----------------------------------------------
210 ENDDO
211 900 CONTINUE
212C-----------------------------------------------
213 DO n=1,nbf
214 r4 = func(n)
215 CALL write_r_c(r4,1)
216 ENDDO
217C-----------
218 RETURN
219 END
my_real
#define my_real
Definition cppsort.cpp:32
dfuncs
subroutine dfuncs(elbuf_tab, func, ifunc, iparg, ixs, pm, el2fa, nbf, isph3d)
Definition dfuncs.F:33
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
initbuf
subroutine initbuf(iparg, ng, mtn, llt, nft, iad, ity, npt, jale, ismstr, jeul, jtur, jthe, jlag, jmult, jhbe, jivf, jpor, jcvt, jclose, jpla, irep, iint, igtyp, israt, isrot, icsen, isorth, isorthg, ifailure)
Definition initbuf.F:38
write_r_c
void write_r_c(float *w, int *len)
Definition write_routtines.c:519