OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s4cumg3.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!|| s4cumg3 ../engine/source/elements/solid/solide4/s4cumg3.f
25!||--- called by ------------------------------------------------------
26!|| s4ke3 ../engine/source/elements/solid/solide4/s4ke3.F
27!||--- calls -----------------------------------------------------
28!|| s8zkebg1 ../engine/source/elements/solid/solide8z/s8zkebg1.F
29!|| s8zked3 ../engine/source/elements/solid/solide8z/s8zked3.F
30!|| s8zkeg3 ../engine/source/elements/solid/solide8z/s8zkeg3.F
31!|| s8zkeg33 ../engine/source/elements/solid/solide8z/s8zkeg33.F
32!||====================================================================
33 SUBROUTINE s4cumg3(
34 1 PX1, PX2, PX3, PX4,
35 2 PY1, PY2, PY3, PY4,
36 3 PZ1, PZ2, PZ3, PZ4,
37 4 K11, K12, K13, K14,
38 5 K22, K23, K24, K33,
39 6 K34, K44, DD, GG,
40 7 DG, G33, IKSUP, NEL)
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45C-----------------------------------------------
46C C o m m o n B l o c k s
47C-----------------------------------------------
48C-----------------------------------------------
49C D u m m y A r g u m e n t s
50C-----------------------------------------------
51 INTEGER, INTENT(IN) :: NEL
52 INTEGER IKSUP
53 my_real
54 . PX1(*), PX2(*), PX3(*), PX4(*),
55 . PY1(*), PY2(*), PY3(*), PY4(*),
56 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
57 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K22(9,*) ,
58 . k23(9,*) ,k24(9,*) ,k33(9,*) ,k34(9,*) ,k44(9,*) ,
59 . dd(3,3,*),gg(*),dg(9,*) ,g33(9,*)
60C-----------------------------------------------
61C L o c a l V a r i a b l e s
62C-----------------------------------------------
63 INTEGER I, J,IS
64 my_real
65 . TT,TV
66C-----------------------------------------------
67C -------symm(diag)------
68 IS =1
69 call s8zked3(1,nel,px1, py1 ,pz1 ,px1, py1, pz1, dd, k11,is)
70 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px2, py2, pz2, dd, k22,is)
71 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px3, py3, pz3, dd, k33,is)
72 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px4, py4, pz4, dd, k44,is)
73C
74 IF (iksup==0) THEN
75 CALL s8zkeg3(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
76 . py1, px1 ,pz1,px1 ,pz1,py1,gg,k11,is)
77 CALL s8zkeg3(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
78 . py2, px2 ,pz2,px2 ,pz2,py2,gg,k22,is)
79 CALL s8zkeg3(1,nel,py3, px3 ,pz3,px3 ,pz3,py3,
80 . py3, px3 ,pz3,px3 ,pz3,py3,gg,k33,is)
81 CALL s8zkeg3(1,nel,py4, px4 ,pz4,px4 ,pz4,py4,
82 . py4, px4 ,pz4,px4 ,pz4,py4,gg,k44,is)
83 ENDIF
84C -------non symm------
85 is =0
86 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px2, py2, pz2, dd, k12,is)
87 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px3, py3, pz3, dd, k13,is)
88 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px4, py4, pz4, dd, k14,is)
89 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px3, py3, pz3, dd, k23,is)
90 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px4, py4, pz4, dd, k24,is)
91 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px4, py4, pz4, dd, k34,is)
92C
93 IF (iksup==0) THEN
94 CALL s8zkeg3(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
95 . py2, px2 ,pz2,px2 ,pz2,py2,gg,k12,is)
96 CALL s8zkeg3(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
97 . py3, px3 ,pz3,px3 ,pz3,py3,gg,k13,is)
98 CALL s8zkeg3(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
99 . py4, px4 ,pz4,px4 ,pz4,py4,gg,k14,is)
100 CALL s8zkeg3(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
101 . py3, px3 ,pz3,px3 ,pz3,py3,gg,k23,is)
102 CALL s8zkeg3(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
103 . py4, px4 ,pz4,px4 ,pz4,py4,gg,k24,is)
104 CALL s8zkeg3(1,nel,py3, px3 ,pz3,px3 ,pz3,py3,
105 . py4, px4 ,pz4,px4 ,pz4,py4,gg,k34,is)
106 ELSE
107 CALL s8zkeg33(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
108 . py2, px2 ,pz2,px2 ,pz2,py2,g33,k12,is)
109 CALL s8zkeg33(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
110 . py3, px3 ,pz3,px3 ,pz3,py3,g33,k13,is)
111 CALL s8zkeg33(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
112 . py4, px4 ,pz4,px4 ,pz4,py4,g33,k14,is)
113 CALL s8zkeg33(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
114 . py3, px3 ,pz3,px3 ,pz3,py3,g33,k23,is)
115 CALL s8zkeg33(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
116 . py4, px4 ,pz4,px4 ,pz4,py4,g33,k24,is)
117 CALL s8zkeg33(1,nel,py3, px3 ,pz3,px3 ,pz3,py3,
118 . py4, px4 ,pz4,px4 ,pz4,py4,g33,k34,is)
119C
120 CALL s8zkebg1(1,nel,px1, py1 ,pz1 ,px2, py2, pz2, dg, k12,is)
121 CALL s8zkebg1(1,nel,px1, py1 ,pz1 ,px3, py3, pz3, dg, k13,is)
122 CALL s8zkebg1(1,nel,px1, py1 ,pz1 ,px4, py4, pz4, dg, k14,is)
123 CALL s8zkebg1(1,nel,px2, py2 ,pz2 ,px3, py3, pz3, dg, k23,is)
124 CALL s8zkebg1(1,nel,px2, py2 ,pz2 ,px4, py4, pz4, dg, k24,is)
125 CALL s8zkebg1(1,nel,px3, py3 ,pz3 ,px4, py4, pz4, dg, k34,is)
126 is =1
127 CALL s8zkebg1(1,nel,px1, py1 ,pz1 ,px1, py1, pz1, dg, k11,is)
128 CALL s8zkebg1(1,nel,px2, py2 ,pz2 ,px2, py2, pz2, dg, k22,is)
129 CALL s8zkebg1(1,nel,px3, py3 ,pz3 ,px3, py3, pz3, dg, k33,is)
130 CALL s8zkebg1(1,nel,px4, py4 ,pz4 ,px4, py4, pz4, dg, k44,is)
131C
132 CALL s8zkeg33(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
133 . py1, px1 ,pz1,px1 ,pz1,py1,g33,k11,is)
134 CALL s8zkeg33(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
135 . py2, px2 ,pz2,px2 ,pz2,py2,g33,k22,is)
136 CALL s8zkeg33(1,nel,py3, px3 ,pz3,px3 ,pz3,py3,
137 . py3, px3 ,pz3,px3 ,pz3,py3,g33,k33,is)
138 CALL s8zkeg33(1,nel,py4, px4 ,pz4,px4 ,pz4,py4,
139 . py4, px4 ,pz4,px4 ,pz4,py4,g33,k44,is)
140 ENDIF
141C
142 RETURN
143 END
s4cumg3
subroutine s4cumg3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, dd, gg, dg, g33, iksup, nel)
Definition s4cumg3.F:41
s8zkebg1
subroutine s8zkebg1(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, dg, kij, is)
Definition s8zkebg1.F:32
s8zked3
subroutine s8zked3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, dd, kij, is)
Definition s8zked3.F:34
s8zkeg33
subroutine s8zkeg33(lft, llt, pxyi, pyxi, pxzi, pzxi, pyzi, pzyi, pxyj, pyxj, pxzj, pzxj, pyzj, pzyj, g33, kij, is)
Definition s8zkeg33.F:34
s8zkeg3
subroutine s8zkeg3(lft, llt, pxyi, pyxi, pxzi, pzxi, pyzi, pzyi, pxyj, pyxj, pxzj, pzxj, pyzj, pzyj, g, kij, is)
Definition s8zkeg3.F:34