OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
assem_s20.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_s20 ../engine/source/implicit/assem_s20.F
25!||--- called by ------------------------------------------------------
26!|| s20ke3 ../engine/source/elements/solid/solide20/s20ke3.F
27!||--- calls -----------------------------------------------------
28!|| assem_kii ../engine/source/implicit/imp_glob_k.F
29!|| assem_kij ../engine/source/implicit/imp_glob_k.F
30!||====================================================================
31 SUBROUTINE assem_s20(
32 1 IXS ,IXS20 , NEL ,IDDL ,NDOF ,K_DIAG,
33 2 K_LT ,IADK ,JDIK ,OFF ,
34 . K11 ,K12 ,K13 ,K14 ,K15 ,K16 ,K17 ,K18 ,K19 ,K110 ,
35 . K111 ,K112 ,K113 ,K114 ,K115 ,K116 ,K117 ,K118 ,K119 ,K120 ,
36 . K22 ,K23 ,K24 ,K25 ,K26 ,K27 ,K28 ,K29 ,K210 ,K211 ,
37 . K212 ,K213 ,K214 ,K215 ,K216 ,K217 ,K218 ,K219 ,K220 ,K33 ,
38 . K34 ,K35 ,K36 ,K37 ,K38 ,K39 ,K310 ,K311 ,K312 ,K313 ,
39 . K314 ,K315 ,K316 ,K317 ,K318 ,K319 ,K320 ,K44 ,K45 ,K46 ,
40 . K47 ,K48 ,K49 ,K410 ,K411 ,K412 ,K413 ,K414 ,K415 ,K416 ,
41 . K417 ,K418 ,K419 ,K420 ,K55 ,K56 ,K57 ,K58 ,K59 ,K510 ,
42 . K511 ,K512 ,K513 ,K514 ,K515 ,K516 ,K517 ,K518 ,K519 ,K520 ,
43 . K66 ,K67 ,K68 ,K69 ,K610 ,K611 ,K612 ,K613 ,K614 ,K615 ,
44 . K616 ,K617 ,K618 ,K619 ,K620 ,K77 ,K78 ,K79 ,K710 ,K711 ,
45 . K712 ,K713 ,K714 ,K715 ,K716 ,K717 ,K718 ,K719 ,K720 ,K88 ,
46 . K89 ,K810 ,K811 ,K812 ,K813 ,K814 ,K815 ,K816 ,K817 ,K818 ,
47 . K819 ,K820 ,K99 ,K910 ,K911 ,K912 ,K913 ,K914 ,K915 ,K916 ,
48 . K917 ,K918 ,K919 ,K920 ,K1010,K1011,K1012,K1013,K1014,K1015,
49 . K1016,K1017,K1018,K1019,K1020,K1111,K1112,K1113,K1114,K1115,
50 . K1116,K1117,K1118,K1119,K1120,K1212,K1213,K1214,K1215,K1216,
51 . K1217,K1218,K1219,K1220,K1313,K1314,K1315,K1316,K1317,K1318,
52 . K1319,K1320,K1414,K1415,K1416,K1417,K1418,K1419,K1420,K1515,
53 . K1516,K1517,K1518,K1519,K1520,K1616,K1617,K1618,K1619,K1620,
54 . K1717,K1718,K1719,K1720,K1818,K1819,K1820,K1919,K1920,K2020)
55C----6---------------------------------------------------------------7---------8
56C I m p l i c i t T y p e s
57C-----------------------------------------------
58#include "implicit_f.inc"
59C-----------------------------------------------
60C G l o b a l P a r a m e t e r s
61C-----------------------------------------------
62#include "mvsiz_p.inc"
63C-----------------------------------------------
64C C o m m o n B l o c k s
65C-----------------------------------------------
66C-----------------------------------------------------------------
67C D u m m y A r g u m e n t s
68C-----------------------------------------------
69 INTEGER IXS(NIXS,*),IXS20(12,*),NEL ,IDDL(*) ,NDOF(*) ,
70 . IADK(*) ,JDIK(*)
71C REAL
72 my_real
73 . K_DIAG(*) ,K_LT(*) ,OFF(*)
74 my_real
75 . K11(3,3,*) ,K12(3,3,*) ,K13(3,3,*) ,K14(3,3,*) ,K15(3,3,*) ,
76 . K16(3,3,*) ,K17(3,3,*) ,K18(3,3,*) ,K19(3,3,*) ,K110(3,3,*),
77 . K111(3,3,*),K112(3,3,*),K113(3,3,*),K114(3,3,*),K115(3,3,*),
78 . K116(3,3,*),K117(3,3,*),K118(3,3,*),K119(3,3,*),K120(3,3,*),
79 . K22(3,3,*) ,K23(3,3,*) ,K24(3,3,*) ,K25(3,3,*) ,K26(3,3,*) ,
80 . K27(3,3,*) ,K28(3,3,*) ,K29(3,3,*) ,K210(3,3,*),K211(3,3,*),
81 . K212(3,3,*),K213(3,3,*),K214(3,3,*),K215(3,3,*),K216(3,3,*),
82 . K217(3,3,*),K218(3,3,*),K219(3,3,*),K220(3,3,*),K33(3,3,*) ,
83 . K34(3,3,*) ,K35(3,3,*) ,K36(3,3,*) ,K37(3,3,*) ,K38(3,3,*) ,
84 . K39(3,3,*) ,K310(3,3,*),K311(3,3,*),K312(3,3,*),K313(3,3,*),
85 . K314(3,3,*),K315(3,3,*),K316(3,3,*),K317(3,3,*),K318(3,3,*),
86 . K319(3,3,*),K320(3,3,*),K44(3,3,*) ,K45(3,3,*) ,K46(3,3,*) ,
87 . K47(3,3,*) ,K48(3,3,*) ,K49(3,3,*) ,K410(3,3,*),K411(3,3,*),
88 . K412(3,3,*),K413(3,3,*),K414(3,3,*),K415(3,3,*),K416(3,3,*),
89 . K417(3,3,*),K418(3,3,*),K419(3,3,*),K420(3,3,*),K55(3,3,*) ,
90 . K56(3,3,*) ,K57(3,3,*) ,K58(3,3,*) ,K59(3,3,*) ,K510(3,3,*),
91 . K511(3,3,*),K512(3,3,*),K513(3,3,*),K514(3,3,*),K515(3,3,*),
92 . k516(3,3,*),k517(3,3,*),k518(3,3,*),k519(3,3,*),k520(3,3,*),
93 . k66(3,3,*) ,k67(3,3,*) ,k68(3,3,*) ,k69(3,3,*) ,k610(3,3,*),
94 . k611(3,3,*),k612(3,3,*),k613(3,3,*),k614(3,3,*),k615(3,3,*),
95 . k616(3,3,*),k617(3,3,*),k618(3,3,*),k619(3,3,*),k620(3,3,*),
96 . k77(3,3,*) ,k78(3,3,*) ,k79(3,3,*) ,k710(3,3,*),k711(3,3,*),
97 . k712(3,3,*),k713(3,3,*),k714(3,3,*),k715(3,3,*),k716(3,3,*),
98 . k717(3,3,*),k718(3,3,*),k719(3,3,*),k720(3,3,*),k88(3,3,*) ,
99 . k89(3,3,*) ,k810(3,3,*),k811(3,3,*),k812(3,3,*),k813(3,3,*),
100 . k814(3,3,*),k815(3,3,*),k816(3,3,*),k817(3,3,*),k818(3,3,*),
101 . k819(3,3,*),k820(3,3,*),k99(3,3,*) ,k910(3,3,*),k911(3,3,*),
102 . k912(3,3,*),k913(3,3,*),k914(3,3,*),k915(3,3,*),k916(3,3,*),
103 . k917(3,3,*),k918(3,3,*),k919(3,3,*),k920(3,3,*),k1010(3,3,*),
104 . k1011(3,3,*),k1012(3,3,*),k1013(3,3,*),k1014(3,3,*),k1015(3,3,*),
105 . k1016(3,3,*),k1017(3,3,*),k1018(3,3,*),k1019(3,3,*),k1020(3,3,*),
106 . k1111(3,3,*),k1112(3,3,*),k1113(3,3,*),k1114(3,3,*),k1115(3,3,*),
107 . k1116(3,3,*),k1117(3,3,*),k1118(3,3,*),k1119(3,3,*),k1120(3,3,*),
108 . k1212(3,3,*),k1213(3,3,*),k1214(3,3,*),k1215(3,3,*),k1216(3,3,*),
109 . k1217(3,3,*),k1218(3,3,*),k1219(3,3,*),k1220(3,3,*),k1313(3,3,*),
110 . k1314(3,3,*),k1315(3,3,*),k1316(3,3,*),k1317(3,3,*),k1318(3,3,*),
111 . k1319(3,3,*),k1320(3,3,*),k1414(3,3,*),k1415(3,3,*),k1416(3,3,*),
112 . k1417(3,3,*),k1418(3,3,*),k1419(3,3,*),k1420(3,3,*),k1515(3,3,*),
113 . k1516(3,3,*),k1517(3,3,*),k1518(3,3,*),k1519(3,3,*),k1520(3,3,*),
114 . k1616(3,3,*),k1617(3,3,*),k1618(3,3,*),k1619(3,3,*),k1620(3,3,*),
115 . k1717(3,3,*),k1718(3,3,*),k1719(3,3,*),k1720(3,3,*),k1818(3,3,*),
116 . k1819(3,3,*),k1820(3,3,*),k1919(3,3,*),k1920(3,3,*),k2020(3,3,*)
117C-----------------------------------------------
118C L o c a l V a r i a b l e s
119C-----------------------------------------------
120 INTEGER K,EP,NN(MVSIZ,20),J,I
121C----6---------------------------------------------------------------7---------8
122 DO K = 1,8
123 do ep = 1,nel
124 nn(ep,k)=ixs(k+1,ep)
125 ENDDO
126 ENDDO
127
128 DO k = 9,20
129 DO ep = 1,nel
130 nn(ep,k)=ixs20(k-8,ep)
131 ENDDO
132 ENDDO
133
134C-------block diagonal k11,k22...k88-----
135 CALL assem_kii(nn(1,1),nel,iddl,iadk,k_diag,k_lt ,k11,3,off)
136 CALL assem_kii(nn(1,2),nel,iddl,iadk,k_diag,k_lt ,k22,3,off)
137 CALL assem_kii(nn(1,3),nel,iddl,iadk,k_diag,k_lt ,k33,3,off)
138 CALL assem_kii(nn(1,4),nel,iddl,iadk,k_diag,k_lt ,k44,3,off)
139 CALL assem_kii(nn(1,5),nel,iddl,iadk,k_diag,k_lt ,k55,3,off)
140 CALL assem_kii(nn(1,6),nel,iddl,iadk,k_diag,k_lt ,k66,3,off)
141 CALL assem_kii(nn(1,7),nel,iddl,iadk,k_diag,k_lt ,k77,3,off)
142 CALL assem_kii(nn(1,8),nel,iddl,iadk,k_diag,k_lt ,k88,3,off)
143 CALL assem_kii(nn(1,9),nel,iddl,iadk,k_diag,k_lt ,k99,3,off)
144 CALL assem_kii(nn(1,10),nel,iddl,iadk,k_diag,k_lt ,k1010,3,off)
145 CALL assem_kii(nn(1,11),nel,iddl,iadk,k_diag,k_lt ,k1111,3,off)
146 CALL assem_kii(nn(1,12),nel,iddl,iadk,k_diag,k_lt ,k1212,3,off)
147 CALL assem_kii(nn(1,13),nel,iddl,iadk,k_diag,k_lt ,k1313,3,off)
148 CALL assem_kii(nn(1,14),nel,iddl,iadk,k_diag,k_lt ,k1414,3,off)
149 CALL assem_kii(nn(1,15),nel,iddl,iadk,k_diag,k_lt ,k1515,3,off)
150 CALL assem_kii(nn(1,16),nel,iddl,iadk,k_diag,k_lt ,k1616,3,off)
151 CALL assem_kii(nn(1,17),nel,iddl,iadk,k_diag,k_lt ,k1717,3,off)
152 CALL assem_kii(nn(1,18),nel,iddl,iadk,k_diag,k_lt ,k1818,3,off)
153 CALL assem_kii(nn(1,19),nel,iddl,iadk,k_diag,k_lt ,k1919,3,off)
154 CALL assem_kii(nn(1,20),nel,iddl,iadk,k_diag,k_lt ,k2020,3,off)
155C-------non diag kij-----
156 CALL assem_kij(nn(1,1),nn(1,2),nel,iddl,iadk,jdik,
157 . k_diag,k_lt,k12,3,off)
158 CALL assem_kij(nn(1,1),nn(1,3),nel,iddl,iadk,jdik,
159 . k_diag,k_lt,k13,3,off)
160 CALL assem_kij(nn(1,1),nn(1,4),nel,iddl,iadk,jdik,
161 . k_diag,k_lt,k14,3,off)
162 CALL assem_kij(nn(1,1),nn(1,5),nel,iddl,iadk,jdik,
163 . k_diag,k_lt,k15,3,off)
164 CALL assem_kij(nn(1,1),nn(1,6),nel,iddl,iadk,jdik,
165 . k_diag,k_lt,k16,3,off)
166 CALL assem_kij(nn(1,1),nn(1,7),nel,iddl,iadk,jdik,
167 . k_diag,k_lt,k17,3,off)
168 CALL assem_kij(nn(1,1),nn(1,8),nel,iddl,iadk,jdik,
169 . k_diag,k_lt,k18,3,off)
170 CALL assem_kij(nn(1,1),nn(1,9),nel,iddl,iadk,jdik,
171 . k_diag,k_lt,k19,3,off)
172 CALL assem_kij(nn(1,1),nn(1,10),nel,iddl,iadk,jdik,
173 . k_diag,k_lt,k110,3,off)
174 CALL assem_kij(nn(1,1),nn(1,11),nel,iddl,iadk,jdik,
175 . k_diag,k_lt,k111,3,off)
176 CALL assem_kij(nn(1,1),nn(1,12),nel,iddl,iadk,jdik,
177 . k_diag,k_lt,k112,3,off)
178 CALL assem_kij(nn(1,1),nn(1,13),nel,iddl,iadk,jdik,
179 . k_diag,k_lt,k113,3,off)
180 CALL assem_kij(nn(1,1),nn(1,14),nel,iddl,iadk,jdik,
181 . k_diag,k_lt,k114,3,off)
182 CALL assem_kij(nn(1,1),nn(1,15),nel,iddl,iadk,jdik,
183 . k_diag,k_lt,k115,3,off)
184 CALL assem_kij(nn(1,1),nn(1,16),nel,iddl,iadk,jdik,
185 . k_diag,k_lt,k116,3,off)
186 CALL assem_kij(nn(1,1),nn(1,17),nel,iddl,iadk,jdik,
187 . k_diag,k_lt,k117,3,off)
188 CALL assem_kij(nn(1,1),nn(1,18),nel,iddl,iadk,jdik,
189 . k_diag,k_lt,k118,3,off)
190 CALL assem_kij(nn(1,1),nn(1,19),nel,iddl,iadk,jdik,
191 . k_diag,k_lt,k119,3,off)
192 CALL assem_kij(nn(1,1),nn(1,20),nel,iddl,iadk,jdik,
193 . k_diag,k_lt,k120,3,off)
194
195
196
197 CALL assem_kij(nn(1,2),nn(1,3),nel,iddl,iadk,jdik,
198 . k_diag,k_lt,k23,3,off)
199 CALL assem_kij(nn(1,2),nn(1,4),nel,iddl,iadk,jdik,
200 . k_diag,k_lt,k24,3,off)
201 CALL assem_kij(nn(1,2),nn(1,5),nel,iddl,iadk,jdik,
202 . k_diag,k_lt,k25,3,off)
203 CALL assem_kij(nn(1,2),nn(1,6),nel,iddl,iadk,jdik,
204 . k_diag,k_lt,k26,3,off)
205 CALL assem_kij(nn(1,2),nn(1,7),nel,iddl,iadk,jdik,
206 . k_diag,k_lt,k27,3,off)
207 CALL assem_kij(nn(1,2),nn(1,8),nel,iddl,iadk,jdik,
208 . k_diag,k_lt,k28,3,off)
209 CALL assem_kij(nn(1,2),nn(1,9),nel,iddl,iadk,jdik,
210 . k_diag,k_lt,k29,3,off)
211 CALL assem_kij(nn(1,2),nn(1,10),nel,iddl,iadk,jdik,
212 . k_diag,k_lt,k210,3,off)
213 CALL assem_kij(nn(1,2),nn(1,11),nel,iddl,iadk,jdik,
214 . k_diag,k_lt,k211,3,off)
215 CALL assem_kij(nn(1,2),nn(1,12),nel,iddl,iadk,jdik,
216 . k_diag,k_lt,k212,3,off)
217 CALL assem_kij(nn(1,2),nn(1,13),nel,iddl,iadk,jdik,
218 . k_diag,k_lt,k213,3,off)
219 CALL assem_kij(nn(1,2),nn(1,14),nel,iddl,iadk,jdik,
220 . k_diag,k_lt,k214,3,off)
221 CALL assem_kij(nn(1,2),nn(1,15),nel,iddl,iadk,jdik,
222 . k_diag,k_lt,k215,3,off)
223 CALL assem_kij(nn(1,2),nn(1,16),nel,iddl,iadk,jdik,
224 . k_diag,k_lt,k216,3,off)
225 CALL assem_kij(nn(1,2),nn(1,17),nel,iddl,iadk,jdik,
226 . k_diag,k_lt,k217,3,off)
227 CALL assem_kij(nn(1,2),nn(1,18),nel,iddl,iadk,jdik,
228 . k_diag,k_lt,k218,3,off)
229 CALL assem_kij(nn(1,2),nn(1,19),nel,iddl,iadk,jdik,
230 . k_diag,k_lt,k219,3,off)
231 CALL assem_kij(nn(1,2),nn(1,20),nel,iddl,iadk,jdik,
232 . k_diag,k_lt,k220,3,off)
233
234
235 CALL assem_kij(nn(1,3),nn(1,4),nel,iddl,iadk,jdik,
236 . k_diag,k_lt,k34,3,off)
237 CALL assem_kij(nn(1,3),nn(1,5),nel,iddl,iadk,jdik,
238 . k_diag,k_lt,k35,3,off)
239 CALL assem_kij(nn(1,3),nn(1,6),nel,iddl,iadk,jdik,
240 . k_diag,k_lt,k36,3,off)
241 CALL assem_kij(nn(1,3),nn(1,7),nel,iddl,iadk,jdik,
242 . k_diag,k_lt,k37,3,off)
243 CALL assem_kij(nn(1,3),nn(1,8),nel,iddl,iadk,jdik,
244 . k_diag,k_lt,k38,3,off)
245 CALL assem_kij(nn(1,3),nn(1,9),nel,iddl,iadk,jdik,
246 . k_diag,k_lt,k39,3,off)
247 CALL assem_kij(nn(1,3),nn(1,10),nel,iddl,iadk,jdik,
248 . k_diag,k_lt,k310,3,off)
249 CALL assem_kij(nn(1,3),nn(1,11),nel,iddl,iadk,jdik,
250 . k_diag,k_lt,k311,3,off)
251 CALL assem_kij(nn(1,3),nn(1,12),nel,iddl,iadk,jdik,
252 . k_diag,k_lt,k312,3,off)
253 CALL assem_kij(nn(1,3),nn(1,13),nel,iddl,iadk,jdik,
254 . k_diag,k_lt,k313,3,off)
255 CALL assem_kij(nn(1,3),nn(1,14),nel,iddl,iadk,jdik,
256 . k_diag,k_lt,k314,3,off)
257 CALL assem_kij(nn(1,3),nn(1,15),nel,iddl,iadk,jdik,
258 . k_diag,k_lt,k315,3,off)
259 CALL assem_kij(nn(1,3),nn(1,16),nel,iddl,iadk,jdik,
260 . k_diag,k_lt,k316,3,off)
261 CALL assem_kij(nn(1,3),nn(1,17),nel,iddl,iadk,jdik,
262 . k_diag,k_lt,k317,3,off)
263 CALL assem_kij(nn(1,3),nn(1,18),nel,iddl,iadk,jdik,
264 . k_diag,k_lt,k318,3,off)
265 CALL assem_kij(nn(1,3),nn(1,19),nel,iddl,iadk,jdik,
266 . k_diag,k_lt,k319,3,off)
267 CALL assem_kij(nn(1,3),nn(1,20),nel,iddl,iadk,jdik,
268 . k_diag,k_lt,k320,3,off)
269
270
271
272 CALL assem_kij(nn(1,4),nn(1,5),nel,iddl,iadk,jdik,
273 . k_diag,k_lt,k45,3,off)
274 CALL assem_kij(nn(1,4),nn(1,6),nel,iddl,iadk,jdik,
275 . k_diag,k_lt,k46,3,off)
276 CALL assem_kij(nn(1,4),nn(1,7),nel,iddl,iadk,jdik,
277 . k_diag,k_lt,k47,3,off)
278 CALL assem_kij(nn(1,4),nn(1,8),nel,iddl,iadk,jdik,
279 . k_diag,k_lt,k48,3,off)
280 CALL assem_kij(nn(1,4),nn(1,9),nel,iddl,iadk,jdik,
281 . k_diag,k_lt,k49,3,off)
282 CALL assem_kij(nn(1,4),nn(1,10),nel,iddl,iadk,jdik,
283 . k_diag,k_lt,k410,3,off)
284 CALL assem_kij(nn(1,4),nn(1,11),nel,iddl,iadk,jdik,
285 . k_diag,k_lt,k411,3,off)
286 CALL assem_kij(nn(1,4),nn(1,12),nel,iddl,iadk,jdik,
287 . k_diag,k_lt,k412,3,off)
288 CALL assem_kij(nn(1,4),nn(1,13),nel,iddl,iadk,jdik,
289 . k_diag,k_lt,k413,3,off)
290 CALL assem_kij(nn(1,4),nn(1,14),nel,iddl,iadk,jdik,
291 . k_diag,k_lt,k414,3,off)
292 CALL assem_kij(nn(1,4),nn(1,15),nel,iddl,iadk,jdik,
293 . k_diag,k_lt,k415,3,off)
294 CALL assem_kij(nn(1,4),nn(1,16),nel,iddl,iadk,jdik,
295 . k_diag,k_lt,k416,3,off)
296 CALL assem_kij(nn(1,4),nn(1,17),nel,iddl,iadk,jdik,
297 . k_diag,k_lt,k417,3,off)
298 CALL assem_kij(nn(1,4),nn(1,18),nel,iddl,iadk,jdik,
299 . k_diag,k_lt,k418,3,off)
300 CALL assem_kij(nn(1,4),nn(1,19),nel,iddl,iadk,jdik,
301 . k_diag,k_lt,k419,3,off)
302 CALL assem_kij(nn(1,4),nn(1,20),nel,iddl,iadk,jdik,
303 . k_diag,k_lt,k420,3,off)
304
305
306
307 CALL assem_kij(nn(1,5),nn(1,6),nel,iddl,iadk,jdik,
308 . k_diag,k_lt,k56,3,off)
309 CALL assem_kij(nn(1,5),nn(1,7),nel,iddl,iadk,jdik,
310 . k_diag,k_lt,k57,3,off)
311 CALL assem_kij(nn(1,5),nn(1,8),nel,iddl,iadk,jdik,
312 . k_diag,k_lt,k58,3,off)
313 CALL assem_kij(nn(1,5),nn(1,9),nel,iddl,iadk,jdik,
314 . k_diag,k_lt,k59,3,off)
315 CALL assem_kij(nn(1,5),nn(1,10),nel,iddl,iadk,jdik,
316 . k_diag,k_lt,k510,3,off)
317 CALL assem_kij(nn(1,5),nn(1,11),nel,iddl,iadk,jdik,
318 . k_diag,k_lt,k511,3,off)
319 CALL assem_kij(nn(1,5),nn(1,12),nel,iddl,iadk,jdik,
320 . k_diag,k_lt,k512,3,off)
321 CALL assem_kij(nn(1,5),nn(1,13),nel,iddl,iadk,jdik,
322 . k_diag,k_lt,k513,3,off)
323 CALL assem_kij(nn(1,5),nn(1,14),nel,iddl,iadk,jdik,
324 . k_diag,k_lt,k514,3,off)
325 CALL assem_kij(nn(1,5),nn(1,15),nel,iddl,iadk,jdik,
326 . k_diag,k_lt,k515,3,off)
327 CALL assem_kij(nn(1,5),nn(1,16),nel,iddl,iadk,jdik,
328 . k_diag,k_lt,k516,3,off)
329 CALL assem_kij(nn(1,5),nn(1,17),nel,iddl,iadk,jdik,
330 . k_diag,k_lt,k517,3,off)
331 CALL assem_kij(nn(1,5),nn(1,18),nel,iddl,iadk,jdik,
332 . k_diag,k_lt,k518,3,off)
333 CALL assem_kij(nn(1,5),nn(1,19),nel,iddl,iadk,jdik,
334 . k_diag,k_lt,k519,3,off)
335 CALL assem_kij(nn(1,5),nn(1,20),nel,iddl,iadk,jdik,
336 . k_diag,k_lt,k520,3,off)
337
338
339
340 CALL assem_kij(nn(1,6),nn(1,7),nel,iddl,iadk,jdik,
341 . k_diag,k_lt,k67,3,off)
342 CALL assem_kij(nn(1,6),nn(1,8),nel,iddl,iadk,jdik,
343 . k_diag,k_lt,k68,3,off)
344 CALL assem_kij(nn(1,6),nn(1,9),nel,iddl,iadk,jdik,
345 . k_diag,k_lt,k69,3,off)
346 CALL assem_kij(nn(1,6),nn(1,10),nel,iddl,iadk,jdik,
347 . k_diag,k_lt,k610,3,off)
348 CALL assem_kij(nn(1,6),nn(1,11),nel,iddl,iadk,jdik,
349 . k_diag,k_lt,k611,3,off)
350 CALL assem_kij(nn(1,6),nn(1,12),nel,iddl,iadk,jdik,
351 . k_diag,k_lt,k612,3,off)
352 CALL assem_kij(nn(1,6),nn(1,13),nel,iddl,iadk,jdik,
353 . k_diag,k_lt,k613,3,off)
354 CALL assem_kij(nn(1,6),nn(1,14),nel,iddl,iadk,jdik,
355 . k_diag,k_lt,k614,3,off)
356 CALL assem_kij(nn(1,6),nn(1,15),nel,iddl,iadk,jdik,
357 . k_diag,k_lt,k615,3,off)
358 CALL assem_kij(nn(1,6),nn(1,16),nel,iddl,iadk,jdik,
359 . k_diag,k_lt,k616,3,off)
360 CALL assem_kij(nn(1,6),nn(1,17),nel,iddl,iadk,jdik,
361 . k_diag,k_lt,k617,3,off)
362 CALL assem_kij(nn(1,6),nn(1,18),nel,iddl,iadk,jdik,
363 . k_diag,k_lt,k618,3,off)
364 CALL assem_kij(nn(1,6),nn(1,19),nel,iddl,iadk,jdik,
365 . k_diag,k_lt,k619,3,off)
366 CALL assem_kij(nn(1,6),nn(1,20),nel,iddl,iadk,jdik,
367 . k_diag,k_lt,k620,3,off)
368
369
370
371
372 CALL assem_kij(nn(1,7),nn(1,8),nel,iddl,iadk,jdik,
373 . k_diag, k_lt,k78,3,off)
374 CALL assem_kij(nn(1,7),nn(1,9),nel,iddl,iadk,jdik,
375 . k_diag,k_lt,k79,3,off)
376 CALL assem_kij(nn(1,7),nn(1,10),nel,iddl,iadk,jdik,
377 . k_diag,k_lt,k710,3,off)
378 CALL assem_kij(nn(1,7),nn(1,11),nel,iddl,iadk,jdik,
379 . k_diag,k_lt,k711,3,off)
380 CALL assem_kij(nn(1,7),nn(1,12),nel,iddl,iadk,jdik,
381 . k_diag,k_lt,k712,3,off)
382 CALL assem_kij(nn(1,7),nn(1,13),nel,iddl,iadk,jdik,
383 . k_diag,k_lt,k713,3,off)
384 CALL assem_kij(nn(1,7),nn(1,14),nel,iddl,iadk,jdik,
385 . k_diag,k_lt,k714,3,off)
386 CALL assem_kij(nn(1,7),nn(1,15),nel,iddl,iadk,jdik,
387 . k_diag,k_lt,k715,3,off)
388 CALL assem_kij(nn(1,7),nn(1,16),nel,iddl,iadk,jdik,
389 . k_diag,k_lt,k716,3,off)
390 CALL assem_kij(nn(1,7),nn(1,17),nel,iddl,iadk,jdik,
391 . k_diag,k_lt,k717,3,off)
392 CALL assem_kij(nn(1,7),nn(1,18),nel,iddl,iadk,jdik,
393 . k_diag,k_lt,k718,3,off)
394 CALL assem_kij(nn(1,7),nn(1,19),nel,iddl,iadk,jdik,
395 . k_diag,k_lt,k719,3,off)
396 CALL assem_kij(nn(1,7),nn(1,20),nel,iddl,iadk,jdik,
397 . k_diag,k_lt,k720,3,off)
398
399
400 CALL assem_kij(nn(1,8),nn(1,9),nel,iddl,iadk,jdik,
401 . k_diag,k_lt,k89,3,off)
402 CALL assem_kij(nn(1,8),nn(1,10),nel,iddl,iadk,jdik,
403 . k_diag,k_lt,k810,3,off)
404 CALL assem_kij(nn(1,8),nn(1,11),nel,iddl,iadk,jdik,
405 . k_diag,k_lt,k811,3,off)
406 CALL assem_kij(nn(1,8),nn(1,12),nel,iddl,iadk,jdik,
407 . k_diag,k_lt,k812,3,off)
408 CALL assem_kij(nn(1,8),nn(1,13),nel,iddl,iadk,jdik,
409 . k_diag,k_lt,k813,3,off)
410 CALL assem_kij(nn(1,8),nn(1,14),nel,iddl,iadk,jdik,
411 . k_diag,k_lt,k814,3,off)
412 CALL assem_kij(nn(1,8),nn(1,15),nel,iddl,iadk,jdik,
413 . k_diag,k_lt,k815,3,off)
414 CALL assem_kij(nn(1,8),nn(1,16),nel,iddl,iadk,jdik,
415 . k_diag,k_lt,k816,3,off)
416 CALL assem_kij(nn(1,8),nn(1,17),nel,iddl,iadk,jdik,
417 . k_diag,k_lt,k817,3,off)
418 CALL assem_kij(nn(1,8),nn(1,18),nel,iddl,iadk,jdik,
419 . k_diag,k_lt,k818,3,off)
420 CALL assem_kij(nn(1,8),nn(1,19),nel,iddl,iadk,jdik,
421 . k_diag,k_lt,k819,3,off)
422 CALL assem_kij(nn(1,8),nn(1,20),nel,iddl,iadk,jdik,
423 . k_diag,k_lt,k820,3,off)
424
425
426 CALL assem_kij(nn(1,9),nn(1,10),nel,iddl,iadk,jdik,
427 . k_diag,k_lt,k910,3,off)
428 CALL assem_kij(nn(1,9),nn(1,11),nel,iddl,iadk,jdik,
429 . k_diag,k_lt,k911,3,off)
430 CALL assem_kij(nn(1,9),nn(1,12),nel,iddl,iadk,jdik,
431 . k_diag,k_lt,k912,3,off)
432 CALL assem_kij(nn(1,9),nn(1,13),nel,iddl,iadk,jdik,
433 . k_diag,k_lt,k913,3,off)
434 CALL assem_kij(nn(1,9),nn(1,14),nel,iddl,iadk,jdik,
435 . k_diag,k_lt,k914,3,off)
436 CALL assem_kij(nn(1,9),nn(1,15),nel,iddl,iadk,jdik,
437 . k_diag,k_lt,k915,3,off)
438 CALL assem_kij(nn(1,9),nn(1,16),nel,iddl,iadk,jdik,
439 . k_diag,k_lt,k916,3,off)
440 CALL assem_kij(nn(1,9),nn(1,17),nel,iddl,iadk,jdik,
441 . k_diag,k_lt,k917,3,off)
442 CALL assem_kij(nn(1,9),nn(1,18),nel,iddl,iadk,jdik,
443 . k_diag,k_lt,k918,3,off)
444 CALL assem_kij(nn(1,9),nn(1,19),nel,iddl,iadk,jdik,
445 . k_diag,k_lt,k919,3,off)
446 CALL assem_kij(nn(1,9),nn(1,20),nel,iddl,iadk,jdik,
447 . k_diag,k_lt,k920,3,off)
448
449
450 CALL assem_kij(nn(1,10),nn(1,11),nel,iddl,iadk,jdik,
451 . k_diag,k_lt,k1011,3,off)
452 CALL assem_kij(nn(1,10),nn(1,12),nel,iddl,iadk,jdik,
453 . k_diag,k_lt,k1012,3,off)
454 CALL assem_kij(nn(1,10),nn(1,13),nel,iddl,iadk,jdik,
455 . k_diag,k_lt,k1013,3,off)
456 CALL assem_kij(nn(1,10),nn(1,14),nel,iddl,iadk,jdik,
457 . k_diag,k_lt,k1014,3,off)
458 CALL assem_kij(nn(1,10),nn(1,15),nel,iddl,iadk,jdik,
459 . k_diag,k_lt,k1015,3,off)
460 CALL assem_kij(nn(1,10),nn(1,16),nel,iddl,iadk,jdik,
461 . k_diag,k_lt,k1016,3,off)
462 CALL assem_kij(nn(1,10),nn(1,17),nel,iddl,iadk,jdik,
463 . k_diag,k_lt,k1017,3,off)
464 CALL assem_kij(nn(1,10),nn(1,18),nel,iddl,iadk,jdik,
465 . k_diag,k_lt,k1018,3,off)
466 CALL assem_kij(nn(1,10),nn(1,19),nel,iddl,iadk,jdik,
467 . k_diag,k_lt,k1019,3,off)
468 CALL assem_kij(nn(1,10),nn(1,20),nel,iddl,iadk,jdik,
469 . k_diag,k_lt,k1020,3,off)
470
471
472 CALL assem_kij(nn(1,11),nn(1,12),nel,iddl,iadk,jdik,
473 . k_diag,k_lt,k1112,3,off)
474 CALL assem_kij(nn(1,11),nn(1,13),nel,iddl,iadk,jdik,
475 . k_diag,k_lt,k1113,3,off)
476 CALL assem_kij(nn(1,11),nn(1,14),nel,iddl,iadk,jdik,
477 . k_diag,k_lt,k1114,3,off)
478 CALL assem_kij(nn(1,11),nn(1,15),nel,iddl,iadk,jdik,
479 . k_diag,k_lt,k1115,3,off)
480 CALL assem_kij(nn(1,11),nn(1,16),nel,iddl,iadk,jdik,
481 . k_diag,k_lt,k1116,3,off)
482 CALL assem_kij(nn(1,11),nn(1,17),nel,iddl,iadk,jdik,
483 . k_diag,k_lt,k1117,3,off)
484 CALL assem_kij(nn(1,11),nn(1,18),nel,iddl,iadk,jdik,
485 . k_diag,k_lt,k1118,3,off)
486 CALL assem_kij(nn(1,11),nn(1,19),nel,iddl,iadk,jdik,
487 . k_diag,k_lt,k1119,3,off)
488 CALL assem_kij(nn(1,11),nn(1,20),nel,iddl,iadk,jdik,
489 . k_diag,k_lt,k1120,3,off)
490
491
492
493 CALL assem_kij(nn(1,12),nn(1,13),nel,iddl,iadk,jdik,
494 . k_diag,k_lt,k1213,3,off)
495 CALL assem_kij(nn(1,12),nn(1,14),nel,iddl,iadk,jdik,
496 . k_diag,k_lt,k1214,3,off)
497 CALL assem_kij(nn(1,12),nn(1,15),nel,iddl,iadk,jdik,
498 . k_diag,k_lt,k1215,3,off)
499 CALL assem_kij(nn(1,12),nn(1,16),nel,iddl,iadk,jdik,
500 . k_diag,k_lt,k1216,3,off)
501 CALL assem_kij(nn(1,12),nn(1,17),nel,iddl,iadk,jdik,
502 . k_diag,k_lt,k1217,3,off)
503 CALL assem_kij(nn(1,12),nn(1,18),nel,iddl,iadk,jdik,
504 . k_diag,k_lt,k1218,3,off)
505 CALL assem_kij(nn(1,12),nn(1,19),nel,iddl,iadk,jdik,
506 . k_diag,k_lt,k1219,3,off)
507 CALL assem_kij(nn(1,12),nn(1,20),nel,iddl,iadk,jdik,
508 . k_diag,k_lt,k1220,3,off)
509
510
511 CALL assem_kij(nn(1,13),nn(1,14),nel,iddl,iadk,jdik,
512 . k_diag,k_lt,k1314,3,off)
513 CALL assem_kij(nn(1,13),nn(1,15),nel,iddl,iadk,jdik,
514 . k_diag,k_lt,k1315,3,off)
515 CALL assem_kij(nn(1,13),nn(1,16),nel,iddl,iadk,jdik,
516 . k_diag,k_lt,k1316,3,off)
517 CALL assem_kij(nn(1,13),nn(1,17),nel,iddl,iadk,jdik,
518 . k_diag,k_lt,k1317,3,off)
519 CALL assem_kij(nn(1,13),nn(1,18),nel,iddl,iadk,jdik,
520 . k_diag,k_lt,k1318,3,off)
521 CALL assem_kij(nn(1,13),nn(1,19),nel,iddl,iadk,jdik,
522 . k_diag,k_lt,k1319,3,off)
523 CALL assem_kij(nn(1,13),nn(1,20),nel,iddl,iadk,jdik,
524 . k_diag,k_lt,k1320,3,off)
525
526
527 CALL assem_kij(nn(1,14),nn(1,15),nel,iddl,iadk,jdik,
528 . k_diag,k_lt,k1415,3,off)
529 CALL assem_kij(nn(1,14),nn(1,16),nel,iddl,iadk,jdik,
530 . k_diag,k_lt,k1416,3,off)
531 CALL assem_kij(nn(1,14),nn(1,17),nel,iddl,iadk,jdik,
532 . k_diag,k_lt,k1417,3,off)
533 CALL assem_kij(nn(1,14),nn(1,18),nel,iddl,iadk,jdik,
534 . k_diag,k_lt,k1418,3,off)
535 CALL assem_kij(nn(1,14),nn(1,19),nel,iddl,iadk,jdik,
536 . k_diag,k_lt,k1419,3,off)
537 CALL assem_kij(nn(1,14),nn(1,20),nel,iddl,iadk,jdik,
538 . k_diag,k_lt,k1420,3,off)
539
540
541 CALL assem_kij(nn(1,15),nn(1,16),nel,iddl,iadk,jdik,
542 . k_diag,k_lt,k1516,3,off)
543 CALL assem_kij(nn(1,15),nn(1,17),nel,iddl,iadk,jdik,
544 . k_diag,k_lt,k1517,3,off)
545 CALL assem_kij(nn(1,15),nn(1,18),nel,iddl,iadk,jdik,
546 . k_diag,k_lt,k1518,3,off)
547 CALL assem_kij(nn(1,15),nn(1,19),nel,iddl,iadk,jdik,
548 . k_diag,k_lt,k1519,3,off)
549 CALL assem_kij(nn(1,15),nn(1,20),nel,iddl,iadk,jdik,
550 . k_diag,k_lt,k1520,3,off)
551
552
553 CALL assem_kij(nn(1,16),nn(1,17),nel,iddl,iadk,jdik,
554 . k_diag,k_lt,k1617,3,off)
555 CALL assem_kij(nn(1,16),nn(1,18),nel,iddl,iadk,jdik,
556 . k_diag,k_lt,k1618,3,off)
557 CALL assem_kij(nn(1,16),nn(1,19),nel,iddl,iadk,jdik,
558 . k_diag,k_lt,k1619,3,off)
559 CALL assem_kij(nn(1,16),nn(1,20),nel,iddl,iadk,jdik,
560 . k_diag,k_lt,k1620,3,off)
561
562
563 CALL assem_kij(nn(1,17),nn(1,18),nel,iddl,iadk,jdik,
564 . k_diag,k_lt,k1718,3,off)
565 CALL assem_kij(nn(1,17),nn(1,19),nel,iddl,iadk,jdik,
566 . k_diag,k_lt,k1719,3,off)
567 CALL assem_kij(nn(1,17),nn(1,20),nel,iddl,iadk,jdik,
568 . k_diag,k_lt,k1720,3,off)
569
570
571 CALL assem_kij(nn(1,18),nn(1,19),nel,iddl,iadk,jdik,
572 . k_diag,k_lt,k1819,3,off)
573 CALL assem_kij(nn(1,18),nn(1,20),nel,iddl,iadk,jdik,
574 . k_diag,k_lt,k1820,3,off)
575
576
577 CALL assem_kij(nn(1,19),nn(1,20),nel,iddl,iadk,jdik,
578 . k_diag,k_lt,k1920,3,off)
579
580C
581C----6---------------------------------------------------------------7---------8
582 RETURN
583 END
assem_s20
subroutine assem_s20(ixs, ixs20, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, off, k11, k12, k13, k14, k15, k16, k17, k18, k19, k110, k111, k112, k113, k114, k115, k116, k117, k118, k119, k120, k22, k23, k24, k25, k26, k27, k28, k29, k210, k211, k212, k213, k214, k215, k216, k217, k218, k219, k220, k33, k34, k35, k36, k37, k38, k39, k310, k311, k312, k313, k314, k315, k316, k317, k318, k319, k320, k44, k45, k46, k47, k48, k49, k410, k411, k412, k413, k414, k415, k416, k417, k418, k419, k420, k55, k56, k57, k58, k59, k510, k511, k512, k513, k514, k515, k516, k517, k518, k519, k520, k66, k67, k68, k69, k610, k611, k612, k613, k614, k615, k616, k617, k618, k619, k620, k77, k78, k79, k710, k711, k712, k713, k714, k715, k716, k717, k718, k719, k720, k88, k89, k810, k811, k812, k813, k814, k815, k816, k817, k818, k819, k820, k99, k910, k911, k912, k913, k914, k915, k916, k917, k918, k919, k920, k1010, k1011, k1012, k1013, k1014, k1015, k1016, k1017, k1018, k1019, k1020, k1111, k1112, k1113, k1114, k1115, k1116, k1117, k1118, k1119, k1120, k1212, k1213, k1214, k1215, k1216, k1217, k1218, k1219, k1220, k1313, k1314, k1315, k1316, k1317, k1318, k1319, k1320, k1414, k1415, k1416, k1417, k1418, k1419, k1420, k1515, k1516, k1517, k1518, k1519, k1520, k1616, k1617, k1618, k1619, k1620, k1717, k1718, k1719, k1720, k1818, k1819, k1820, k1919, k1920, k2020)
Definition assem_s20.F:55
assem_kii
subroutine assem_kii(ni, nel, iddl, iadk, k_diag, k_lt, kii, nd, off)
Definition imp_glob_k.F:962
assem_kij
subroutine assem_kij(ni, nj, nel, iddl, iadk, jdik, k_diag, k_lt, kij, nd, off)
Definition imp_glob_k.F:1041