OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
thickvar.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!|| thickvar ../engine/source/elements/shell/coque/thickvar.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- calls -----------------------------------------------------
28!|| startimeg ../engine/source/system/timer.F
29!|| stoptimeg ../engine/source/system/timer.F
30!||--- uses -----------------------------------------------------
31!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
32!|| element_mod ../common_source/modules/elements/element_mod.F90
33!||====================================================================
34 SUBROUTINE thickvar(IPARG,ELBUF_TAB,IXC,IXTG,THKSH4_VAR,
35 . THKSH3_VAR,THKNOD,THKE,SH4TREE,SH3TREE)
36C-----------------------------------------------
37C M o d u l e s
38C-----------------------------------------------
39 USE elbufdef_mod
40 use element_mod , only : nixc,nixtg
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45#include "comlock.inc"
46C-----------------------------------------------
47C C o m m o n B l o c k s
48C-----------------------------------------------
49#include "com01_c.inc"
50#include "com04_c.inc"
51#include "param_c.inc"
52#include "remesh_c.inc"
53#include "task_c.inc"
54#include "vect01_c.inc"
55C-----------------------------------------------
56C D u m m y A r g u m e n t s
57C-----------------------------------------------
58 INTEGER IPARG(NPARG,*), IXC(NIXC,*), IXTG(NIXTG,*),
59 . sh4tree(ksh4tree,*), sh3tree(ksh3tree,*)
60C REAL
61 my_real
62 . thksh4_var(*), thksh3_var(*), thknod(*),
63 . thke(*)
64 TYPE(elbuf_struct_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
65C-----------------------------------------------
66C L o c a l V a r i a b l e s
67C-----------------------------------------------
68 INTEGER NG, I, J, N, NEL, ISH3N, IHBE, ITHK,
69 . le,lelt,nelt(2*(4**levelmax)),lev,ne,son,lelt1,lelt2
70 my_real
71 . thk
72 TYPE(g_bufel_) ,POINTER :: GBUF
73C-----------------------------------------------
74 IF(nadmesh==0)THEN
75 250 CONTINUE
76#include "lockon.inc"
77 IF(nsgdone>ngroup) THEN
78#include "lockoff.inc"
79 GOTO 252
80 ENDIF
81 ng=nsgdone
82 nsgdone = ng + 1
83#include "lockoff.inc"
84 IF (iparg(8,ng) == 1)GOTO 250
85 ity =iparg(5,ng)
86 IF(ity/=3.AND.ity/=7)GOTO 250
87
88 IF (iddw>0) CALL startimeg(ng)
89 mtn=iparg(1,ng)
90 IF(mtn == 0) GOTO 250
91 nel=iparg(2,ng)
92 nft=iparg(3,ng)
93 iad=iparg(4,ng)
94 ity=iparg(5,ng)
95 ihbe =iparg(23,ng)
96 ish3n =iparg(23,ng)
97 ithk = iparg(28,ng)
98 gbuf => elbuf_tab(ng)%GBUF
99c
100 IF (ity == 3) THEN
101 IF(ithk==0)THEN
102 DO i=1,nel
103 n=nft+i
104 thksh4_var(n)=thke(n)
105 DO j=2,5
106 thknod(ixc(j,n))=max(thknod(ixc(j,n)),thksh4_var(n))
107 END DO
108 END DO
109 ELSE
110 DO i=1,nel
111 n=nft+i
112 thksh4_var(n)=gbuf%THK(i)
113 DO j=2,5
114 thknod(ixc(j,n))=max(thknod(ixc(j,n)),thksh4_var(n))
115 END DO
116 END DO
117 END IF
118 ELSEIF(ity == 7)THEN
119 IF(ithk==0)THEN
120 DO i=1,nel
121 n=nft+i
122 thksh3_var(n)=thke(numelc+n)
123 DO j=2,4
124 thknod(ixtg(j,n))=max(thknod(ixtg(j,n)),thksh3_var(n))
125 END DO
126 END DO
127 ELSE
128 DO i=1,nel
129 n=nft+i
130 thksh3_var(n)=gbuf%THK(i)
131 DO j=2,4
132 thknod(ixtg(j,n))=max(thknod(ixtg(j,n)),thksh3_var(n))
133 END DO
134 END DO
135 END IF
136 END IF
137
138 IF (iddw>0) CALL stoptimeg(ng)
139 GOTO 250
140 252 CONTINUE
141C-----------------------------------------------
142 ELSE
143 450 CONTINUE
144#include "lockon.inc"
145 IF(nsgdone>ngroup) THEN
146#include "lockoff.inc"
147 GOTO 452
148 ENDIF
149 ng=nsgdone
150 nsgdone = ng + 1
151#include "lockoff.inc"
152 IF(iparg(8,ng) == 1)GOTO 450
153 ity =iparg(5,ng)
154 IF(ity/=3.AND.ity/=7)GOTO 450
155
156 IF (iddw>0) CALL startimeg(ng)
157 mtn=iparg(1,ng)
158 IF(mtn == 0) GOTO 450
159 nel=iparg(2,ng)
160 nft=iparg(3,ng)
161 iad=iparg(4,ng)
162 ity=iparg(5,ng)
163 ihbe =iparg(23,ng)
164 ish3n =iparg(23,ng)
165 ithk = iparg(28,ng)
166 gbuf => elbuf_tab(ng)%GBUF
167c
168 IF (ity == 3) THEN
169 DO i=1,nel
170 n=nft+i
171
172 lev=sh4tree(3,n)
173 IF(lev < 0 .OR. gbuf%OFF(i)==zero) cycle
174
175 IF(ithk==0)THEN
176 thk=thke(n)
177 ELSE
178 thk=gbuf%THK(i)
179 END IF
180
181 lelt =1
182 nelt(1)=n
183
184 lelt1 =0
185 lelt2 =1
186
187 DO WHILE (lev < levelmax)
188
189 DO le=lelt1+1,lelt2
190
191 ne =nelt(le)
192 son=sh4tree(2,ne)
193
194 IF(son==0) cycle
195
196 lelt=lelt+1
197 nelt(lelt)=son
198
199 lelt=lelt+1
200 nelt(lelt)=son+1
201
202 lelt=lelt+1
203 nelt(lelt)=son+2
204
205 lelt=lelt+1
206 nelt(lelt)=son+3
207
208 END DO
209
210 lev =lev+1
211 lelt1 =lelt2
212 lelt2 =lelt
213
214 END DO
215
216 DO le=1,lelt
217 ne =nelt(le)
218 thksh4_var(ne)=thk
219 DO j=2,5
220 thknod(ixc(j,ne))=max(thknod(ixc(j,ne)),thk)
221 END DO
222 END DO
223 END DO
224 ELSEIF(ity == 7)THEN
225 DO i=1,nel
226 n=nft+i
227
228 lev=sh3tree(3,n)
229 IF(lev < 0 .OR. gbuf%OFF(i)==zero) cycle
230
231 IF(ithk==0)THEN
232 thk=thke(n)
233 ELSE
234 thk=gbuf%THK(i)
235 END IF
236
237 lelt =1
238 nelt(1)=n
239
240 lelt1 =0
241 lelt2 =1
242
243 DO WHILE (lev < levelmax)
244
245 DO le=lelt1+1,lelt2
246
247 ne =nelt(le)
248 son=sh3tree(2,ne)
249
250 IF(son==0) cycle
251
252 lelt=lelt+1
253 nelt(lelt)=son
254
255 lelt=lelt+1
256 nelt(lelt)=son+1
257
258 lelt=lelt+1
259 nelt(lelt)=son+2
260
261 lelt=lelt+1
262 nelt(lelt)=son+3
263
264 END DO
265
266 lev =lev+1
267 lelt1 =lelt2
268 lelt2 =lelt
269
270 END DO
271 DO le=1,lelt
272 ne =nelt(le)
273 thksh3_var(ne)=thk
274 DO j=2,4
275 thknod(ixtg(j,ne))=max(thknod(ixtg(j,ne)),thk)
276 END DO
277 END DO
278 END DO
279 END IF
280
281 IF (iddw>0) CALL stoptimeg(ng)
282 GOTO 450
283 452 CONTINUE
284 END IF
285C-----------------------------------------------
286 RETURN
287 END
my_real
#define my_real
Definition cppsort.cpp:32
startimeg
subroutine startimeg(ng)
Definition timer.F:1371
stoptimeg
subroutine stoptimeg(ng)
Definition timer.F:1419
max
#define max(a, b)
Definition macros.h:21
thickvar
subroutine thickvar(iparg, elbuf_tab, ixc, ixtg, thksh4_var, thksh3_var, thknod, thke, sh4tree, sh3tree)
Definition thickvar.F:36