OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
scortho31.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!|| scortho31 ../engine/source/elements/thickshell/solidec/scortho31.F
25!||--- called by ------------------------------------------------------
26!|| srota6 ../engine/source/output/anim/generate/srota6.F
27!|| srotorth ../engine/source/elements/solid/srotorth.F
28!|| thsol ../engine/source/output/th/thsol.F
29!||====================================================================
30 SUBROUTINE scortho31(
31 . X1, X2, X3, X4, X5, X6, X7, X8,
32 . Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8,
33 . Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8,
34 . E1X, E2X, E3X, E1Y, E2Y, E3Y, E1Z, E2Z, E3Z )
35C-----------------------------------------------
36C I m p l i c i t T y p e s
37C-----------------------------------------------
38#include "implicit_f.inc"
39C-----------------------------------------------
40C D u m m y A r g u m e n t s
41C-----------------------------------------------
42C REAL
43 my_real
44 . x1, x2, x3, x4, x5, x6, x7, x8,
45 . y1, y2, y3, y4, y5, y6, y7, y8,
46 . z1, z2, z3, z4, z5, z6, z7, z8,
47 . e1x, e1y, e1z,
48 . e2x, e2y, e2z,
49 . e3x, e3y, e3z
50C-----------------------------------------------
51C L o c a l V a r i a b l e s
52C-----------------------------------------------
53 INTEGER I
54C REAL
55C 12
56 my_real
57 . x17 , x28 , x35 , x46,
58 . y17 , y28 , y35 , y46,
59 . z17 , z28 , z35 , z46,
60 . a17 , a28 ,
61 . b17 , b28 ,
62 . c17 , c28 ,
63 . rx , ry , rz ,
64 . tx , ty , tz ,
65 . det,c1,c2, sign,
66 . sx , sy , sz ,
67 . rsx , rsy , rsz
68C-----------------------------------------------
69 x17=x7-x1
70 x28=x8-x2
71 x35=x5-x3
72 x46=x6-x4
73 y17=y7-y1
74 y28=y8-y2
75 y35=y5-y3
76 y46=y6-y4
77 z17=z7-z1
78 z28=z8-z2
79 z35=z5-z3
80 z46=z6-z4
81C
82 rx=x17+x28-x35-x46
83 ry=y17+y28-y35-y46
84 rz=z17+z28-z35-z46
85 a17=x17+x46
86 a28=x28+x35
87 b17=y17+y46
88 b28=y28+y35
89 c17=z17+z46
90 c28=z28+z35
91C
92 tx=a17-a28
93 ty=b17-b28
94 tz=c17-c28
95C
96c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
97C
98 e3x = ty * rz - tz * ry
99 e3y = tz * rx - tx * rz
100 e3z = tx * ry - ty * rx
101C
102 det = sqrt(e3x*e3x + e3y*e3y + e3z*e3z)
103 IF ( det/=zero) det = one / det
104 e3x = e3x * det
105 e3y = e3y * det
106 e3z = e3z * det
107C
108 c1=sqrt(tx*tx+ty*ty+tz*tz)
109 c2=sqrt(rx*rx+ry*ry+rz*rz)
110 e1x=tx*c2 +(ry * e3z - rz * e3y)*c1
111 e1y=ty*c2 +(rz * e3x - rx * e3z)*c1
112 e1z=tz*c2 +(rx * e3y - ry * e3x)*c1
113 det = sqrt(e1x*e1x + e1y*e1y + e1z*e1z)
114 IF ( det/=zero) det = one / det
115 e1x = e1x*det
116 e1y = e1y*det
117 e1z = e1z*det
118C
119 e2x = e3y * e1z - e3z * e1y
120 e2y = e3z * e1x - e3x * e1z
121 e2z = e3x * e1y - e3y * e1x
122c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
123 RETURN
124 END
my_real
#define my_real
Definition cppsort.cpp:32
scortho31
subroutine scortho31(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition scortho31.F:35