OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
srepiso3.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!|| srepiso3 ../engine/source/elements/solid/solide/srepiso3.F
25!||--- called by ------------------------------------------------------
26!|| scoor3 ../engine/source/elements/solid/solide/scoor3.F
27!|| scoor3_fvm ../engine/source/ale/alefvm/scoor3_fvm.F
28!|| srcoor3 ../engine/source/elements/solid/solide/srcoor3.F
29!||====================================================================
30 SUBROUTINE srepiso3(
31 1 X1, X2, X3, X4,
32 2 X5, X6, X7, X8,
33 3 Y1, Y2, Y3, Y4,
34 4 Y5, Y6, Y7, Y8,
35 5 Z1, Z2, Z3, Z4,
36 6 Z5, Z6, Z7, Z8,
37 7 RX, RY, RZ, SX,
38 8 SY, SZ, TX, TY,
39 9 TZ, NEL)
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C G l o b a l P a r a m e t e r s
46C-----------------------------------------------
47#include "mvsiz_p.inc"
48C-----------------------------------------------
49C D u m m y A r g u m e n t s
50C-----------------------------------------------
51 INTEGER, INTENT(IN) :: NEL
52C REAL
53 DOUBLE PRECISION, INTENT(IN), DIMENSION(MVSIZ) ::
54 . X1, X2, X3, X4, X5, X6, X7, X8,
55 . Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8,
56 . Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8
57
58 my_real, INTENT(OUT), DIMENSION(MVSIZ) ::
59 . RX, RY, RZ, SX, SY, SZ, TX, TY, TZ
60C-----------------------------------------------
61C C o m m o n B l o c k s
62C-----------------------------------------------
63C-----------------------------------------------
64C L o c a l V a r i a b l e s
65C-----------------------------------------------
66 INTEGER I
67C REAL
68 my_real
69 . x17 , x28 , x35 , x46,
70 . y17 , y28 , y35 , y46,
71 . z17 , z28 , z35 , z46,
72 . a17 , a28 , b17 , b28 ,
73 . c17 , c28
74C=======================================================================
75 DO i=1,nel
76 x17 = x7(i) - x1(i)
77 x28 = x8(i) - x2(i)
78 x35 = x5(i) - x3(i)
79 x46 = x6(i) - x4(i)
80 y17 = y7(i) - y1(i)
81 y28 = y8(i) - y2(i)
82 y35 = y5(i) - y3(i)
83 y46 = y6(i) - y4(i)
84 z17 = z7(i) - z1(i)
85 z28 = z8(i) - z2(i)
86 z35 = z5(i) - z3(i)
87 z46 = z6(i) - z4(i)
88
89 a17 = x17 + x46
90 a28 = x28 + x35
91 b17 = y17 + y46
92 b28 = y28 + y35
93 c17 = z17 + z46
94 c28 = z28 + z35
95C
96 rx(i) = x17 + x28 - x35 - x46
97 ry(i) = y17 + y28 - y35 - y46
98 rz(i) = z17 + z28 - z35 - z46
99 sx(i) = a17 + a28
100 sy(i) = b17 + b28
101 sz(i) = c17 + c28
102 tx(i) = a17 - a28
103 ty(i) = b17 - b28
104 tz(i) = c17 - c28
105 ENDDO
106C-------------
107 RETURN
108 END
subroutine srepiso3(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, nel)
Definition srepiso3.F:40