OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8xref_imp.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!||====================================================================
25!|| s8xref_imp ../engine/source/elements/solid/solide8s/s8xref_imp.F
26!||--- called by ------------------------------------------------------
27!|| s8sforc3 ../engine/source/elements/solid/solide8s/s8sforc3.F
28!|| srcoor3_imp ../engine/source/elements/solid/solide8s/srcoor3_imp.F
29!||====================================================================
30 SUBROUTINE s8xref_imp(
31 1 OFFG, XREF, XD1, XD2,
32 2 XD3, XD4, XD5, XD6,
33 3 XD7, XD8, YD1, YD2,
34 4 YD3, YD4, YD5, YD6,
35 5 YD7, YD8, ZD1, ZD2,
36 6 ZD3, ZD4, ZD5, ZD6,
37 7 ZD7, ZD8, R, NEL)
38C-----------------------------------------------
39C I m p l i c i t T y p e s
40C-----------------------------------------------
41#include "implicit_f.inc"
42C-----------------------------------------------
43C G l o b a l P a r a m e t e r s
44C-----------------------------------------------
45#include "mvsiz_p.inc"
46C-----------------------------------------------
47C C o m m o n B l o c k s
48C-----------------------------------------------
49C-----------------------------------------------
50C D u m m y A r g u m e n t s
51C-----------------------------------------------
52 INTEGER NEL
53C REAL
54 my_real
55 . OFFG(*)
56
57 DOUBLE PRECISION
58 . XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ),
59 . xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
60 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
61 . yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
62 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
63 . zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz),
64 . r(3,3,mvsiz), xref(nel,21)
65C-----------------------------------------------
66C L o c a l V a r i a b l e s
67C-----------------------------------------------
68 INTEGER I
69 DOUBLE PRECISION
70 . U,V,W,UL,VL,WL
71C=======================================================================
72 DO I=1,nel
73 !IF(ABS(OFFG(I)) <= ONE )THEN !sb voir comment faire avec OFFG
74 u = xd2(i)-xd1(i)
75 v = yd2(i)-yd1(i)
76 w = zd2(i)-zd1(i)
77 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
78 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
79 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
80 xref(i,1) = ul
81 xref(i,2) = vl
82 xref(i,3) = wl
83
84 u = xd3(i)-xd1(i)
85 v = yd3(i)-yd1(i)
86 w = zd3(i)-zd1(i)
87 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
88 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
89 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
90 xref(i,4) = ul
91 xref(i,5) = vl
92 xref(i,6) = wl
93
94 u = xd4(i)-xd1(i)
95 v = yd4(i)-yd1(i)
96 w = zd4(i)-zd1(i)
97 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
98 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
99 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
100 xref(i,7) = ul
101 xref(i,8) = vl
102 xref(i,9) = wl
103
104 u = xd5(i)-xd1(i)
105 v = yd5(i)-yd1(i)
106 w = zd5(i)-zd1(i)
107 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
108 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
109 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
110 xref(i,10) = ul
111 xref(i,11) = vl
112 xref(i,12) = wl
113
114 u = xd6(i)-xd1(i)
115 v = yd6(i)-yd1(i)
116 w = zd6(i)-zd1(i)
117 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
118 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
119 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
120 xref(i,13) = ul
121 xref(i,14) = vl
122 xref(i,15) = wl
123
124 u = xd7(i)-xd1(i)
125 v = yd7(i)-yd1(i)
126 w = zd7(i)-zd1(i)
127 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
128 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
129 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
130 xref(i,16) = ul
131 xref(i,17) = vl
132 xref(i,18) = wl
133
134 u = xd8(i)-xd1(i)
135 v = yd8(i)-yd1(i)
136 w = zd8(i)-zd1(i)
137 ul = r(1,1,i)*u+r(2,1,i)*v+r(3,1,i)*w
138 vl = r(1,2,i)*u+r(2,2,i)*v+r(3,2,i)*w
139 wl = r(1,3,i)*u+r(2,3,i)*v+r(3,3,i)*w
140 xref(i,19) = ul
141 xref(i,20) = vl
142 xref(i,21) = wl
143
144 ENDDO
145C-----------
146 RETURN
147 END
s8xref_imp
subroutine s8xref_imp(offg, xref, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, r, nel)
Definition s8xref_imp.F:38