OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s20ke3.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!|| s20ke3 ../engine/source/elements/solid/solide20/s20ke3.F
25!||--- called by ------------------------------------------------------
26!|| imp_glob_k ../engine/source/implicit/imp_glob_k.F
27!|| imp_glob_k0 ../engine/source/implicit/imp_glob_k.F
28!||--- calls -----------------------------------------------------
29!|| assem_s20 ../engine/source/implicit/assem_s20.F
30!|| mmats ../engine/source/elements/solid/solide8z/mmats.F
31!|| mmstifs ../engine/source/elements/solid/solide8z/mmstifs.F
32!|| s20coork ../engine/source/elements/solid/solide20/s20coork.F
33!|| s20cumg3 ../engine/source/elements/solid/solide20/s20cumg3.F
34!|| s20deri3 ../engine/source/elements/solid/solide20/s20deri3.F
35!|| s20eoff ../engine/source/elements/solid/solide20/s20eoff.F
36!|| s20kgeo3 ../engine/source/elements/solid/solide20/s20kgeo3.f
37!|| s20rst ../engine/source/elements/solid/solide20/s20rst.F
38!|| sorthdir3 ../engine/source/elements/solid/solide/sorthdir3.F
39!|| sreploc3 ../engine/source/elements/solid/solide/sreploc3.F
40!||--- uses -----------------------------------------------------
41!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
42!||====================================================================
43 SUBROUTINE s20ke3(
44 1 PM, GEO, IXS, IXS20,
45 2 X, ELBUF_STR,ETAG, IDDL,
46 3 NDOF, K_DIAG, K_LT, IADK,
47 4 JDIK, NEL, IPM, IGEO,
48 5 IKGEO, BUFMAT, NFT, MTN,
49 6 ISMSTR, JHBE, IREP, IGTYP,
50 7 ISORTH)
51C-----------------------------------------------
52C M o d u l e s
53C-----------------------------------------------
54 USE elbufdef_mod
55C-----------------------------------------------
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-----------------------------------------------
66#include "com04_c.inc"
67#include "param_c.inc"
68C-----------------------------------------------
69C D u m m y A r g u m e n t s
70C-----------------------------------------------
71 INTEGER, INTENT(IN) :: NFT
72 INTEGER, INTENT(IN) :: MTN
73 INTEGER, INTENT(IN) :: ISMSTR
74 INTEGER, INTENT(IN) :: JHBE
75 INTEGER, INTENT(IN) :: IREP
76 INTEGER, INTENT(IN) :: IGTYP
77 INTEGER, INTENT(IN) :: ISORTH
78 INTEGER IXS(NIXS,*),IXS20(12,*), IKGEO
79 INTEGER NEL ,IPM(NPROPMI,*),IGEO(NPROPGI,*)
80C
81 my_real
82 . pm(npropm,*), geo(npropg,*), x(*),
83 . k11(9,mvsiz) ,k12(9,mvsiz) ,k13(9,mvsiz) ,
84 . k14(9,mvsiz) ,k15(9,mvsiz) ,
85 . k16(9,mvsiz) ,k17(9,mvsiz) ,k18(9,mvsiz) ,
86 . k19(9,mvsiz) ,k110(9,mvsiz) ,
87 . k111(9,mvsiz) ,k112(9,mvsiz) ,k113(9,mvsiz) ,
88 . k114(9,mvsiz) ,k115(9,mvsiz) ,
89 . k116(9,mvsiz) ,k117(9,mvsiz) ,k118(9,mvsiz) ,
90 . k119(9,mvsiz) ,k120(9,mvsiz) ,
91 . k22(9,mvsiz) ,k23(9,mvsiz) ,k24(9,mvsiz) ,
92 . k25(9,mvsiz) ,k26(9,mvsiz) ,
93 . k27(9,mvsiz) ,k28(9,mvsiz) ,k29(9,mvsiz) ,
94 . k210(9,mvsiz) ,k211(9,mvsiz) ,
95 . k212(9,mvsiz) ,k213(9,mvsiz) ,k214(9,mvsiz) ,
96 . k215(9,mvsiz) ,k216(9,mvsiz) ,
97 . k217(9,mvsiz) ,k218(9,mvsiz) ,k219(9,mvsiz) ,
98 . k220(9,mvsiz) ,k33(9,mvsiz) ,
99 . k34(9,mvsiz) ,k35(9,mvsiz) ,k36(9,mvsiz) ,
100 . k37(9,mvsiz) ,k38(9,mvsiz) ,
101 . k39(9,mvsiz) ,k310(9,mvsiz) ,k311(9,mvsiz) ,
102 . k312(9,mvsiz) ,k313(9,mvsiz) ,
103 . k314(9,mvsiz) ,k315(9,mvsiz) ,k316(9,mvsiz) ,
104 . k317(9,mvsiz) ,k318(9,mvsiz) ,
105 . k319(9,mvsiz) ,k320(9,mvsiz) ,k44(9,mvsiz) ,
106 . k45(9,mvsiz) ,k46(9,mvsiz) ,
107 . k47(9,mvsiz) ,k48(9,mvsiz) ,k49(9,mvsiz) ,
108 . k410(9,mvsiz) ,k411(9,mvsiz) ,
109 . k412(9,mvsiz) ,k413(9,mvsiz) ,k414(9,mvsiz) ,
110 . k415(9,mvsiz) ,k416(9,mvsiz) ,
111 . k417(9,mvsiz) ,k418(9,mvsiz) ,k419(9,mvsiz) ,
112 . k420(9,mvsiz) ,k55(9,mvsiz) ,
113 . k56(9,mvsiz) ,k57(9,mvsiz) ,k58(9,mvsiz) ,
114 . k59(9,mvsiz) ,k510(9,mvsiz) ,
115 . k511(9,mvsiz) ,k512(9,mvsiz) ,k513(9,mvsiz) ,
116 . k514(9,mvsiz) ,k515(9,mvsiz) ,
117 . k516(9,mvsiz) ,k517(9,mvsiz) ,k518(9,mvsiz) ,
118 . k519(9,mvsiz) ,k520(9,mvsiz) ,
119 . k66(9,mvsiz) ,k67(9,mvsiz) ,k68(9,mvsiz) ,
120 . k69(9,mvsiz) ,k610(9,mvsiz) ,
121 . k611(9,mvsiz) ,k612(9,mvsiz) ,k613(9,mvsiz) ,
122 . k614(9,mvsiz) ,k615(9,mvsiz) ,
123 . k616(9,mvsiz) ,k617(9,mvsiz) ,k618(9,mvsiz) ,
124 . k619(9,mvsiz) ,k620(9,mvsiz) ,
125 . k77(9,mvsiz) ,k78(9,mvsiz) ,k79(9,mvsiz) ,
126 . k710(9,mvsiz) ,k711(9,mvsiz) ,
127 . k712(9,mvsiz) ,k713(9,mvsiz) ,k714(9,mvsiz) ,
128 . k715(9,mvsiz) ,k716(9,mvsiz) ,
129 . k717(9,mvsiz) ,k718(9,mvsiz) ,k719(9,mvsiz) ,
130 . k720(9,mvsiz) ,k88(9,mvsiz) ,
131 . k89(9,mvsiz) ,k810(9,mvsiz) ,k811(9,mvsiz) ,
132 . k812(9,mvsiz) ,k813(9,mvsiz) ,
133 . k814(9,mvsiz) ,k815(9,mvsiz) ,k816(9,mvsiz) ,
134 . k817(9,mvsiz) ,k818(9,mvsiz) ,
135 . k819(9,mvsiz) ,k820(9,mvsiz) ,k99(9,mvsiz) ,
136 . k910(9,mvsiz) ,k911(9,mvsiz) ,
137 . k912(9,mvsiz) ,k913(9,mvsiz) ,k914(9,mvsiz) ,
138 . k915(9,mvsiz) ,k916(9,mvsiz) ,
139 . k917(9,mvsiz) ,k918(9,mvsiz) ,k919(9,mvsiz) ,
140 . k920(9,mvsiz) ,k1010(9,mvsiz) ,
141 . k1011(9,mvsiz) ,k1012(9,mvsiz) ,k1013(9,mvsiz) ,
142 . k1014(9,mvsiz) ,k1015(9,mvsiz) ,
143 . k1016(9,mvsiz) ,k1017(9,mvsiz) ,k1018(9,mvsiz) ,
144 . k1019(9,mvsiz) ,k1020(9,mvsiz) ,
145 . k1111(9,mvsiz) ,k1112(9,mvsiz) ,k1113(9,mvsiz) ,
146 . k1114(9,mvsiz) ,k1115(9,mvsiz) ,
147 . k1116(9,mvsiz) ,k1117(9,mvsiz) ,k1118(9,mvsiz) ,
148 . k1119(9,mvsiz) ,k1120(9,mvsiz) ,
149 . k1212(9,mvsiz) ,k1213(9,mvsiz) ,k1214(9,mvsiz) ,
150 . k1215(9,mvsiz) ,k1216(9,mvsiz) ,
151 . k1217(9,mvsiz) ,k1218(9,mvsiz) ,k1219(9,mvsiz) ,
152 . k1220(9,mvsiz) ,k1313(9,mvsiz) ,
153 . k1314(9,mvsiz) ,k1315(9,mvsiz) ,k1316(9,mvsiz) ,
154 . k1317(9,mvsiz) ,k1318(9,mvsiz) ,
155 . k1319(9,mvsiz) ,k1320(9,mvsiz) ,k1414(9,mvsiz) ,
156 . k1415(9,mvsiz) ,k1416(9,mvsiz) ,
157 . k1417(9,mvsiz) ,k1418(9,mvsiz) ,k1419(9,mvsiz) ,
158 . k1420(9,mvsiz) ,k1515(9,mvsiz) ,
159 . k1516(9,mvsiz) ,k1517(9,mvsiz) ,k1518(9,mvsiz) ,
160 . k1519(9,mvsiz) ,k1520(9,mvsiz) ,
161 . k1616(9,mvsiz) ,k1617(9,mvsiz) ,k1618(9,mvsiz) ,
162 . k1619(9,mvsiz) ,k1620(9,mvsiz) ,
163 . k1717(9,mvsiz) ,k1718(9,mvsiz) ,k1719(9,mvsiz) ,
164 . k1720(9,mvsiz) ,k1818(9,mvsiz) ,
165 . k1819(9,mvsiz) ,k1820(9,mvsiz) ,k1919(9,mvsiz) ,
166 . k1920(9,mvsiz) ,k2020(9,mvsiz) ,
167 . bufmat(*)
168
169 INTEGER ETAG(*) , IDDL(*) ,NDOF(*) , IADK(*) ,JDIK(*)
170
171 my_real
172 . K_DIAG(*) , K_LT(*)
173 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
174C-----------------------------------------------
175C L o c a l V a r i a b l e s
176C-----------------------------------------------
177 INTEGER NIPMAX,NPE, NF1, NF2, IADBUF,IBID,
178 . j,i,n,k
179 parameter(nipmax=81)
180 parameter(npe=20)
181
182
183 INTEGER NGL(MVSIZ), IR, IS, IT, ILAY, NPTT ,NPTS ,NPTR, NLAY, IP,
184 . mxt(mvsiz),ikorth, nc(mvsiz,20), ngeo(mvsiz),nptg
185 my_real
186 . xx(mvsiz,20), yy(mvsiz,20), zz(mvsiz,20),hh(2,mvsiz),
187 . dnidr(mvsiz,npe),dnids(mvsiz,npe),dnidt(mvsiz,npe),
188 . rx(mvsiz,nipmax) , ry(mvsiz,nipmax) , rz(mvsiz,nipmax) ,
189 . sx(mvsiz,nipmax) , sy(mvsiz,nipmax) , sz(mvsiz,nipmax) ,
190 . tx(mvsiz,nipmax) , ty(mvsiz,nipmax) , tz(mvsiz,nipmax) ,
191 , e1x(mvsiz,nipmax) ,e2x(mvsiz,nipmax) ,e3x(mvsiz,nipmax) ,
192 . e1y(mvsiz,nipmax) ,e2y(mvsiz,nipmax) ,e3y(mvsiz,nipmax) ,
193 . e1z(mvsiz,nipmax) ,e2z(mvsiz,nipmax) ,e3z(mvsiz,nipmax) ,
194 . px(mvsiz,npe,nipmax),py(mvsiz,npe,nipmax),pz(mvsiz,npe,nipmax),
195 . off(mvsiz), wi, volnp(mvsiz,nipmax),stin(mvsiz,npe),
196 . ni(npe,nipmax), volg(mvsiz),ul(mvsiz,npe),gama(mvsiz,6),
197 . dm(9,mvsiz), dgm(9,mvsiz),gm(9,mvsiz), deltax(mvsiz),bid(1),
198 . dd(9,mvsiz),gg(mvsiz),dg(9,mvsiz),g33(9,mvsiz),offg(mvsiz)
199 TYPE(g_bufel_) ,POINTER :: GBUF
200 TYPE(L_BUFEL_) ,POINTER :: LBUF
201 DOUBLE PRECISION
202 . VOLDP(MVSIZ)
203C----------------
204
205 my_real
206 . w_gauss(9,9),a_gauss(9,9)
207 DATA w_gauss /
208 1 2. ,0. ,0. ,
209 1 0. ,0. ,0. ,
210 1 0. ,0. ,0. ,
211 2 1. ,1. ,0. ,
212 2 0. ,0. ,0. ,
213 2 0. ,0. ,0. ,
214 3 0.555555555555556,0.888888888888889,0.555555555555556,
215 3 0. ,0. ,0. ,
216 3 0. ,0. ,0. ,
217 4 0.347854845137454,0.652145154862546,0.652145154862546,
218 4 0.347854845137454,0. ,0. ,
219 4 0. ,0. ,0. ,
220 5 0.236926885056189,0.478628670499366,0.568888888888889,
221 5 0.478628670499366,0.236926885056189,0. ,
222 5 0. ,0. ,0. ,
223 6 0.171324492379170,0.360761573048139,0.467913934572691,
224 6 0.467913934572691,0.360761573048139,0.171324492379170,
225 6 0. ,0. ,0. ,
226 7 0.129484966168870,0.279705391489277,0.381830050505119,
227 7 0.417959183673469,0.381830050505119,0.279705391489277,
228 7 0.129484966168870,0. ,0. ,
229 8 0.101228536290376,0.222381034453374,0.313706645877887,
230 8 0.362683783378362,0.362683783378362,0.313706645877887,
231 8 0.222381034453374,0.101228536290376,0. ,
232 9 0.081274388361574,0.180648160694857,0.260610696402935,
233 9 0.312347077040003,0.330239355001260,0.312347077040003,
234 9 0.260610696402935,0.180648160694857,0.081274388361574/
235 DATA a_gauss /
236 1 0. ,0. ,0. ,
237 1 0. ,0. ,0. ,
238 1 0. ,0. ,0. ,
239 2 -.577350269189626,0.577350269189626,0. ,
240 2 0. ,0. ,0. ,
241 2 0. ,0. ,0. ,
242 3 -.774596669241483,0. ,0.774596669241483,
243 3 0. ,0. ,0. ,
244 3 0. ,0. ,0. ,
245 4 -.861136311594053,-.339981043584856,0.339981043584856,
246 4 0.861136311594053,0. ,0. ,
247 4 0. ,0. ,0. ,
248 5 -.906179845938664,-.538469310105683,0. ,
249 5 0.538469310105683,0.906179845938664,0. ,
250 5 0. ,0. ,0. ,
251 6 -.932469514203152,-.661209386466265,-.238619186083197,
252 6 0.238619186083197,0.661209386466265,0.932469514203152,
253 6 0. ,0. ,0. ,
254 7 -.949107912342759,-.741531185599394,-.405845151377397,
255 7 0. ,0.405845151377397,0.741531185599394,
256 7 0.949107912342759,0. ,0. ,
257 8 -.960289856497536,-.796666477413627,-.525532409916329,
258 8 -.183434642495650,0.183434642495650,0.525532409916329,
259 8 0.796666477413627,0.960289856497536,0. ,
260 9 -.968160239507626,-.836031107326636,-.613371432700590,
261 9 -.324253423403809,0. ,0.324253423403809,
262 9 0.613371432700590,0.836031107326636,0.968160239507626/
263
264C-----------------------------------------------
265C S o u r c e L i n e s
266C-----------------------------------------------
267 gbuf => elbuf_str%GBUF
268 nf1=nft+1
269 nf2=nf1-(numels8+numels10)
270
271
272 IF (igtyp == 21.OR.igtyp == 22) THEN
273 ikorth=2
274 ELSEIF (isorth>0) THEN
275 ikorth=1
276 ELSE
277 ikorth=0
278 ENDIF
279C-----------
280 CALL s20coork(
281 1 x, ixs(1,nf1), ixs20(1,nf2),npe,
282 2 xx, yy, zz, gbuf%OFF,
283 3 offg, gbuf%SMSTR, nc, ngl,
284 4 mxt, ngeo, nel, k11,
285 5 k12, k13, k14, k15,
286 6 k16, k17, k18, k19,
287 7 k110, k111, k112, k113,
288 8 k114, k115, k116, k117,
289 9 k118, k119, k120, k22,
290 a k23, k24, k25, k26,
291 b k27, k28, k29, k210,
292 c k211, k212, k213, k214,
293 d k215, k216, k217, k218,
294 e k219, k220, k33, k34,
295 f k35, k36, k37, k38,
296 g k39, k310, k311, k312,
297 h k313, k314, k315, k316,
298 i k317, k318, k319, k320,
299 j k44, k45, k46, k47,
300 k k48, k49, k410, k411,
301 l k412, k413, k414, k415,
302 m k416, k417, k418, k419,
303 n k420, k55, k56, k57,
304 o k58, k59, k510, k511,
305 p k512, k513, k514, k515,
306 q k516, k517, k518, k519,
307 r k520, k66, k67, k68,
308 s k69, k610, k611, k612,
309 t k613, k614, k615, k616,
310 u k617, k618, k619, k620,
311 v k77, k78, k79, k710,
312 w k711, k712, k713, k714,
313 x k715, k716, k717, k718,
314 y k719, k720, k88, k89,
315 z k810, k811, k812, k813,
316 1 k814, k815, k816, k817,
317 2 k818, k819, k820, k99,
318 3 k910, k911, k912, k913,
319 4 k914, k915, k916, k917,
320 5 k918, k919, k920, k1010,
321 6 k1011, k1012, k1013, k1014,
322 7 k1015, k1016, k1017, k1018,
323 8 k1019, k1020, k1111, k1112,
324 9 k1113, k1114, k1115, k1116,
325 a k1117, k1118, k1119, k1120,
326 b k1212, k1213, k1214, k1215,
327 c k1216, k1217, k1218, k1219,
328 d k1220, k1313, k1314, k1315,
329 e k1316, k1317, k1318, k1319,
330 f k1320, k1414, k1415, k1416,
331 g k1417, k1418, k1419, k1420,
332 h k1515, k1516, k1517, k1518,
333 i k1519, k1520, k1616, k1617,
334 j k1618, k1619, k1620, k1717,
335 k k1718, k1719, k1720, k1818,
336 l k1819, k1820, k1919, k1920,
337 m k2020, ismstr)
338C-----------------------------
339
340 DO n=1,npe
341 DO i=1,nel
342 ul(i,n) = zero
343 ENDDO
344 ENDDO
345 DO i=1,nel
346 volg(i) = zero
347 ENDDO
348
349C-----------------------------
350C POINTS D' INTEGRATION
351C-----------------------------
352 nlay = elbuf_str%NLAY
353 nptr = elbuf_str%NPTR
354 npts = elbuf_str%NPTS
355 nptt = elbuf_str%NPTT
356 nptg=nptt*npts*nptr*nlay
357 IF(nptg>nipmax)stop 933
358 DO it=1,nptt
359 DO is=1,npts
360 DO ir=1,nptr
361C-----------
362 ip = ir + ( (is-1) + (it-1)*npts )*nptr
363 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*w_gauss(it,nptt)
364 CALL s20rst(
365 1 a_gauss(ir,nptr),a_gauss(is,npts),a_gauss(it,nptt),ni(1,ip),
366 2 dnidr ,dnids ,dnidt )
367C
368C
369 CALL s20deri3(
370 1 ngl, offg, a_gauss(ir,nptr),a_gauss(is,npts),
371 2 a_gauss(it,nptt),wi, dnidr, dnids,
372 3 dnidt, rx(1,ip), ry(1,ip), rz(1,ip),
373 4 sx(1,ip), sy(1,ip), sz(1,ip), tx(1,ip),
374 5 ty(1,ip), tz(1,ip), xx, yy,
375 6 zz, px(1,1,ip), py(1,1,ip), pz(1,1,ip),
376 7 volnp(1,ip), deltax, stin, ni(1,ip),
377 8 volg, ul, ir, is,
378 9 it, voldp, nel)
379
380 CALL sreploc3(
381 1 rx(1,ip), ry(1,ip), rz(1,ip), sx(1,ip),
382 2 sy(1,ip), sz(1,ip), tx(1,ip), ty(1,ip),
383 3 tz(1,ip), e1x(1,ip),e2x(1,ip),e3x(1,ip),
384 4 e1y(1,ip),e2y(1,ip),e3y(1,ip),e1z(1,ip),
385 5 e2z(1,ip),e3z(1,ip),nel)
386 IF (isorth == 0) THEN
387 DO i=1,nel
388 gama(i,1) = one
389 gama(i,2) = zero
390 gama(i,3) = zero
391 gama(i,4) = zero
392 gama(i,5) = one
393 gama(i,6) = zero
394 ENDDO
395 ELSE
396 CALL sorthdir3(
397 1 rx(1,ip), ry(1,ip), rz(1,ip), sx(1,ip),
398 2 sy(1,ip), sz(1,ip), tx(1,ip), ty(1,ip),
399 3 tz(1,ip), e1x(1,ip),e2x(1,ip),e3x(1,ip),
400 4 e1y(1,ip),e2y(1,ip),e3y(1,ip),e1z(1,ip),
401 5 e2z(1,ip),e3z(1,ip),gbuf%GAMA,gama,
402 6 nel, irep)
403
404 ENDIF
405C
406 ENDDO
407 ENDDO
408 ENDDO
409
410
411C-----------
412 IF (mtn>=28) THEN
413 iadbuf = ipm(7,mxt(1))
414 ELSE
415 iadbuf = 1
416 ENDIF
417
418 CALL mmats(1 ,nel ,pm ,mxt ,hh ,
419 . mtn ,ikorth ,ipm ,igeo ,gama ,
420 . bufmat(iadbuf) ,dm ,dgm ,gm ,
421 . jhbe ,gbuf%SIG ,bid ,nptg ,nel )
422
423
424 ibid = 0
425C-----------Begin integrating points-----
426 ilay = 1
427 DO ir=1,nptr
428 DO is=1,npts
429 DO it=1,nptt
430 lbuf => elbuf_str%BUFLY(ilay)%LBUF(it,is,it)
431C-----------
432 ip = ir + ( (is-1) + (it-1)*npts )*nptr
433 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*w_gauss(it,nptt)
434C
435 CALL mmstifs(
436 1 pm, mxt, hh, volnp(1,ip),
437 2 ibid, dd, gg, dg,
438 3 g33, dm, gm, dgm,
439 4 ikorth, gbuf%SIG, ir, is,
440 5 it, nel, jhbe, mtn)
441
442
443
444 CALL s20cumg3(
445 1 px(1,1,ip),py(1,1,ip),pz(1,1,ip),dd,
446 2 gg, dg, g33, ikorth,
447 3 k11, k12, k13, k14,
448 4 k15, k16, k17, k18,
449 5 k19, k110, k111, k112,
450 6 k113, k114, k115, k116,
451 7 k117, k118, k119, k120,
452 8 k22, k23, k24, k25,
453 9 k26, k27, k28, k29,
454 a k210, k211, k212, k213,
455 b k214, k215, k216, k217,
456 c k218, k219, k220, k33,
457 d k34, k35, k36, k37,
458 e k38, k39, k310, k311,
459 f k312, k313, k314, k315,
460 g k316, k317, k318, k319,
461 h k320, k44, k45, k46,
462 i k47, k48, k49, k410,
463 j k411, k412, k413, k414,
464 k k415, k416, k417, k418,
465 l k419, k420, k55, k56,
466 m k57, k58, k59, k510,
467 n k511, k512, k513, k514,
468 o k515, k516, k517, k518,
469 p k519, k520, k66, k67,
470 q k68, k69, k610, k611,
471 r k612, k613, k614, k615,
472 s k616, k617, k618, k619,
473 t k620, k77, k78, k79,
474 u k710, k711, k712, k713,
475 v k714, k715, k716, k717,
476 w k718, k719, k720, k88,
477 x k89, k810, k811, k812,
478 y k813, k814, k815, k816,
479 z k817, k818, k819, k820,
480 1 k99, k910, k911, k912,
481 2 k913, k914, k915, k916,
482 3 k917, k918, k919, k920,
483 4 k1010, k1011, k1012, k1013,
484 5 k1014, k1015, k1016, k1017,
485 6 k1018, k1019, k1020, k1111,
486 7 k1112, k1113, k1114, k1115,
487 8 k1116, k1117, k1118, k1119,
488 9 k1120, k1212, k1213, k1214,
489 a k1215, k1216, k1217, k1218,
490 b k1219, k1220, k1313, k1314,
491 c k1315, k1316, k1317, k1318,
492 d k1319, k1320, k1414, k1415,
493 e k1416, k1417, k1418, k1419,
494 f k1420, k1515, k1516, k1517,
495 g k1518, k1519, k1520, k1616,
496 h k1617, k1618, k1619, k1620,
497 i k1717, k1718, k1719, k1720,
498 j k1818, k1819, k1820, k1919,
499 k k1920, k2020, nel)
500
501 IF (ikgeo>0) THEN
502 CALL s20kgeo3(
503 1 lbuf%SIG, volnp(1,ip),px(1,1,ip), py(1,1,ip),
504 2 pz(1,1,ip), k11, k12, k13,
505 3 k14, k15, k16, k17,
506 4 k18, k19, k110, k111,
507 5 k112, k113, k114, k115,
508 6 k116, k117, k118, k119,
509 7 k120, k22, k23, k24,
510 8 k25, k26, k27, k28,
511 9 k29, k210, k211, k212,
512 a k213, k214, k215, k216,
513 b k217, k218, k219, k220,
514 c k33, k34, k35, k36,
515 d k37, k38, k39, k310,
516 e k311, k312, k313, k314,
517 f k315, k316, k317, k318,
518 g k319, k320, k44, k45,
519 h k46, k47, k48, k49,
520 i k410, k411, k412, k413,
521 j k414, k415, k416, k417,
522 k k418, k419, k420, k55,
523 l k56, k57, k58, k59,
524 m k510, k511, k512, k513,
525 n k514, k515, k516, k517,
526 o k518, k519, k520, k66,
527 p k67, k68, k69, k610,
528 q k611, k612, k613, k614,
529 r k615, k616, k617, k618,
530 s k619, k620, k77, k78,
531 t k79, k710, k711, k712,
532 u k713, k714, k715, k716,
533 v k717, k718, k719, k720,
534 w k88, k89, k810, k811,
535 x k812, k813, k814, k815,
536 y k816, k817, k818, k819,
537 z k820, k99, k910, k911,
538 1 k912, k913, k914, k915,
539 2 k916, k917, k918, k919,
540 3 k920, k1010, k1011, k1012,
541 4 k1013, k1014, k1015, k1016,
542 5 k1017, k1018, k1019, k1020,
543 6 k1111, k1112, k1113, k1114,
544 7 k1115, k1116, k1117, k1118,
545 8 k1119, k1120, k1212, k1213,
546 9 k1214, k1215, k1216, k1217,
547 a k1218, k1219, k1220, k1313,
548 b k1314, k1315, k1316, k1317,
549 c k1318, k1319, k1320, k1414,
550 d k1415, k1416, k1417, k1418,
551 e k1419, k1420, k1515, k1516,
552 f k1517, k1518, k1519, k1520,
553 g k1616, k1617, k1618, k1619,
554 h k1620, k1717, k1718, k1719,
555 i k1720, k1818, k1819, k1820,
556 j k1919, k1920, k2020, nel)
557 ENDIF
558
559 ENDDO
560 ENDDO
561 ENDDO
562
563
564 IF (neig>0) CALL s20eoff(
565 1 1, nel, ixs(1,nf1),ixs20(1,nf2), etag, offg)
566 CALL assem_s20(
567 1 ixs(1,nf1),ixs20(1,nf2),nel,iddl ,ndof ,
568 2 k_diag,k_lt ,iadk ,jdik ,offg ,
569 . k11 ,k12 ,k13 ,k14 ,k15 ,k16 ,k17 ,k18 ,k19 ,k110 ,
570 . k111 ,k112 ,k113 ,k114 ,k115 ,k116 ,k117 ,k118 ,k119 ,k120 ,
571 . k22 ,k23 ,k24 ,k25 ,k26 ,k27 ,k28 ,k29 ,k210 ,k211 ,
572 . k212 ,k213 ,k214 ,k215 ,k216 ,k217 ,k218 ,k219 ,k220 ,k33 ,
573 . k34 ,k35 ,k36 ,k37 ,k38 ,k39 ,k310 ,k311 ,k312 ,k313 ,
574 . k314 ,k315 ,k316 ,k317 ,k318 ,k319 ,k320 ,k44 ,k45 ,k46 ,
575 . k47 ,k48 ,k49 ,k410 ,k411 ,k412 ,k413 ,k414 ,k415 ,k416 ,
576 . k417 ,k418 ,k419 ,k420 ,k55 ,k56 ,k57 ,k58 ,k59 ,k510 ,
577 . k511 ,k512 ,k513 ,k514 ,k515 ,k516 ,k517 ,k518 ,k519 ,k520 ,
578 . k66 ,k67 ,k68 ,k69 ,k610 ,k611 ,k612 ,k613 ,k614 ,k615 ,
579 . k616 ,k617 ,k618 ,k619 ,k620 ,k77 ,k78 ,k79 ,k710 ,k711 ,
580 . k712 ,k713 ,k714 ,k715 ,k716 ,k717 ,k718 ,k719 ,k720 ,k88 ,
581 . k89 ,k810 ,k811 ,k812 ,k813 ,k814 ,k815 ,k816 ,k817 ,k818 ,
582 . k819 ,k820 ,k99 ,k910 ,k911 ,k912 ,k913 ,k914 ,k915 ,k916 ,
583 . k917 ,k918 ,k919 ,k920 ,k1010,k1011,k1012,k1013,k1014,k1015,
584 . k1016,k1017,k1018,k1019,k1020,k1111,k1112,k1113,k1114,k1115,
585 . k1116,k1117,k1118,k1119,k1120,k1212,k1213,k1214,k1215,k1216,
586 . k1217,k1218,k1219,k1220,k1313,k1314,k1315,k1316,k1317,k1318,
587 . k1319,k1320,k1414,k1415,k1416,k1417,k1418,k1419,k1420,k1515,
588 . k1516,k1517,k1518,k1519,k1520,k1616,k1617,k1618,k1619,k1620,
589 . k1717,k1718,k1719,k1720,k1818,k1819,k1820,k1919,k1920,k2020)
590
591
592C
593 RETURN
594 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
my_real
#define my_real
Definition cppsort.cpp:32
mmats
subroutine mmats(jft, jlt, pm, mat, hh, mtn, iorth, ipm, igeo, gama, uparam, cc, cg, g33, jhbe, sig, eps, nppt, nel)
Definition mmats.F:45
mmstifs
subroutine mmstifs(pm, mat, hh, vol, icsig, dd, gg, dg, g33, dm, gm, dgm, iorth, sig, ir, is, it, nel, jhbe, mtn)
Definition mmstifs.F:44
s20coork
subroutine s20coork(x, ixs, ixs20, npe, xx, yy, zz, offg, off, sav, nc, ngl, mxt, ngeo, nel, 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, ismstr)
Definition s20coork.F:86
s20cumg3
subroutine s20cumg3(px, py, pz, dd, gg, dg, g33, iksup, 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, nel)
Definition s20cumg3.F:89
s20eoff
subroutine s20eoff(jft, jlt, ixs, ixs20, etag, off)
Definition s20eoff.F:29
s20ke3
subroutine s20ke3(pm, geo, ixs, ixs20, x, elbuf_str, etag, iddl, ndof, k_diag, k_lt, iadk, jdik, nel, ipm, igeo, ikgeo, bufmat, nft, mtn, ismstr, jhbe, irep, igtyp, isorth)
Definition s20ke3.F:51
s20kgeo3
subroutine s20kgeo3(sig, vol, px, py, pz, 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, nel)
Definition s20kgeo3.F:86
sorthdir3
subroutine sorthdir3(rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, gama0, gama, nel, irep)
Definition sorthdir3.F:42
s20deri3
subroutine s20deri3(ngl, off, r, s, t, w, dnidr, dnids, dnidt, dxdr, dydr, dzdr, dxds, dyds, dzds, dxdt, dydt, dzdt, xx, yy, zz, px, py, pz, vol, deltax, deltax2, ip, nip, ul, volg, voldp)
Definition s20deri3.F:38
s20rst
subroutine s20rst(r, s, t, ni, dnidr, dnids, dnidt)
Definition s20deri3.F:505
sreploc3
subroutine sreploc3(rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition sreploc3.F:32