OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
assem_s8.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!|| assem_s8 ../engine/source/implicit/assem_s8.F
25!||--- called by ------------------------------------------------------
26!|| s8cke3 ../engine/source/elements/thickshell/solide8c/s8cke3.F
27!|| s8ske3 ../engine/source/elements/solid/solide8s/s8ske3.F
28!|| s8zke3 ../engine/source/elements/solid/solide8z/s8zke3.F
29!||--- calls -----------------------------------------------------
30!|| assem_kii ../engine/source/implicit/imp_glob_k.F
31!|| assem_kij ../engine/source/implicit/imp_glob_k.F
32!||--- uses -----------------------------------------------------
33!|| element_mod ../common_source/modules/elements/element_mod.F90
34!||====================================================================
35 SUBROUTINE assem_s8(
36 1 IXS ,NEL ,IDDL ,NDOF ,K_DIAG,
37 2 K_LT ,IADK ,JDIK ,K11 ,K12 ,
38 3 K13 ,K14 ,K15 ,K16 ,K17 ,
39 4 K18 ,K22 ,K23 ,K24 ,K25 ,
40 5 K26 ,K27 ,K28 ,K33 ,K34 ,
41 6 K35 ,K36 ,K37 ,K38 ,K44 ,
42 7 K45 ,K46 ,K47 ,K48 ,K55 ,
43 8 K56 ,K57 ,K58 ,K66 ,K67 ,
44 9 K68 ,K77 ,K78 ,K88 ,OFF )
45 use element_mod , only : nixs
46C----6---------------------------------------------------------------7---------8
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 IXS(NIXS,*),NEL ,IDDL(*) ,NDOF(*) ,
61 . IADK(*) ,JDIK(*)
62C REAL
63 my_real
64 . K_DIAG(*) ,K_LT(*) ,
65 . K11(3,3,*),K12(3,3,*),K13(3,3,*),K14(3,3,*),
66 . K15(3,3,*),K16(3,3,*),K17(3,3,*),K18(3,3,*),
67 . K22(3,3,*),K23(3,3,*),K24(3,3,*),K25(3,3,*),
68 . K26(3,3,*),K27(3,3,*),K28(3,3,*),
69 . k33(3,3,*),k34(3,3,*),k35(3,3,*),k36(3,3,*),
70 . k37(3,3,*),k38(3,3,*),k44(3,3,*),k45(3,3,*),
71 . k46(3,3,*),k47(3,3,*),k48(3,3,*),
72 . k55(3,3,*),k56(3,3,*),k57(3,3,*),k58(3,3,*),
73 . k66(3,3,*),k67(3,3,*),k68(3,3,*),
74 . k77(3,3,*),k78(3,3,*),k88(3,3,*),off(*)
75C-----------------------------------------------
76C L o c a l V a r i a b l e s
77C-----------------------------------------------
78 INTEGER K,EP,NN(MVSIZ,8)
79C----6---------------------------------------------------------------7---------8
80 DO K = 1,8
81 do ep = 1,nel
82 nn(ep,k)=ixs(k+1,ep)
83 ENDDO
84 ENDDO
85C-------block diagonal k11,k22...k88-----
86 CALL assem_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt ,k11,3,off)
87 CALL assem_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt ,k22,3,off)
88 CALL assem_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt ,k33,3,off)
89 CALL assem_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt ,k44,3,off)
90 CALL assem_kii(nn(1,5),nel,iddl,iadk,k_diag,k_lt ,k55,3,off)
91 CALL assem_kii(nn(1,6),nel,iddl,iadk,k_diag,k_lt ,k66,3,off)
92 CALL assem_kii(nn(1,7),nel,iddl,iadk,k_diag,k_lt ,k77,3,off)
93 CALL assem_kii(nn(1,8),nel,iddl,iadk,k_diag,k_lt ,k88,3,off)
94C-------non diag kij-----
95 CALL assem_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
96 . k_diag,k_lt,k12,3,off)
97 CALL assem_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
98 . k_diag,k_lt,k13,3,off)
99 CALL assem_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
100 . k_diag,k_lt,k14,3,off)
101 CALL assem_kij(nn(1,1),nn(1,5),nel,iddl,iadk,jdik,
102 . k_diag,k_lt,k15,3,off)
103 CALL assem_kij(nn(1,1),nn(1,6),nel,iddl,iadk,jdik,
104 . k_diag,k_lt,k16,3,off)
105 CALL assem_kij(nn(1,1),nn(1,7),nel,iddl,iadk,jdik,
106 . k_diag,k_lt,k17,3,off)
107 CALL assem_kij(nn(1,1),nn(1,8),nel,iddl,iadk,jdik,
108 . k_diag,k_lt,k18,3,off)
109 CALL assem_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
110 . k_diag,k_lt,k23,3,off)
111 CALL assem_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
112 . k_diag,k_lt,k24,3,off)
113 CALL assem_kij(nn(1,2),nn(1,5),nel,iddl,iadk,jdik,
114 . k_diag,k_lt,k25,3,off)
115 CALL assem_kij(nn(1,2),nn(1,6),nel,iddl,iadk,jdik,
116 . k_diag,k_lt,k26,3,off)
117 CALL assem_kij(nn(1,2),nn(1,7),nel,iddl,iadk,jdik,
118 . k_diag,k_lt,k27,3,off)
119 CALL assem_kij(nn(1,2),nn(1,8),nel,iddl,iadk,jdik,
120 . k_diag,k_lt,k28,3,off)
121 CALL assem_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
122 . k_diag,k_lt,k34,3,off)
123 CALL assem_kij(nn(1,3),nn(1,5),nel,iddl,iadk,jdik,
124 . k_diag,k_lt,k35,3,off)
125 CALL assem_kij(nn(1,3),nn(1,6),nel,iddl,iadk,jdik,
126 . k_diag,k_lt,k36,3,off)
127 CALL assem_kij(nn(1,3),nn(1,7),nel,iddl,iadk,jdik,
128 . k_diag,k_lt,k37,3,off)
129 CALL assem_kij(nn(1,3),nn(1,8),nel,iddl,iadk,jdik,
130 . k_diag,k_lt,k38,3,off)
131 CALL assem_kij(nn(1,4),nn(1,5),nel,iddl,iadk,jdik,
132 . k_diag,k_lt,k45,3,off)
133 CALL assem_kij(nn(1,4),nn(1,6),nel,iddl,iadk,jdik,
134 . k_diag,k_lt,k46,3,off)
135 CALL assem_kij(nn(1,4),nn(1,7),nel,iddl,iadk,jdik,
136 . k_diag,k_lt,k47,3,off)
137 CALL assem_kij(nn(1,4),nn(1,8),nel,iddl,iadk,jdik,
138 . k_diag,k_lt,k48,3,off)
139 CALL assem_kij(nn(1,5),nn(1,6),nel,iddl,iadk,jdik,
140 . k_diag,k_lt,k56,3,off)
141 CALL assem_kij(nn(1,5),nn(1,7),nel,iddl,iadk,jdik,
142 . k_diag,k_lt,k57,3,off)
143 CALL assem_kij(nn(1,5),nn(1,8),nel,iddl,iadk,jdik,
144 . k_diag,k_lt,k58,3,off)
145 CALL assem_kij(nn(1,6),nn(1,7),nel,iddl,iadk,jdik,
146 . k_diag,k_lt,k67,3,off)
147 CALL assem_kij(nn(1,6),nn(1,8),nel,iddl,iadk,jdik,
148 . k_diag,k_lt,k68,3,off)
149 CALL assem_kij(nn(1,7),nn(1,8),nel,iddl,iadk,jdik,
150 . k_diag, k_lt,k78,3,off)
151C
152C----6---------------------------------------------------------------7---------8
153 RETURN
154 END
assem_s8
subroutine assem_s8(ixs, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, 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, off)
Definition assem_s8.F:45
assem_kii
subroutine assem_kii(ni, nel, iddl, iadk, k_diag, k_lt, kii, nd, off)
Definition imp_glob_k.F:964
assem_kij
subroutine assem_kij(ni, nj, nel, iddl, iadk, jdik, k_diag, k_lt, kij, nd, off)
Definition imp_glob_k.F:1043