OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s10kgeo3.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!|| s10kgeo3 ../engine/source/elements/solid/solide10/s10kgeo3.f
25!||--- called by ------------------------------------------------------
26!|| s10ke3 ../engine/source/elements/solid/solide10/s10ke3.F
27!||--- calls -----------------------------------------------------
28!|| s4keh3 ../engine/source/elements/solid/solide4/s4keh3.F
29!||====================================================================
30 SUBROUTINE s10kgeo3(
31 1 SIG, VOL, PX, PY,
32 2 PZ, K11, K12, K13,
33 3 K14, K15, K16, K17,
34 4 K18, K19, K10, K22,
35 5 K23, K24, K25, K26,
36 6 K27, K28, K29, K20,
37 7 K33, K34, K35, K36,
38 8 K37, K38, K39, K30,
39 9 K44, K45, K46, K47,
40 A K48, K49, K40, K55,
41 B K56, K57, K58, K59,
42 C K50, K66, K67, K68,
43 D K69, K60, K77, K78,
44 E K79, K70, K88, K89,
45 F K80, K99, K90, K00,
46 G NEL)
47C-----------------------------------------------
48C I m p l i c i t T y p e s
49C-----------------------------------------------
50#include "implicit_f.inc"
51C-----------------------------------------------
52C G l o b a l P a r a m e t e r s
53C-----------------------------------------------
54#include "mvsiz_p.inc"
55C-----------------------------------------------
56C C o m m o n B l o c k s
57C-----------------------------------------------
58C-----------------------------------------------
59C D u m m y A r g u m e n t s
60C-----------------------------------------------
61 INTEGER NEL
62C REAL
63 my_real
64 . PX(MVSIZ,10),PY(MVSIZ,10),PZ(MVSIZ,10),
65 . VOL(*),SIG(NEL,6)
66 my_real
67 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K15(9,*) ,
68 . K16(9,*) ,K17(9,*) ,K18(9,*) ,K19(9,*) ,K10(9,*) ,
69 . K22(9,*) ,K23(9,*) ,K24(9,*) ,K25(9,*) ,K26(9,*) ,
70 . K27(9,*) ,K28(9,*) ,K29(9,*) ,K20(9,*) ,K33(9,*) ,
71 . K34(9,*) ,K35(9,*) ,K36(9,*) ,K37(9,*) ,K38(9,*) ,
72 . K39(9,*) ,K30(9,*) ,K44(9,*) ,K45(9,*) ,K46(9,*) ,
73 . K47(9,*) ,K48(9,*) ,K49(9,*) ,K40(9,*) ,K55(9,*) ,
74 . K56(9,*) ,K57(9,*) ,K58(9,*) ,K59(9,*) ,K50(9,*) ,
75 . K66(9,*) ,K67(9,*) ,K68(9,*) ,K69(9,*) ,K60(9,*) ,
76 . K77(9,*) ,K78(9,*) ,K79(9,*) ,K70(9,*) ,K88(9,*) ,
77 . k89(9,*) ,k80(9,*) ,k99(9,*) ,k90(9,*) ,k00(9,*)
78C-----------------------------------------------
79C L o c a l V a r i a b l e s
80C-----------------------------------------------
81 INTEGER I,J
82 my_real
83 . fxx(mvsiz),fyy(mvsiz),fzz(mvsiz),fxy(mvsiz),fxz(mvsiz),
84 . fyz(mvsiz)
85C-----------------------------------------------
86 DO i=1,nel
87 fxx(i)=sig(i,1)*vol(i)
88 fyy(i)=sig(i,2)*vol(i)
89 fzz(i)=sig(i,3)*vol(i)
90 fxy(i)=sig(i,4)*vol(i)
91 fyz(i)=sig(i,5)*vol(i)
92 fxz(i)=sig(i,6)*vol(i)
93 ENDDO
94C -----------------------------------------------
95 CALL s4keh3(1,nel,
96 . px(1,1),py(1,1),pz(1,1),px(1,1),py(1,1),pz(1,1),
97 . fxx, fyy ,fzz ,fxy, fyz, fxz,k11 )
98 CALL s4keh3(1,nel,
99 . px(1,2),py(1,2),pz(1,2),px(1,2),py(1,2),pz(1,2),
100 . fxx, fyy ,fzz ,fxy, fyz, fxz,k22 )
101 CALL s4keh3(1,nel,
102 . px(1,3),py(1,3),pz(1,3),px(1,3),py(1,3),pz(1,3),
103 . fxx, fyy ,fzz ,fxy, fyz, fxz,k33 )
104 CALL s4keh3(1,nel,
105 . px(1,4),py(1,4),pz(1,4),px(1,4),py(1,4),pz(1,4),
106 . fxx, fyy ,fzz ,fxy, fyz, fxz,k44 )
107 CALL s4keh3(1,nel,
108 . px(1,5),py(1,5),pz(1,5),px(1,5),py(1,5),pz(1,5),
109 . fxx, fyy ,fzz ,fxy, fyz, fxz,k55 )
110 CALL s4keh3(1,nel,
111 . px(1,6),py(1,6),pz(1,6),px(1,6),py(1,6),pz(1,6),
112 . fxx, fyy ,fzz ,fxy, fyz, fxz,k66 )
113 CALL s4keh3(1,nel,
114 . px(1,7),py(1,7),pz(1,7),px(1,7),py(1,7),pz(1,7),
115 . fxx, fyy ,fzz ,fxy, fyz, fxz,k77 )
116 CALL s4keh3(1,nel,
117 . px(1,8),py(1,8),pz(1,8),px(1,8),py(1,8),pz(1,8),
118 . fxx, fyy ,fzz ,fxy, fyz, fxz,k88 )
119 CALL s4keh3(1,nel,
120 . px(1,9),py(1,9),pz(1,9),px(1,9),py(1,9),pz(1,9),
121 . fxx, fyy ,fzz ,fxy, fyz, fxz,k99 )
122 CALL s4keh3(1,nel,
123 . px(1,10),py(1,10),pz(1,10),px(1,10),py(1,10),pz(1,10),
124 . fxx, fyy ,fzz ,fxy, fyz, fxz,k00 )
125C
126 CALL s4keh3(1,nel,
127 . px(1,1),py(1,1),pz(1,1),px(1,2),py(1,2),pz(1,2),
128 . fxx, fyy ,fzz ,fxy, fyz, fxz,k12 )
129 CALL s4keh3(1,nel,
130 . px(1,1),py(1,1),pz(1,1),px(1,3),py(1,3),pz(1,3),
131 . fxx, fyy ,fzz ,fxy, fyz, fxz,k13 )
132 CALL s4keh3(1,nel,
133 . px(1,1),py(1,1),pz(1,1),px(1,4),py(1,4),pz(1,4),
134 . fxx, fyy ,fzz ,fxy, fyz, fxz,k14 )
135 CALL s4keh3(1,nel,
136 . px(1,1),py(1,1),pz(1,1),px(1,5),py(1,5),pz(1,5),
137 . fxx, fyy ,fzz ,fxy, fyz, fxz,k15 )
138 CALL s4keh3(1,nel,
139 . px(1,1),py(1,1),pz(1,1),px(1,6),py(1,6),pz(1,6),
140 . fxx, fyy ,fzz ,fxy, fyz, fxz,k16 )
141 CALL s4keh3(1,nel,
142 . px(1,1),py(1,1),pz(1,1),px(1,7),py(1,7),pz(1,7),
143 . fxx, fyy ,fzz ,fxy, fyz, fxz,k17 )
144 CALL s4keh3(1,nel,
145 . px(1,1),py(1,1),pz(1,1),px(1,8),py(1,8),pz(1,8),
146 . fxx, fyy ,fzz ,fxy, fyz, fxz,k18 )
147 CALL s4keh3(1,nel,
148 . px(1,1),py(1,1),pz(1,1),px(1,9),py(1,9),pz(1,9),
149 . fxx, fyy ,fzz ,fxy, fyz, fxz,k19 )
150 CALL s4keh3(1,nel,
151 . px(1,1),py(1,1),pz(1,1),px(1,10),py(1,10),pz(1,10),
152 . fxx, fyy ,fzz ,fxy, fyz, fxz,k10 )
153 CALL s4keh3(1,nel,
154 . px(1,2),py(1,2),pz(1,2),px(1,3),py(1,3),pz(1,3),
155 . fxx, fyy ,fzz ,fxy, fyz, fxz,k23 )
156 CALL s4keh3(1,nel,
157 . px(1,2),py(1,2),pz(1,2),px(1,4),py(1,4),pz(1,4),
158 . fxx, fyy ,fzz ,fxy, fyz, fxz,k24 )
159 CALL s4keh3(1,nel,
160 . px(1,2),py(1,2),pz(1,2),px(1,5),py(1,5),pz(1,5),
161 . fxx, fyy ,fzz ,fxy, fyz, fxz,k25 )
162 CALL s4keh3(1,nel,
163 . px(1,2),py(1,2),pz(1,2),px(1,6),py(1,6),pz(1,6),
164 . fxx, fyy ,fzz ,fxy, fyz, fxz,k26 )
165 CALL s4keh3(1,nel,
166 . px(1,2),py(1,2),pz(1,2),px(1,7),py(1,7),pz(1,7),
167 . fxx, fyy ,fzz ,fxy, fyz, fxz,k27 )
168 CALL s4keh3(1,nel,
169 . px(1,2),py(1,2),pz(1,2),px(1,8),py(1,8),pz(1,8),
170 . fxx, fyy ,fzz ,fxy, fyz, fxz,k28 )
171 CALL s4keh3(1,nel,
172 . px(1,2),py(1,2),pz(1,2),px(1,9),py(1,9),pz(1,9),
173 . fxx, fyy ,fzz ,fxy, fyz, fxz,k29 )
174 CALL s4keh3(1,nel,
175 . px(1,2),py(1,2),pz(1,2),px(1,10),py(1,10),pz(1,10),
176 . fxx, fyy ,fzz ,fxy, fyz, fxz,k20 )
177 CALL s4keh3(1,nel,
178 . px(1,3),py(1,3),pz(1,3),px(1,4),py(1,4),pz(1,4),
179 . fxx, fyy ,fzz ,fxy, fyz, fxz,k34 )
180 CALL s4keh3(1,nel,
181 . px(1,3),py(1,3),pz(1,3),px(1,5),py(1,5),pz(1,5),
182 . fxx, fyy ,fzz ,fxy, fyz, fxz,k35 )
183 CALL s4keh3(1,nel,
184 . px(1,3),py(1,3),pz(1,3),px(1,6),py(1,6),pz(1,6),
185 . fxx, fyy ,fzz ,fxy, fyz, fxz,k36 )
186 CALL s4keh3(1,nel,
187 . px(1,3),py(1,3),pz(1,3),px(1,7),py(1,7),pz(1,7),
188 . fxx, fyy ,fzz ,fxy, fyz, fxz,k37 )
189 CALL s4keh3(1,nel,
190 . px(1,3),py(1,3),pz(1,3),px(1,8),py(1,8),pz(1,8),
191 . fxx, fyy ,fzz ,fxy, fyz, fxz,k38 )
192 CALL s4keh3(1,nel,
193 . px(1,3),py(1,3),pz(1,3),px(1,9),py(1,9),pz(1,9),
194 . fxx, fyy ,fzz ,fxy, fyz, fxz,k39 )
195 CALL s4keh3(1,nel,
196 . px(1,3),py(1,3),pz(1,3),px(1,10),py(1,10),pz(1,10),
197 . fxx, fyy ,fzz ,fxy, fyz, fxz,k30 )
198 CALL s4keh3(1,nel,
199 . px(1,4),py(1,4),pz(1,4),px(1,5),py(1,5),pz(1,5),
200 . fxx, fyy ,fzz ,fxy, fyz, fxz,k45 )
201 CALL s4keh3(1,nel,
202 . px(1,4),py(1,4),pz(1,4),px(1,6),py(1,6),pz(1,6),
203 . fxx, fyy ,fzz ,fxy, fyz, fxz,k46 )
204 CALL s4keh3(1,nel,
205 . px(1,4),py(1,4),pz(1,4),px(1,7),py(1,7),pz(1,7),
206 . fxx, fyy ,fzz ,fxy, fyz, fxz,k47 )
207 CALL s4keh3(1,nel,
208 . px(1,4),py(1,4),pz(1,4),px(1,8),py(1,8),pz(1,8),
209 . fxx, fyy ,fzz ,fxy, fyz, fxz,k48 )
210 CALL s4keh3(1,nel,
211 . px(1,4),py(1,4),pz(1,4),px(1,9),py(1,9),pz(1,9),
212 . fxx, fyy ,fzz ,fxy, fyz, fxz,k49 )
213 CALL s4keh3(1,nel,
214 . px(1,4),py(1,4),pz(1,4),px(1,10),py(1,10),pz(1,10),
215 . fxx, fyy ,fzz ,fxy, fyz, fxz,k40 )
216 CALL s4keh3(1,nel,
217 . px(1,5),py(1,5),pz(1,5),px(1,6),py(1,6),pz(1,6),
218 . fxx, fyy ,fzz ,fxy, fyz, fxz,k56 )
219 CALL s4keh3(1,nel,
220 . px(1,5),py(1,5),pz(1,5),px(1,7),py(1,7),pz(1,7),
221 . fxx, fyy ,fzz ,fxy, fyz, fxz,k57 )
222 CALL s4keh3(1,nel,
223 . px(1,5),py(1,5),pz(1,5),px(1,8),py(1,8),pz(1,8),
224 . fxx, fyy ,fzz ,fxy, fyz, fxz,k58 )
225 CALL s4keh3(1,nel,
226 . px(1,5),py(1,5),pz(1,5),px(1,9),py(1,9),pz(1,9),
227 . fxx, fyy ,fzz ,fxy, fyz, fxz,k59 )
228 CALL s4keh3(1,nel,
229 . px(1,5),py(1,5),pz(1,5),px(1,10),py(1,10),pz(1,10),
230 . fxx, fyy ,fzz ,fxy, fyz, fxz,k50 )
231 CALL s4keh3(1,nel,
232 . px(1,6),py(1,6),pz(1,6),px(1,7),py(1,7),pz(1,7),
233 . fxx, fyy ,fzz ,fxy, fyz, fxz,k67 )
234 CALL s4keh3(1,nel,
235 . px(1,6),py(1,6),pz(1,6),px(1,8),py(1,8),pz(1,8),
236 . fxx, fyy ,fzz ,fxy, fyz, fxz,k68 )
237 CALL s4keh3(1,nel,
238 . px(1,6),py(1,6),pz(1,6),px(1,9),py(1,9),pz(1,9),
239 . fxx, fyy ,fzz ,fxy, fyz, fxz,k69 )
240 CALL s4keh3(1,nel,
241 . px(1,6),py(1,6),pz(1,6),px(1,10),py(1,10),pz(1,10),
242 . fxx, fyy ,fzz ,fxy, fyz, fxz,k60 )
243 CALL s4keh3(1,nel,
244 . px(1,7),py(1,7),pz(1,7),px(1,8),py(1,8),pz(1,8),
245 . fxx, fyy ,fzz ,fxy, fyz, fxz,k78 )
246 CALL s4keh3(1,nel,
247 . px(1,7),py(1,7),pz(1,7),px(1,9),py(1,9),pz(1,9),
248 . fxx, fyy ,fzz ,fxy, fyz, fxz,k79 )
249 CALL s4keh3(1,nel,
250 . px(1,7),py(1,7),pz(1,7),px(1,10),py(1,10),pz(1,10),
251 . fxx, fyy ,fzz ,fxy, fyz, fxz,k70 )
252 CALL s4keh3(1,nel,
253 . px(1,8),py(1,8),pz(1,8),px(1,9),py(1,9),pz(1,9),
254 . fxx, fyy ,fzz ,fxy, fyz, fxz,k89 )
255 CALL s4keh3(1,nel,
256 . px(1,8),py(1,8),pz(1,8),px(1,10),py(1,10),pz(1,10),
257 . fxx, fyy ,fzz ,fxy, fyz, fxz,k80 )
258 CALL s4keh3(1,nel,
259 . px(1,9),py(1,9),pz(1,9),px(1,10),py(1,10),pz(1,10),
260 . fxx, fyy ,fzz ,fxy, fyz, fxz,k90 )
261C
262 RETURN
263 END
s10kgeo3
subroutine s10kgeo3(sig, vol, px, py, pz, k11, k12, k13, k14, k15, k16, k17, k18, k19, k10, k22, k23, k24, k25, k26, k27, k28, k29, k20, k33, k34, k35, k36, k37, k38, k39, k30, k44, k45, k46, k47, k48, k49, k40, k55, k56, k57, k58, k59, k50, k66, k67, k68, k69, k60, k77, k78, k79, k70, k88, k89, k80, k99, k90, k00, nel)
Definition s10kgeo3.F:47
s4keh3
subroutine s4keh3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, fxx, fyy, fzz, fxy, fyz, fxz, kij)
Definition s4keh3.F:35