OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8zkgei3.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!|| s8zkgei3 ../engine/source/elements/solid/solide8z/s8zkgei3.f
25!||--- called by ------------------------------------------------------
26!|| s8cke3 ../engine/source/elements/thickshell/solide8c/s8cke3.F
27!|| s8zke3 ../engine/source/elements/solid/solide8z/s8zke3.F
28!||--- calls -----------------------------------------------------
29!|| s8zkehij ../engine/source/elements/solid/solide8z/s8zkehij.F
30!||====================================================================
31 SUBROUTINE s8zkgei3(
32 1 SIG, VOL, PX1, PX2,
33 2 PX3, PX4, PY1, PY2,
34 3 PY3, PY4, PZ1, PZ2,
35 4 PZ3, PZ4, PX5, PX6,
36 5 PX7, PX8, PY5, PY6,
37 6 PY7, PY8, PZ5, PZ6,
38 7 PZ7, PZ8, K11, K12,
39 8 K13, K14, K15, K16,
40 9 K17, K18, K22, K23,
41 A K24, K25, K26, K27,
42 B K28, K33, K34, K35,
43 C K36, K37, K38, K44,
44 D K45, K46, K47, K48,
45 E K55, K56, K57, K58,
46 F K66, K67, K68, K77,
47 G K78, K88, NEL)
48C-----------------------------------------------
49C I m p l i c i t T y p e s
50C-----------------------------------------------
51#include "implicit_f.inc"
52C-----------------------------------------------
53C G l o b a l P a r a m e t e r s
54C-----------------------------------------------
55#include "mvsiz_p.inc"
56C-----------------------------------------------
57C C o m m o n B l o c k s
58C-----------------------------------------------
59C-----------------------------------------------
60C D u m m y A r g u m e n t s
61C-----------------------------------------------
62 INTEGER NEL
63C REAL
64 my_real
65 . PX1(*), PX2(*), PX3(*), PX4(*),
66 . PY1(*), PY2(*), PY3(*), PY4(*),
67 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
68 . PX5(*), PX6(*), PX7(*), PX8(*),
69 . PY5(*), PY6(*), PY7(*), PY8(*),
70 . PZ5(*), PZ6(*), PZ7(*), PZ8(*),
71 . K11(3,3,*),K12(3,3,*),K13(3,3,*),K14(3,3,*) ,K15(3,3,*) ,
72 . K16(3,3,*),K17(3,3,*),K18(3,3,*),K22(3,3,*) ,K23(3,3,*) ,
73 . K24(3,3,*),K25(3,3,*),K26(3,3,*),K27(3,3,*) ,K28(3,3,*) ,
74 . K33(3,3,*),K34(3,3,*),K35(3,3,*),K36(3,3,*) ,K37(3,3,*) ,
75 . K38(3,3,*),K44(3,3,*),K45(3,3,*),K46(3,3,*) ,K47(3,3,*) ,
76 . K48(3,3,*),K55(3,3,*),K56(3,3,*),K57(3,3,*) ,K58(3,3,*) ,
77 . K66(3,3,*),K67(3,3,*),K68(3,3,*),K77(3,3,*) ,K78(3,3,*) ,
78 . k88(3,3,*),vol(*),sig(nel,6)
79C-----------------------------------------------
80C L o c a l V a r i a b l e s
81C-----------------------------------------------
82 INTEGER I,J
83C REAL
84 my_real
85 . FXX(MVSIZ),FYY(MVSIZ),FZZ(MVSIZ),FXY(MVSIZ),FXZ(MVSIZ),
86 . FYZ(MVSIZ)
87C-----------------------------------------------
88 DO I=1,nel
89 fxx(i)=sig(i,1)*vol(i)
90 fyy(i)=sig(i,2)*vol(i)
91 fzz(i)=sig(i,3)*vol(i)
92 fxy(i)=sig(i,4)*vol(i)
93 fyz(i)=sig(i,5)*vol(i)
94 fxz(i)=sig(i,6)*vol(i)
95 ENDDO
96C------------K11,K17,K77
97 CALL s8zkehij(1,nel,
98 . px1, py1 ,pz1 ,px1, py1, pz1, k11 ,
99 . fxx, fyy ,fzz ,fxy, fyz, fxz)
100 CALL s8zkehij(1,nel,
101 . px1, py1 ,pz1 ,px7, py7, pz7, k17 ,
102 . fxx, fyy ,fzz ,fxy, fyz, fxz)
103 CALL s8zkehij(1,nel,
104 . px7, py7, pz7 ,px7, py7, pz7, k77 ,
105 . fxx, fyy ,fzz ,fxy, fyz, fxz)
106C------------K12,K18,K27,K78
107 CALL s8zkehij(1,nel,
108 . px1, py1 ,pz1 ,px2, py2, pz2, k12 ,
109 . fxx, fyy ,fzz ,fxy, fyz, fxz)
110 CALL s8zkehij(1,nel,
111 . px1, py1 ,pz1 ,px8, py8, pz8, k18 ,
112 . fxx, fyy ,fzz ,fxy, fyz, fxz)
113 CALL s8zkehij(1,nel,
114 . px2, py2 ,pz2 ,px7, py7, pz7, k27 ,
115 . fxx, fyy ,fzz ,fxy, fyz, fxz)
116 CALL s8zkehij(1,nel,
117 . px7, py7 ,pz7 ,px8, py8, pz8, k78 ,
118 . fxx, fyy ,fzz ,fxy, fyz, fxz)
119C------------K13,K15,K37,K57
120 CALL s8zkehij(1,nel,
121 . px1, py1 ,pz1 ,px3, py3, pz3, k13 ,
122 . fxx, fyy ,fzz ,fxy, fyz, fxz)
123 CALL s8zkehij(1,nel,
124 . px1, py1 ,pz1 ,px5, py5, pz5, k15 ,
125 . fxx, fyy ,fzz ,fxy, fyz, fxz)
126 CALL s8zkehij(1,nel,
127 . px3, py3 ,pz3 ,px7, py7, pz7, k37 ,
128 . fxx, fyy ,fzz ,fxy, fyz, fxz)
129 CALL s8zkehij(1,nel,
130 . px5, py5 ,pz5 ,px7, py7, pz7, k57 ,
131 . fxx, fyy ,fzz ,fxy, fyz, fxz)
132C------------K14,K16,K47,K67
133 CALL s8zkehij(1,nel,
134 . px1, py1 ,pz1 ,px4, py4, pz4, k14 ,
135 . fxx, fyy ,fzz ,fxy, fyz, fxz)
136 CALL s8zkehij(1,nel,
137 . px1, py1 ,pz1 ,px6, py6, pz6, k16 ,
138 . fxx, fyy ,fzz ,fxy, fyz, fxz)
139 CALL s8zkehij(1,nel,
140 . px4, py4 ,pz4 ,px7, py7, pz7, k47 ,
141 . fxx, fyy ,fzz ,fxy, fyz, fxz)
142 CALL s8zkehij(1,nel,
143 . px6, py6 ,pz6 ,px7, py7, pz7, k67 ,
144 . fxx, fyy ,fzz ,fxy, fyz, fxz)
145C------------K22,K28,K88
146 CALL s8zkehij(1,nel,
147 . px2, py2 ,pz2 ,px2, py2, pz2, k22 ,
148 . fxx, fyy ,fzz ,fxy, fyz, fxz)
149 CALL s8zkehij(1,nel,
150 . px2, py2 ,pz2 ,px8, py8, pz8, k28 ,
151 . fxx, fyy ,fzz ,fxy, fyz, fxz)
152 CALL s8zkehij(1,nel,
153 . px8, py8 ,pz8 ,px8, py8, pz8, k88 ,
154 . fxx, fyy ,fzz ,fxy, fyz, fxz)
155C------------K23,K25,K38,K58
156 CALL s8zkehij(1,nel,
157 . px2, py2 ,pz2 ,px3, py3, pz3, k23 ,
158 . fxx, fyy ,fzz ,fxy, fyz, fxz)
159 CALL s8zkehij(1,nel,
160 . px2, py2 ,pz2 ,px5, py5, pz5, k25 ,
161 . fxx, fyy ,fzz ,fxy, fyz, fxz)
162 CALL s8zkehij(1,nel,
163 . px3, py3 ,pz3 ,px8, py8, pz8, k38 ,
164 . fxx, fyy ,fzz ,fxy, fyz, fxz)
165 CALL s8zkehij(1,nel,
166 . px5, py5 ,pz5 ,px8, py8, pz8, k58 ,
167 . fxx, fyy ,fzz ,fxy, fyz, fxz)
168C------------K24,K26,K48,K68
169 CALL s8zkehij(1,nel,
170 . px2, py2 ,pz2 ,px4, py4, pz4, k24 ,
171 . fxx, fyy ,fzz ,fxy, fyz, fxz)
172 CALL s8zkehij(1,nel,
173 . px2, py2 ,pz2 ,px6, py6, pz6, k26 ,
174 . fxx, fyy ,fzz ,fxy, fyz, fxz)
175 CALL s8zkehij(1,nel,
176 . px4, py4 ,pz4 ,px8, py8, pz8, k48 ,
177 . fxx, fyy ,fzz ,fxy, fyz, fxz)
178 CALL s8zkehij(1,nel,
179 . px6, py6 ,pz6 ,px8, py8, pz8, k68 ,
180 . fxx, fyy ,fzz ,fxy, fyz, fxz)
181C------------K33,K35,K55
182 CALL s8zkehij(1,nel,
183 . px3, py3 ,pz3 ,px3, py3, pz3, k33 ,
184 . fxx, fyy ,fzz ,fxy, fyz, fxz)
185 CALL s8zkehij(1,nel,
186 . px3, py3 ,pz3 ,px5, py5, pz5, k35 ,
187 . fxx, fyy ,fzz ,fxy, fyz, fxz)
188 CALL s8zkehij(1,nel,
189 . px5, py5 ,pz5 ,px5, py5, pz5, k55 ,
190 . fxx, fyy ,fzz ,fxy, fyz, fxz)
191C------------K34,K36,K56,K45
192 CALL s8zkehij(1,nel,
193 . px3, py3 ,pz3 ,px4, py4, pz4, k34 ,
194 . fxx, fyy ,fzz ,fxy, fyz, fxz)
195 CALL s8zkehij(1,nel,
196 . px3, py3 ,pz3 ,px6, py6, pz6, k36 ,
197 . fxx, fyy ,fzz ,fxy, fyz, fxz)
198 CALL s8zkehij(1,nel,
199 . px5, py5 ,pz5 ,px6, py6, pz6, k56 ,
200 . fxx, fyy ,fzz ,fxy, fyz, fxz)
201 CALL s8zkehij(1,nel,
202 . px4, py4 ,pz4 ,px5, py5, pz5, k45 ,
203 . fxx, fyy ,fzz ,fxy, fyz, fxz)
204C------------K44,K46,K66
205 CALL s8zkehij(1,nel,
206 . px4, py4 ,pz4 ,px4, py4, pz4, k44 ,
207 . fxx, fyy ,fzz ,fxy, fyz, fxz)
208 CALL s8zkehij(1,nel,
209 . px4, py4 ,pz4 ,px6, py6, pz6, k46 ,
210 . fxx, fyy ,fzz ,fxy, fyz, fxz)
211 CALL s8zkehij(1,nel,
212 . px6, py6 ,pz6 ,px6, py6, pz6, k66 ,
213 . fxx, fyy ,fzz ,fxy, fyz, fxz)
214C
215 RETURN
216 END
s8zkehij
subroutine s8zkehij(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, kij, fxx, fyy, fzz, fxy, fyz, fxz)
Definition s8zkehij.F:31
s8zkgei3
subroutine s8zkgei3(sig, vol, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px5, px6, px7, px8, py5, py6, py7, py8, pz5, pz6, pz7, pz8, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, nel)
Definition s8zkgei3.F:48