OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8zlkep3.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!|| s8zlkep3 ../engine/source/elements/solid/solide8z/s8zlkep3.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!|| s8zked3 ../engine/source/elements/solid/solide8z/s8zked3.F
30!||====================================================================
31 SUBROUTINE s8zlkep3(
32 1 PX1, PX2, PX3, PX4,
33 2 PY1, PY2, PY3, PY4,
34 3 PZ1, PZ2, PZ3, PZ4,
35 4 K11, K12, K13, K14,
36 5 K15, K16, K17, K18,
37 6 K22, K23, K24, K25,
38 7 K26, K27, K28, K33,
39 8 K34, K35, K36, K37,
40 9 K38, K44, K45, K46,
41 A K47, K48, K55, K56,
42 B K57, K58, K66, K67,
43 C K68, K77, K78, K88,
44 D HH, VOL, FAC, ICP,
45 E NEL)
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C G l o b a l P a r a m e t e r s
52C-----------------------------------------------
53#include "mvsiz_p.inc"
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57C-----------------------------------------------
58C D u m m y A r g u m e n t s
59C-----------------------------------------------
60 INTEGER, INTENT(IN) :: NEL
61 INTEGER ICP
62C REAL
63 my_real
64 . PX1(*), PX2(*), PX3(*), PX4(*),
65 . PY1(*), PY2(*), PY3(*), PY4(*),
66 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
67 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K15(9,*) ,
68 . K16(9,*) ,K17(9,*) ,K18(9,*) ,K22(9,*) ,K23(9,*) ,
69 . K24(9,*) ,K25(9,*) ,K26(9,*) ,K27(9,*) ,K28(9,*) ,
70 . K33(9,*) ,K34(9,*) ,K35(9,*) ,K36(9,*) ,K37(9,*) ,
71 . K38(9,*) ,K44(9,*) ,K45(9,*) ,K46(9,*) ,K47(9,*) ,
72 . K48(9,*) ,K55(9,*) ,K56(9,*) ,K57(9,*) ,K58(9,*) ,
73 . K66(9,*) ,K67(9,*) ,K68(9,*) ,K77(9,*) ,K78(9,*) ,
74 . k88(9,*) ,hh(2,*) ,vol(*) ,fac(*)
75C-----------------------------------------------
76C L o c a l V a r i a b l e s
77C-----------------------------------------------
78 INTEGER I, J,IS
79C REAL
80 my_real
81 . TT,TV,DD(3,3,MVSIZ),
82 . PX5(MVSIZ), PX6(MVSIZ), PX7(MVSIZ), PX8(MVSIZ),
83 . PY5(MVSIZ), PY6(MVSIZ), PY7(MVSIZ), PY8(MVSIZ),
84 . PZ5(MVSIZ), PZ6(MVSIZ), PZ7(MVSIZ), PZ8(MVSIZ)
85C-----------------------------------------------
86 IF (ICP==1) then
87 DO i=1,nel
88 tv=(hh(1,i)+two_third*hh(2,i))*vol(i)
89 tt=tv
90 dd(1,1,i)=tt
91 dd(2,2,i)=tt
92 dd(3,3,i)=tt
93 dd(1,2,i)=tv
94 dd(1,3,i)=tv
95 dd(2,1,i)=tv
96 dd(2,3,i)=tv
97 dd(3,1,i)=tv
98 dd(3,2,i)=tv
99 ENDDO
100 ELSE
101 DO i=1,nel
102 tv=(hh(1,i)*fac(i)+two_third*hh(2,i)*(one-fac(i)))*vol(i)
103 tt=tv
104 dd(1,1,i)=tt
105 dd(2,2,i)=tt
106 dd(3,3,i)=tt
107 dd(1,2,i)=tv
108 dd(1,3,i)=tv
109 dd(2,1,i)=tv
110 dd(2,3,i)=tv
111 dd(3,1,i)=tv
112 dd(3,2,i)=tv
113 ENDDO
114 END IF !(ICP==1) THEN
115C
116 DO i=1,nel
117 px5(i)=-px3(i)
118 py5(i)=-py3(i)
119 pz5(i)=-pz3(i)
120 px6(i)=-px4(i)
121 py6(i)=-py4(i)
122 pz6(i)=-pz4(i)
123 px7(i)=-px1(i)
124 py7(i)=-py1(i)
125 pz7(i)=-pz1(i)
126 px8(i)=-px2(i)
127 py8(i)=-py2(i)
128 pz8(i)=-pz2(i)
129 ENDDO
130C -------symm(diag)------
131 is =1
132C
133 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px1, py1, pz1, dd, k11,is)
134 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px2, py2, pz2, dd, k22,is)
135 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px3, py3, pz3, dd, k33,is)
136 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px4, py4, pz4, dd, k44,is)
137 CALL s8zked3(1,nel,px5, py5 ,pz5 ,px5, py5, pz5, dd, k55,is)
138 CALL s8zked3(1,nel,px6, py6 ,pz6 ,px6, py6, pz6, dd, k66,is)
139 CALL s8zked3(1,nel,px7, py7 ,pz7 ,px7, py7, pz7, dd, k77,is)
140 CALL s8zked3(1,nel,px8, py8 ,pz8 ,px8, py8, pz8, dd, k88,is)
141C -------non symm------
142 is =0
143C
144 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px2, py2, pz2, dd, k12,is)
145 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px3, py3, pz3, dd, k13,is)
146 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px4, py4, pz4, dd, k14,is)
147 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px5, py5, pz5, dd, k15,is)
148 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px6, py6, pz6, dd, k16,is)
149 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px7, py7, pz7, dd, k17,is)
150 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px8, py8, pz8, dd, k18,is)
151 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px3, py3, pz3, dd, k23,is)
152 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px4, py4, pz4, dd, k24,is)
153 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px5, py5, pz5, dd, k25,is)
154 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px6, py6, pz6, dd, k26,is)
155 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px7, py7, pz7, dd, k27,is)
156 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px8, py8, pz8, dd, k28,is)
157 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px4, py4, pz4, dd, k34,is)
158 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px5, py5, pz5, dd, k35,is)
159 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px6, py6, pz6, dd, k36,is)
160 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px7, py7, pz7, dd, k37,is)
161 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px8, py8, pz8, dd, k38,is)
162 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px5, py5, pz5, dd, k45,is)
163 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px6, py6, pz6, dd, k46,is)
164 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px7, py7, pz7, dd, k47,is)
165 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px8, py8, pz8, dd, k48,is)
166 CALL s8zked3(1,nel,px5, py5 ,pz5 ,px6, py6, pz6, dd, k56,is)
167 CALL s8zked3(1,nel,px5, py5 ,pz5 ,px7, py7, pz7, dd, k57,is)
168 CALL s8zked3(1,nel,px5, py5 ,pz5 ,px8, py8, pz8, dd, k58,is)
169 CALL s8zked3(1,nel,px6, py6 ,pz6 ,px7, py7, pz7, dd, k67,is)
170 CALL s8zked3(1,nel,px6, py6 ,pz6 ,px8, py8, pz8, dd, k68,is)
171 CALL s8zked3(1,nel,px7, py7 ,pz7 ,px8, py8, pz8, dd, k78,is)
172C
173C ---------------
174 RETURN
175 END
s8zked3
subroutine s8zked3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, dd, kij, is)
Definition s8zked3.F:34
s8zlkep3
subroutine s8zlkep3(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, hh, vol, fac, icp, nel)
Definition s8zlkep3.F:46