OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s4kgeo3.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!|| s4kgeo3 ../engine/source/elements/solid/solide4/s4kgeo3.F
25!||--- called by ------------------------------------------------------
26!|| s4ke3 ../engine/source/elements/solid/solide4/s4ke3.F
27!||--- calls -----------------------------------------------------
28!|| s4keh3 ../engine/source/elements/solid/solide4/s4keh3.F
29!||====================================================================
30 SUBROUTINE s4kgeo3(
31 1 SIG, VOL, PX1, PX2,
32 2 PX3, PX4, PY1, PY2,
33 3 PY3, PY4, PZ1, PZ2,
34 4 PZ3, PZ4, K11, K12,
35 5 K13, K14, K22, K23,
36 6 K24, K33, K34, K44,
37 7 NEL)
38C-----------------------------------------------
39C I m p l i c i t T y p e s
40C-----------------------------------------------
41#include "implicit_f.inc"
42C-----------------------------------------------
43C G l o b a l P a r a m e t e r s
44C-----------------------------------------------
45#include "mvsiz_p.inc"
46C-----------------------------------------------
47C C o m m o n B l o c k s
48C-----------------------------------------------
49C-----------------------------------------------
50C D u m m y A r g u m e n t s
51C-----------------------------------------------
52 INTEGER NEL
53C REAL
54 my_real
55 . PX1(*), PX2(*), PX3(*), PX4(*),
56 . PY1(*), PY2(*), PY3(*), PY4(*),
57 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
58 . K11(3,3,*),K12(3,3,*),K13(3,3,*),K14(3,3,*) ,
59 . k22(3,3,*) ,k23(3,3,*) ,k24(3,3,*),
60 . k33(3,3,*),k34(3,3,*),
61 . k44(3,3,*),vol(*),sig(nel,6)
62C-----------------------------------------------
63C L o c a l V a r i a b l e s
64C-----------------------------------------------
65 INTEGER I,J
66C REAL
67 my_real
68 . FXX(MVSIZ),FYY(MVSIZ),FZZ(MVSIZ),FXY(MVSIZ),FXZ(MVSIZ),
69 . FYZ(MVSIZ)
70C-----------------------------------------------
71 DO I=1,nel
72 fxx(i)=sig(i,1)*vol(i)
73 fyy(i)=sig(i,2)*vol(i)
74 fzz(i)=sig(i,3)*vol(i)
75 fxy(i)=sig(i,4)*vol(i)
76 fyz(i)=sig(i,5)*vol(i)
77 fxz(i)=sig(i,6)*vol(i)
78 ENDDO
79C -----------------------------------------------
80 CALL s4keh3(1,nel,
81 . px1, py1 ,pz1 ,px1, py1, pz1,
82 . fxx, fyy ,fzz ,fxy, fyz, fxz,k11 )
83 CALL s4keh3(1,nel,
84 . px2, py2 ,pz2 ,px2, py2, pz2,
85 . fxx, fyy ,fzz ,fxy, fyz, fxz,k22 )
86 CALL s4keh3(1,nel,
87 . px3, py3 ,pz3 ,px3, py3, pz3,
88 . fxx, fyy ,fzz ,fxy, fyz, fxz,k33 )
89 CALL s4keh3(1,nel,
90 . px4, py4 ,pz4 ,px4, py4, pz4,
91 . fxx, fyy ,fzz ,fxy, fyz, fxz,k44 )
92 CALL s4keh3(1,nel,
93 . px1, py1 ,pz1 ,px2, py2, pz2,
94 . fxx, fyy ,fzz ,fxy, fyz, fxz,k12 )
95 CALL s4keh3(1,nel,
96 . px1, py1 ,pz1 ,px3, py3, pz3,
97 . fxx, fyy ,fzz ,fxy, fyz, fxz,k13 )
98 CALL s4keh3(1,nel,
99 . px1, py1 ,pz1 ,px4, py4, pz4,
100 . fxx, fyy ,fzz ,fxy, fyz, fxz,k14 )
101 CALL s4keh3(1,nel,
102 . px2, py2 ,pz2 ,px3, py3, pz3,
103 . fxx, fyy ,fzz ,fxy, fyz, fxz,k23 )
104 CALL s4keh3(1,nel,
105 . px2, py2 ,pz2 ,px4, py4, pz4,
106 . fxx, fyy ,fzz ,fxy, fyz, fxz,k24 )
107 CALL s4keh3(1,nel,
108 . px3, py3 ,pz3 ,px4, py4, pz4,
109 . fxx, fyy ,fzz ,fxy, fyz, fxz,k34 )
110C
111 RETURN
112 END
s4keh3
subroutine s4keh3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, fxx, fyy, fzz, fxy, fyz, fxz, kij)
Definition s4keh3.F:35
s4kgeo3
subroutine s4kgeo3(sig, vol, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, nel)
Definition s4kgeo3.F:38