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 IS
64C-----------------------------------------------
65C -------symm(diag)------
66 IS =1
67 call s8zked3(1,nel,px1, py1 ,pz1 ,px1, py1, pz1, dd, k11,is)
68 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px2, py2, pz2, dd, k22,is)
69 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px3, py3, pz3, dd, k33,is)
70 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px4, py4, pz4, dd, k44,is)
71C
72 IF (iksup==0) THEN
73 CALL s8zkeg3(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
74 . py1, px1 ,pz1,px1 ,pz1,py1,gg,k11,is)
75 CALL s8zkeg3(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
76 . py2, px2 ,pz2,px2 ,pz2,py2,gg,k22,is)
77 CALL s8zkeg3(1,nel,py3, px3 ,pz3,px3 ,pz3,py3,
78 . py3, px3 ,pz3,px3 ,pz3,py3,gg,k33,is)
79 CALL s8zkeg3(1,nel,py4, px4 ,pz4,px4 ,pz4,py4,
80 . py4, px4 ,pz4,px4 ,pz4,py4,gg,k44,is)
81 ENDIF
82C -------non symm------
83 is =0
84 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px2, py2, pz2, dd, k12,is)
85 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px3, py3, pz3, dd, k13,is)
86 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px4, py4, pz4, dd, k14,is)
87 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px3, py3, pz3, dd, k23,is)
88 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px4, py4, pz4, dd, k24,is)
89 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px4, py4, pz4, dd, k34,is)
90C
91 IF (iksup==0) THEN
92 CALL s8zkeg3(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
93 . py2, px2 ,pz2,px2 ,pz2,py2,gg,k12,is)
94 CALL s8zkeg3(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
95 . py3, px3 ,pz3,px3 ,pz3,py3,gg,k13,is)
96 CALL s8zkeg3(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
97 . py4, px4 ,pz4,px4 ,pz4,py4,gg,k14,is)
98 CALL s8zkeg3(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
99 . py3, px3 ,pz3,px3 ,pz3,py3,gg,k23,is)
100 CALL s8zkeg3(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
101 . py4, px4 ,pz4,px4 ,pz4,py4,gg,k24,is)
102 CALL s8zkeg3(1,nel,py3, px3 ,pz3,px3 ,pz3,py3,
103 . py4, px4 ,pz4,px4 ,pz4,py4,gg,k34,is)
104 ELSE
105 CALL s8zkeg33(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
106 . py2, px2 ,pz2,px2 ,pz2,py2,g33,k12,is)
107 CALL s8zkeg33(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
108 . py3, px3 ,pz3,px3 ,pz3,py3,g33,k13,is)
109 CALL s8zkeg33(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
110 . py4, px4 ,pz4,px4 ,pz4,py4,g33,k14,is)
111 CALL s8zkeg33(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
112 . py3, px3 ,pz3,px3 ,pz3,py3,g33,k23,is)
113 CALL s8zkeg33(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
114 . py4, px4 ,pz4,px4 ,pz4,py4,g33,k24,is)
115 CALL s8zkeg33(1,nel,py3, px3 ,pz3,px3 ,pz3,py3,
116 . py4, px4 ,pz4,px4 ,pz4,py4,g33,k34,is)
117C
118 CALL s8zkebg1(1,nel,px1, py1 ,pz1 ,px2, py2, pz2, dg, k12,is)
119 CALL s8zkebg1(1,nel,px1, py1 ,pz1 ,px3, py3, pz3, dg, k13,is)
120 CALL s8zkebg1(1,nel,px1, py1 ,pz1 ,px4, py4, pz4, dg, k14,is)
121 CALL s8zkebg1(1,nel,px2, py2 ,pz2 ,px3, py3, pz3, dg, k23,is)
122 CALL s8zkebg1(1,nel,px2, py2 ,pz2 ,px4, py4, pz4, dg, k24,is)
123 CALL s8zkebg1(1,nel,px3, py3 ,pz3 ,px4, py4, pz4, dg, k34,is)
124 is =1
125 CALL s8zkebg1(1,nel,px1, py1 ,pz1 ,px1, py1, pz1, dg, k11,is)
126 CALL s8zkebg1(1,nel,px2, py2 ,pz2 ,px2, py2, pz2, dg, k22,is)
127 CALL s8zkebg1(1,nel,px3, py3 ,pz3 ,px3, py3, pz3, dg, k33,is)
128 CALL s8zkebg1(1,nel,px4, py4 ,pz4 ,px4, py4, pz4, dg, k44,is)
129C
130 CALL s8zkeg33(1,nel,py1, px1 ,pz1,px1 ,pz1,py1,
131 . py1, px1 ,pz1,px1 ,pz1,py1,g33,k11,is)
132 CALL s8zkeg33(1,nel,py2, px2 ,pz2,px2 ,pz2,py2,
133 . py2, px2 ,pz2,px2 ,pz2,py2,g33,k22,is)
134 CALL s8zkeg33(1,nel,py3, px3 ,pz3,px3 ,pz3,py3,
135 . py3, px3 ,pz3,px3 ,pz3,py3,g33,k33,is)
136 CALL s8zkeg33(1,nel,py4, px4 ,pz4,px4 ,pz4,py4,
137 . py4, px4 ,pz4,px4 ,pz4,py4,g33,k44,is)
138 ENDIF
139C
140 RETURN
141 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
s4ke3
subroutine s4ke3(pm, geo, ixs, x, gbuf, etag, iddl, ndof, k_diag, k_lt, iadk, jdik, nel, ipm, igeo, ikgeo, bufmat, nft, mtn, ismstr, jhbe, irep, isorth, iformdt)
Definition s4ke3.F:51
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