OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
tetra4v.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!|| tetra4v ../engine/source/elements/solid/solide/tetra4v.F
25!||--- called by ------------------------------------------------------
26!|| deges4v ../engine/source/elements/solid/solide/deges4v.f
27!||====================================================================
28 SUBROUTINE tetra4v(DET,
29 . X1, X2, X3, X4,
30 . Y1, Y2, Y3, Y4,
31 . Z1, Z2, Z3, Z4)
32C-----------------------------------------------
33C I m p l i c i t T y p e s
34C-----------------------------------------------
35#include "implicit_f.inc"
36C-----------------------------------------------
37C D u m m y A r g u m e n t s
38C-----------------------------------------------
39 my_real
40 . det,
41 . x1, x2, x3, x4, y1, y2, y3, y4,z1, z2, z3, z4
42C-----------------------------------------------
43C L o c a l V a r i a b l e s
44C-----------------------------------------------
45 INTEGER I
46C 12
47 my_real
48 . x41, y41, z41, x42, y42, z42, x43, y43, z43
49 my_real
50 . rx, ry, rz, sx, sy, sz,tx, ty, tz,
51 . tsx, tsy,tsz
52C-----------------------------------------------
53 x43 = x4 - x3
54 y43 = y4 - y3
55 z43 = z4 - z3
56 x41 = x4 - x1
57 y41 = y4 - y1
58 z41 = z4 - z1
59 x42 = x4 - x2
60 y42 = y4 - y2
61 z42 = z4 - z2
62C
63 rx = -x41
64 ry = -y41
65 rz = -z41
66C
67 sx = -x42
68 sy = -y42
69 sz = -z42
70C
71 tx = -x43
72 ty = -y43
73 tz = -z43
74C
75 tsx = y43*z42 - y42*z43
76 tsy = z43*x42 - z42*x43
77 tsz = x43*y42 - x42*y43
78C------V=r*(s^t)/6
79 det = abs(rx*tsx + ry*tsy + rz*tsz)*one_over_6
80C
81 RETURN
82 END
my_real
#define my_real
Definition cppsort.cpp:32
deges4v
subroutine deges4v(det, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8)
Definition deges4v.F:37
tetra4v
subroutine tetra4v(det, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)
Definition tetra4v.F:32