OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s4coor3.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!|| s4coor3 ../engine/source/elements/solid/solide4/s4coor3.F
25!||--- called by ------------------------------------------------------
26!|| s4forc3 ../engine/source/elements/solid/solide4/s4forc3.F
27!||--- uses -----------------------------------------------------
28!|| element_mod ../common_source/modules/elements/element_mod.F90
29!||====================================================================
30 SUBROUTINE s4coor3(
31 1 X, IXS, V, W,
32 2 X1, X2, X3, X4,
33 3 Y1, Y2, Y3, Y4,
34 4 Z1, Z2, Z3, Z4,
35 5 VX1, VX2, VX3, VX4,
36 6 VY1, VY2, VY3, VY4,
37 7 VZ1, VZ2, VZ3, VZ4,
38 8 VDX1, VDX2, VDX3, VDX4,
39 9 VDY1, VDY2, VDY3, VDY4,
40 A VDZ1, VDZ2, VDZ3, VDZ4,
41 B VDX, VDY, VDZ, VD2,
42 C VIS, OFFG, OFF, SAV,
43 D RHO, RHOO, NC1, NC2,
44 E NC3, NC4, NGL, MXT,
45 F NGEO, F11, F21, F31,
46 G F12, F22, F32, F13,
47 H F23, F33, F14, F24,
48 I F34, XD1, XD2, XD3,
49 J XD4, YD1, YD2, YD3,
50 K YD4, ZD1, ZD2, ZD3,
51 L ZD4, XDP, NEL, JALE,
52 M ISMSTR, JEUL, JLAG)
53 use element_mod , only : nixs
54C-----------------------------------------------
55C I m p l i c i t T y p e s
56C-----------------------------------------------
57#include "implicit_f.inc"
58C-----------------------------------------------
59C G l o b a l P a r a m e t e r s
60C-----------------------------------------------
61#include "mvsiz_p.inc"
62C-----------------------------------------------
63C C o m m o n B l o c k s
64C-----------------------------------------------
65#include "scr05_c.inc"
66#include "scr18_c.inc"
67C-----------------------------------------------
68C D u m m y A r g u m e n t s
69C-----------------------------------------------
70 INTEGER, INTENT(IN) :: JALE
71 INTEGER, INTENT(IN) :: ISMSTR
72 INTEGER, INTENT(IN) :: JEUL
73 INTEGER, INTENT(IN) :: JLAG
74 INTEGER, INTENT(IN) :: NEL
75 INTEGER, INTENT(IN) :: IXS(NIXS,*)
76 my_real
77 . X(3,*),V(3,*),W(3,*), VIS(MVSIZ),
78 . X1(MVSIZ), X2(MVSIZ), X3(MVSIZ), X4(MVSIZ),
79 . Y1(MVSIZ), Y2(MVSIZ), Y3(MVSIZ), Y4(MVSIZ),
80 . Z1(MVSIZ), Z2(MVSIZ), Z3(MVSIZ), Z4(MVSIZ),
81 . VX1(MVSIZ), VX2(MVSIZ), VX3(MVSIZ), VX4(MVSIZ),
82 . VY1(MVSIZ), VY2(MVSIZ), VY3(MVSIZ), VY4(MVSIZ),
83 . VZ1(MVSIZ), VZ2(MVSIZ), VZ3(MVSIZ), VZ4(MVSIZ),
84 . VDX1(MVSIZ),VDX2(MVSIZ),VDX3(MVSIZ),VDX4(MVSIZ),
85 . VDY1(MVSIZ),VDY2(MVSIZ),VDY3(MVSIZ),VDY4(MVSIZ),
86 . VDZ1(MVSIZ),VDZ2(MVSIZ),VDZ3(MVSIZ),VDZ4(MVSIZ),
87 . VDX(MVSIZ), VDY(MVSIZ), VDZ(MVSIZ),VD2(MVSIZ),
88 . OFFG(NEL),OFF(MVSIZ),RHO(NEL),
89 . F11(MVSIZ),F21(MVSIZ),F31(MVSIZ),F12(MVSIZ),F22(MVSIZ),F32(MVSIZ),
90 . F13(MVSIZ),F23(MVSIZ),F33(MVSIZ),F14(MVSIZ),F24(MVSIZ),F34(MVSIZ),
91 . RHOO(MVSIZ)
92 INTEGER NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ), NC4(MVSIZ),
93 . mxt(mvsiz), ngl(mvsiz),ngeo(mvsiz)
94
95 DOUBLE PRECISION
96 . xdp(3,*), sav(nel,9),
97 . xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
98 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
99 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz)
100
101 my_real
102 . off_l
103C-----------------------------------------------
104C L o c a l V a r i a b l e s
105C-----------------------------------------------
106 INTEGER I
107C-----------------------------------------------
108#include "vectorize.inc"
109 DO I=1,nel
110 vis(i)=zero
111 ngeo(i)=ixs(10,i)
112 ngl(i)=ixs(11,i)
113 mxt(i)=ixs(1,i)
114 nc1(i)=ixs(2,i)
115 nc2(i)=ixs(4,i)
116 nc3(i)=ixs(7,i)
117 nc4(i)=ixs(6,i)
118 rhoo(i)=rho(i)
119 ENDDO
120#include "vectorize.inc"
121 DO i=1,nel
122 x1(i) =x(1,nc1(i))
123 y1(i) =x(2,nc1(i))
124 z1(i) =x(3,nc1(i))
125 x2(i) =x(1,nc2(i))
126 y2(i) =x(2,nc2(i))
127 z2(i) =x(3,nc2(i))
128 x3(i) =x(1,nc3(i))
129 y3(i) =x(2,nc3(i))
130 z3(i) =x(3,nc3(i))
131 x4(i) =x(1,nc4(i))
132 y4(i) =x(2,nc4(i))
133 z4(i) =x(3,nc4(i))
134 vx1(i)=v(1,nc1(i))
135 vy1(i)=v(2,nc1(i))
136 vz1(i)=v(3,nc1(i))
137 vx2(i)=v(1,nc2(i))
138 vy2(i)=v(2,nc2(i))
139 vz2(i)=v(3,nc2(i))
140 vx3(i)=v(1,nc3(i))
141 vy3(i)=v(2,nc3(i))
142 vz3(i)=v(3,nc3(i))
143 vx4(i)=v(1,nc4(i))
144 vy4(i)=v(2,nc4(i))
145 vz4(i)=v(3,nc4(i))
146 ENDDO
147
148 off_l = zero
149C----------------------------
150C NODAL COORDINATES |
151C----------------------------
152 IF((ismstr<=4.AND.jlag>0).OR.(ismstr==12.AND.idtmin(1)==3)) THEN
153
154 IF(iresp == 1) THEN
155#include "vectorize.inc"
156 DO i=1,nel
157 IF(abs(offg(i))>one)THEN
158 xd1(i)=sav(i,1)
159 yd1(i)=sav(i,2)
160 zd1(i)=sav(i,3)
161 xd2(i)=sav(i,4)
162 yd2(i)=sav(i,5)
163 zd2(i)=sav(i,6)
164 xd3(i)=sav(i,7)
165 yd3(i)=sav(i,8)
166 zd3(i)=sav(i,9)
167 xd4(i)=zero
168 yd4(i)=zero
169 zd4(i)=zero
170 off(i) = abs(offg(i))-one
171C OFF_L = MIN(OFF_L,OFFG(I))
172 ELSE
173 xd1(i)=xdp(1,nc1(i))
174 yd1(i)=xdp(2,nc1(i))
175 zd1(i)=xdp(3,nc1(i))
176 xd2(i)=xdp(1,nc2(i))
177 yd2(i)=xdp(2,nc2(i))
178 zd2(i)=xdp(3,nc2(i))
179 xd3(i)=xdp(1,nc3(i))
180 yd3(i)=xdp(2,nc3(i))
181 zd3(i)=xdp(3,nc3(i))
182 xd4(i)=xdp(1,nc4(i))
183 yd4(i)=xdp(2,nc4(i))
184 zd4(i)=xdp(3,nc4(i))
185 off(i) = abs(offg(i))
186C OFF_L = MIN(OFF_L,OFFG(I))
187 ENDIF
188 ENDDO
189 off_l = min(off_l,minval(offg(1:nel)))
190 ELSE
191#include "vectorize.inc"
192 DO i=1,nel
193 IF(abs(offg(i))>one)THEN
194 xd1(i)=sav(i,1)
195 yd1(i)=sav(i,2)
196 zd1(i)=sav(i,3)
197 xd2(i)=sav(i,4)
198 yd2(i)=sav(i,5)
199 zd2(i)=sav(i,6)
200 xd3(i)=sav(i,7)
201 yd3(i)=sav(i,8)
202 zd3(i)=sav(i,9)
203 xd4(i)=zero
204 yd4(i)=zero
205 zd4(i)=zero
206 off(i) = abs(offg(i))-one
207C OFF_L = MIN(OFF_L,OFFG(I))
208 ELSE
209 xd1(i)=x1(i)
210 yd1(i)=y1(i)
211 zd1(i)=z1(i)
212 xd2(i)=x2(i)
213 yd2(i)=y2(i)
214 zd2(i)=z2(i)
215 xd3(i)=x3(i)
216 yd3(i)=y3(i)
217 zd3(i)=z3(i)
218 xd4(i)=x4(i)
219 yd4(i)=y4(i)
220 zd4(i)=z4(i)
221 off(i) = abs(offg(i))
222C OFF_L = MIN(OFF_L,OFFG(I))
223 ENDIF
224 ENDDO
225 off_l = min(off_l,minval(offg(1:nel)))
226 ENDIF
227C
228 ELSE
229C
230 IF(iresp==1)THEN
231#include "vectorize.inc"
232 DO i=1,nel
233 xd1(i)=xdp(1,nc1(i))
234 yd1(i)=xdp(2,nc1(i))
235 zd1(i)=xdp(3,nc1(i))
236 xd2(i)=xdp(1,nc2(i))
237 yd2(i)=xdp(2,nc2(i))
238 zd2(i)=xdp(3,nc2(i))
239 xd3(i)=xdp(1,nc3(i))
240 yd3(i)=xdp(2,nc3(i))
241 zd3(i)=xdp(3,nc3(i))
242 xd4(i)=xdp(1,nc4(i))
243 yd4(i)=xdp(2,nc4(i))
244 zd4(i)=xdp(3,nc4(i))
245 off(i) = abs(offg(i))
246C OFF_L = MIN(OFF_L,OFFG(I))
247 ENDDO
248 off_l = min(off_l,minval(offg(1:nel)))
249 ELSE
250#include "vectorize.inc"
251 DO i=1,nel
252 xd1(i)=x1(i)
253 yd1(i)=y1(i)
254 zd1(i)=z1(i)
255 xd2(i)=x2(i)
256 yd2(i)=y2(i)
257 zd2(i)=z2(i)
258 xd3(i)=x3(i)
259 yd3(i)=y3(i)
260 zd3(i)=z3(i)
261 xd4(i)=x4(i)
262 yd4(i)=y4(i)
263 zd4(i)=z4(i)
264 off(i) = abs(offg(i))
265C OFF_L = MIN(OFF_L,OFFG(I))
266 ENDDO
267 off_l = min(off_l,minval(offg(1:nel)))
268 ENDIF
269C
270 ENDIF
271C
272 IF(off_l<zero)THEN
273#include "vectorize.inc"
274 DO i=1,nel
275 IF(offg(i)<zero)THEN
276 vx1(i)=zero
277 vy1(i)=zero
278 vz1(i)=zero
279 vx2(i)=zero
280 vy2(i)=zero
281 vz2(i)=zero
282 vx3(i)=zero
283 vy3(i)=zero
284 vz3(i)=zero
285 vx4(i)=zero
286 vy4(i)=zero
287 vz4(i)=zero
288 ENDIF
289 ENDDO
290 ENDIF
291
292 f11(1:nel)=zero
293 f21(1:nel)=zero
294 f31(1:nel)=zero
295 f12(1:nel)=zero
296 f22(1:nel)=zero
297 f32(1:nel)=zero
298 f13(1:nel)=zero
299 f23(1:nel)=zero
300 f33(1:nel)=zero
301 f14(1:nel)=zero
302 f24(1:nel)=zero
303 f34(1:nel)=zero
304C
305 IF (jlag/=0)THEN
306 vd2(1:nel)=zero
307 RETURN
308C
309 ELSEIF(jale/=0)THEN
310#include "vectorize.inc"
311 DO i=1,nel
312 vdx1(i)=vx1(i)-w(1,nc1(i))
313 vdy1(i)=vy1(i)-w(2,nc1(i))
314 vdz1(i)=vz1(i)-w(3,nc1(i))
315 vdx2(i)=vx2(i)-w(1,nc2(i))
316 vdy2(i)=vy2(i)-w(2,nc2(i))
317 vdz2(i)=vz2(i)-w(3,nc2(i))
318 vdx3(i)=vx3(i)-w(1,nc3(i))
319 vdy3(i)=vy3(i)-w(2,nc3(i))
320 vdz3(i)=vz3(i)-w(3,nc3(i))
321 vdx4(i)=vx4(i)-w(1,nc4(i))
322 vdy4(i)=vy4(i)-w(2,nc4(i))
323 vdz4(i)=vz4(i)-w(3,nc4(i))
324 ENDDO
325 ELSEIF(jeul/=0)THEN
326#include "vectorize.inc"
327 DO i=1,nel
328 vdx1(i)=vx1(i)
329 vdy1(i)=vy1(i)
330 vdz1(i)=vz1(i)
331 vdx2(i)=vx2(i)
332 vdy2(i)=vy2(i)
333 vdz2(i)=vz2(i)
334 vdx3(i)=vx3(i)
335 vdy3(i)=vy3(i)
336 vdz3(i)=vz3(i)
337 vdx4(i)=vx4(i)
338 vdy4(i)=vy4(i)
339 vdz4(i)=vz4(i)
340 ENDDO
341 ENDIF
342#include "vectorize.inc"
343 DO i=1,nel
344 vdx(i)=fourth*(vdx1(i)+vdx2(i)+vdx3(i)+vdx4(i))
345 vdy(i)=fourth*(vdy1(i)+vdy2(i)+vdy3(i)+vdy4(i))
346 vdz(i)=fourth*(vdz1(i)+vdz2(i)+vdz3(i)+vdz4(i))
347 vd2(i)=nine*(vdx(i)**2+vdy(i)**2+vdz(i)**2)
348 ENDDO
349C
350 RETURN
351 END
s4coor3
subroutine s4coor3(x, ixs, v, w, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, vdx1, vdx2, vdx3, vdx4, vdy1, vdy2, vdy3, vdy4, vdz1, vdz2, vdz3, vdz4, vdx, vdy, vdz, vd2, vis, offg, off, sav, rho, rhoo, nc1, nc2, nc3, nc4, ngl, mxt, ngeo, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, xd1, xd2, xd3, xd4, yd1, yd2, yd3, yd4, zd1, zd2, zd3, zd4, xdp, nel, jale, ismstr, jeul, jlag)
Definition s4coor3.F:53
min
#define min(a, b)
Definition macros.h:20