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