OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8zke3.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!|| s8zke3 ../engine/source/elements/solid/solide8z/s8zke3.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_s8 ../engine/source/implicit/assem_s8.F
30!|| mmat_h1 ../engine/source/elements/solid/solide8z/mmat_h1.F
31!|| mmats ../engine/source/elements/solid/solide8z/mmats.F
32!|| mmstifs ../engine/source/elements/solid/solide8z/mmstifs.F
33!|| s8eoff ../engine/source/elements/solid/solide8/s8eoff.F
34!|| s8zcumg3 ../engine/source/elements/solid/solide8z/s8zcumg3.F
35!|| s8zderict3 ../engine/source/elements/solid/solide8z/s8zderict3.F
36!|| s8zderii3 ../engine/source/elements/solid/solide8z/s8zderii3.f
37!|| s8zkgei3 ../engine/source/elements/solid/solide8z/s8zkgei3.F
38!|| s8zkgeo3 ../engine/source/elements/solid/solide8z/s8zkgeo3.F
39!|| s8zlke3 ../engine/source/elements/solid/solide8z/s8zlke3.F
40!|| s8zlkep3 ../engine/source/elements/solid/solide8z/s8zlkep3.F
41!|| s8zsigp3 ../engine/source/elements/solid/solide8z/s8zsigp3.F
42!|| srcoork ../engine/source/elements/solid/solide8z/srcoork.F
43!||--- uses -----------------------------------------------------
44!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
45!|| element_mod ../common_source/modules/elements/element_mod.F90
46!||====================================================================
47 SUBROUTINE s8zke3(
48 1 PM, GEO, IXS, X,
49 2 ELBUF_STR,NEL, ICP, ICSIG,
50 3 ETAG, IDDL, NDOF, K_DIAG,
51 4 K_LT, IADK, JDIK, MPT,
52 5 IPM, IGEO, IKGEO, BUFMAT,
53 6 NFT, MTN, ISMSTR, JHBE,
54 7 JCVT, IGTYP, ISORTH, IREP)
55C-----------------------------------------------
56C M o d u l e s
57C-----------------------------------------------
58 USE elbufdef_mod
59 use element_mod , only : nixs
60C-----------------------------------------------
61C I m p l i c i t T y p e s
62C-----------------------------------------------
63#include "implicit_f.inc"
64C-----------------------------------------------
65C G l o b a l P a r a m e t e r s
66C-----------------------------------------------
67#include "mvsiz_p.inc"
68C-----------------------------------------------
69C C o m m o n B l o c k s
70C-----------------------------------------------
71#include "com04_c.inc"
72#include "param_c.inc"
73C-----------------------------------------------
74C D u m m y A r g u m e n t s
75C-----------------------------------------------
76 INTEGER, INTENT(IN) :: IREP
77 INTEGER, INTENT(IN) :: NFT
78 INTEGER, INTENT(IN) :: MTN
79 INTEGER, INTENT(IN) :: ISMSTR
80 INTEGER, INTENT(IN) :: JHBE
81 INTEGER, INTENT(IN) :: JCVT
82 INTEGER, INTENT(IN) :: IGTYP
83 INTEGER, INTENT(IN) :: ISORTH
84 INTEGER NEL ,ICP, ICSIG,IKGEO,MPT
85 INTEGER NDOF(*) ,IADK(*) ,JDIK(*),
86 . ixs(nixs,*),ipm(npropmi,*),igeo(npropgi,*),etag(*),iddl(*)
87C
88 my_real
89 . pm(npropm,*), geo(npropg,*), x(*),
90 . k11(9,mvsiz) , k12(9,mvsiz) ,k13(9,mvsiz) ,k14(9,mvsiz) ,
91 . k15(9,mvsiz) , k16(9,mvsiz) ,k17(9,mvsiz) ,k18(9,mvsiz) ,
92 . k22(9,mvsiz) ,k23(9,mvsiz) ,k24(9,mvsiz) ,k25(9,mvsiz) ,
93 . k26(9,mvsiz) ,k27(9,mvsiz) ,k28(9,mvsiz) ,k33(9,mvsiz) ,
94 . k34(9,mvsiz) ,k35(9,mvsiz) ,k36(9,mvsiz) ,k37(9,mvsiz) ,
95 . k38(9,mvsiz) ,k44(9,mvsiz) ,k45(9,mvsiz) ,k46(9,mvsiz) ,
96 . k47(9,mvsiz) ,k48(9,mvsiz) ,k55(9,mvsiz) ,k56(9,mvsiz) ,
97 . k57(9,mvsiz) ,k58(9,mvsiz) ,k66(9,mvsiz) ,k67(9,mvsiz) ,
98 . k68(9,mvsiz) ,k77(9,mvsiz) ,k78(9,mvsiz) ,k88(9,mvsiz) ,
99 . offg(mvsiz) ,bufmat(*) ,k_diag(*) ,k_lt(*)
100 TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
101C-----------------------------------------------
102C L o c a l V a r i a b l e s
103C-----------------------------------------------
104 INTEGER NF1, I, IKORTH, IADBUF, IPREDU,
105 . IR, IS, IT,IL,IP,ICPG,NPTR,NPTS,NPTT,NLAY,NPPT,
106 . IAD0,PID
107 INTEGER MXT(MVSIZ),NGL(MVSIZ),NGEO(MVSIZ),IDEG(MVSIZ)
108 my_real
109 . VOLN(MVSIZ), DELTAX(MVSIZ),
110 . AJC1(MVSIZ) , AJC2(MVSIZ) , AJC3(MVSIZ) ,
111 . ajc4(mvsiz) , ajc5(mvsiz) , ajc6(mvsiz) ,
112 . ajc7(mvsiz) , ajc8(mvsiz) , ajc9(mvsiz) ,
113 . aj1(mvsiz) , aj2(mvsiz) , aj3(mvsiz) ,
114 . aj4(mvsiz) , aj5(mvsiz) , aj6(mvsiz)
115C----------------
116 INTEGER NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ), NC4(MVSIZ),
117 . NC5(MVSIZ), NC6(MVSIZ), NC7(MVSIZ), NC8(MVSIZ)
118 my_real
119 . OFF(MVSIZ) ,
120 . X1(MVSIZ), X2(MVSIZ), X3(MVSIZ), X4(MVSIZ),
121 . X5(MVSIZ), X6(MVSIZ), X7(MVSIZ), X8(MVSIZ),
122 . Y1(MVSIZ), Y2(MVSIZ), Y3(MVSIZ), Y4(MVSIZ),
123 . y5(mvsiz), y6(mvsiz), y7(mvsiz), y8(mvsiz),
124 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
125 . z5(mvsiz), z6(mvsiz), z7(mvsiz), z8(mvsiz),
126 . px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),
127 . px5(mvsiz),px6(mvsiz),px7(mvsiz),px8(mvsiz),
128 . py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),
129 . py5(mvsiz),py6(mvsiz),py7(mvsiz),py8(mvsiz),
130 . pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),
131 . pz5(mvsiz),pz6(mvsiz),pz7(mvsiz),pz8(mvsiz),
132 . pxy1(mvsiz),pxy2(mvsiz),pxy3(mvsiz),pxy4(mvsiz),
133 . pxy5(mvsiz),pxy6(mvsiz),pxy7(mvsiz),pxy8(mvsiz),
134 . pyx1(mvsiz),pyx2(mvsiz),pyx3(mvsiz),pyx4(mvsiz),
135 . pyx5(mvsiz),pyx6(mvsiz),pyx7(mvsiz),pyx8(mvsiz),
136 . pxz1(mvsiz),pxz2(mvsiz),pxz3(mvsiz),pxz4(mvsiz),
137 . pxz5(mvsiz),pxz6(mvsiz),pxz7(mvsiz),pxz8(mvsiz),
138 . pzx1(mvsiz),pzx2(mvsiz),pzx3(mvsiz),pzx4(mvsiz),
139 . pzx5(mvsiz),pzx6(mvsiz),pzx7(mvsiz),pzx8(mvsiz),
140 . pyz1(mvsiz),pyz2(mvsiz),pyz3(mvsiz),pyz4(mvsiz),
141 . pyz5(mvsiz),pyz6(mvsiz),pyz7(mvsiz),pyz8(mvsiz),
142 . pzy1(mvsiz),pzy2(mvsiz),pzy3(mvsiz),pzy4(mvsiz),
143 . pzy5(mvsiz),pzy6(mvsiz),pzy7(mvsiz),pzy8(mvsiz),
144 . pxc1(mvsiz),pxc2(mvsiz),pxc3(mvsiz),pxc4(mvsiz),
145 . pyc1(mvsiz),pyc2(mvsiz),pyc3(mvsiz),pyc4(mvsiz),
146 . pzc1(mvsiz),pzc2(mvsiz),pzc3(mvsiz),pzc4(mvsiz),
147 . bxy1(mvsiz),bxy2(mvsiz),bxy3(mvsiz),bxy4(mvsiz),
148 . bxy5(mvsiz),bxy6(mvsiz),bxy7(mvsiz),bxy8(mvsiz),
149 . byx1(mvsiz),byx2(mvsiz),byx3(mvsiz),byx4(mvsiz),
150 . byx5(mvsiz),byx6(mvsiz),byx7(mvsiz),byx8(mvsiz),
151 . bxz1(mvsiz),bxz2(mvsiz),bxz3(mvsiz),bxz4(mvsiz),
152 . bxz5(mvsiz),bxz6(mvsiz),bxz7(mvsiz),bxz8(mvsiz),
153 . bzx1(mvsiz),bzx2(mvsiz),bzx3(mvsiz),bzx4(mvsiz),
154 . bzx5(mvsiz),bzx6(mvsiz),bzx7(mvsiz),bzx8(mvsiz),
155 . byz1(mvsiz),byz2(mvsiz),byz3(mvsiz),byz4(mvsiz),
156 . byz5(mvsiz),byz6(mvsiz),byz7(mvsiz),byz8(mvsiz),
157 . bzy1(mvsiz),bzy2(mvsiz),bzy3(mvsiz),bzy4(mvsiz),
158 . bzy5(mvsiz),bzy6(mvsiz),bzy7(mvsiz),bzy8(mvsiz),
159 . px1h1(mvsiz),px2h1(mvsiz),px3h1(mvsiz),px4h1(mvsiz),
160 . px1h2(mvsiz),px2h2(mvsiz),px3h2(mvsiz),px4h2(mvsiz),
161 . px1h3(mvsiz),px2h3(mvsiz),px3h3(mvsiz),px4h3(mvsiz),
162 . px1h4(mvsiz),px2h4(mvsiz),px3h4(mvsiz),px4h4(mvsiz),
163 . hx(mvsiz,4), hy(mvsiz,4), hz(mvsiz,4)
164 my_real
165 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
166 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
167 . r31(mvsiz),r32(mvsiz),r33(mvsiz)
168 my_real
169 . bid,wi,smax(mvsiz),volg(mvsiz),nu(mvsiz),
170 . nu1(mvsiz),gama(mvsiz,6)
171 TYPE(g_bufel_) ,POINTER :: GBUF
172 TYPE(L_BUFEL_) ,POINTER :: LBUF
173 my_real
174 . W_GAUSS(9,9),A_GAUSS(9,9)
175 DATA w_gauss /
176 1 2. ,0. ,0. ,
177 1 0. ,0. ,0. ,
178 1 0. ,0. ,0. ,
179 2 1. ,1. ,0. ,
180 2 0. ,0. ,0. ,
181 2 0. ,0. ,0. ,
182 3 0.555555555555556,0.888888888888889,0.555555555555556,
183 3 0. ,0. ,0. ,
184 3 0. ,0. ,0. ,
185 4 0.347854845137454,0.652145154862546,0.652145154862546,
186 4 0.347854845137454,0. ,0. ,
187 4 0. ,0. ,0. ,
188 5 0.236926885056189,0.478628670499366,0.568888888888889,
189 5 0.478628670499366,0.236926885056189,0. ,
190 5 0. ,0. ,0. ,
191 6 0.171324492379170,0.360761573048139,0.467913934572691,
192 6 0.467913934572691,0.360761573048139,0.171324492379170,
193 6 0. ,0. ,0. ,
194 7 0.129484966168870,0.279705391489277,0.381830050505119,
195 7 0.417959183673469,0.381830050505119,0.279705391489277,
196 7 0.129484966168870,0. ,0. ,
197 8 0.101228536290376,0.222381034453374,0.313706645877887,
198 8 0.362683783378362,0.362683783378362,0.313706645877887,
199 8 0.222381034453374,0.101228536290376,0. ,
200 9 0.081274388361574,0.180648160694857,0.260610696402935,
201 9 0.312347077040003,0.330239355001260,0.312347077040003,
202 9 0.260610696402935,0.180648160694857,0.081274388361574/
203 DATA a_gauss /
204 1 0. ,0. ,0. ,
205 1 0. ,0. ,0. ,
206 1 0. ,0. ,0. ,
207 2 -.577350269189626,0.577350269189626,0. ,
208 2 0. ,0. ,0. ,
209 2 0. ,0. ,0. ,
210 3 -.774596669241483,0. ,0.774596669241483,
211 3 0. ,0. ,0. ,
212 3 0. ,0. ,0. ,
213 4 -.861136311594053,-.339981043584856,0.339981043584856,
214 4 0.861136311594053,0. ,0. ,
215 4 0. ,0. ,0. ,
216 5 -.906179845938664,-.538469310105683,0. ,
217 5 0.538469310105683,0.906179845938664,0. ,
218 5 0. ,0. ,0. ,
219 6 -.932469514203152,-.661209386466265,-.238619186083197,
220 6 0.238619186083197,0.661209386466265,0.932469514203152,
221 6 0. ,0. ,0. ,
222 7 -.949107912342759,-.741531185599394,-.405845151377397,
223 7 0. ,0.405845151377397,0.741531185599394,
224 7 0.949107912342759,0. ,0. ,
225 8 -.960289856497536,-.796666477413627,-.525532409916329,
226 8 -.183434642495650,0.183434642495650,0.525532409916329,
227 8 0.796666477413627,0.960289856497536,0. ,
228 9 -.968160239507626,-.836031107326636,-.613371432700590,
229 9 -.324253423403809,0. ,0.324253423403809,
230 9 0.613371432700590,0.836031107326636,0.968160239507626/
231
232 my_real
233 . hh(2,mvsiz),fac(mvsiz),c1,e0(mvsiz),dd(9,mvsiz),gg(mvsiz),
234 . dm(9,mvsiz),gm(9,mvsiz),dgm(9,mvsiz),dg(9,mvsiz),
235 . g33(9,mvsiz),hh1(2,mvsiz)
236C----HH(1,):lamda,HH(2,):G; HH1 :effective lamda,g for Icpre=1
237C-----------------------------------------------
238C S o u r c e L i n e s
239C-----------------------------------------------
240 gbuf => elbuf_str%GBUF
241 nlay = elbuf_str%NLAY
242 IF (mpt == 222) THEN
243 nptr = 2
244 npts = 2
245 nptt = 2
246 ELSE
247 nptr = elbuf_str%NPTR
248 npts = elbuf_str%NPTS
249 nptt = elbuf_str%NPTT
250 ENDIF
251 iad0 = 1
252 IF (isorth > 0) iad0 = 1 + 6*nel
253 IF (igtyp == 21.OR.igtyp == 22) THEN
254 ikorth=2
255 ELSEIF (isorth>0) THEN
256 ikorth=1
257 ELSE
258 ikorth=0
259 ENDIF
260C-----------
261 nf1=nft+1
262C-----------
263 CALL srcoork(
264 1 x, ixs(1,nf1),x1, x2,
265 2 x3, x4, x5, x6,
266 3 x7, x8, y1, y2,
267 4 y3, y4, y5, y6,
268 5 y7, y8, z1, z2,
269 6 z3, z4, z5, z6,
270 7 z7, z8, gbuf%OFF, offg,
271 8 gbuf%SMSTR,r11, r12, r13,
272 9 r21, r22, r23, r31,
273 a r32, r33, nc1, nc2,
274 b nc3, nc4, nc5, nc6,
275 c nc7, nc8, ngl, mxt,
276 d ngeo, k11, k12, k13,
277 e k14, k15, k16, k17,
278 f k18, k22, k23, k24,
279 g k25, k26, k27, k28,
280 h k33, k34, k35, k36,
281 i k37, k38, k44, k45,
282 j k46, k47, k48, k55,
283 k k56, k57, k58, k66,
284 l k67, k68, k77, k78,
285 m k88, jhbe, gbuf%GAMA, gama,
286 n nel, ismstr, jhbe, jcvt,
287 o irep, igtyp, isorth)
288C
289 pid=ngeo(1)
290C IGTYP=NINT(GEO(12,PID))
291C CALL GET_GAMA(MTN,NEL,NBGAMA,JCVT,PM,IXS(1,NF1),IGTYP)
292 DO i=1,nel
293 nu(i)=min(half,pm(21,mxt(i)))
294 c1 =pm(32,mxt(i))
295 e0(i) =three*(one-two*nu(i))*c1
296C USB(I)=EM01/C1
297 deltax(i)=ep30
298 ideg(i)=0
299 ENDDO
300 IF (jhbe/=14.AND.jhbe<20) THEN
301 DO i=1,nel
302 nu(i)=zero
303 ENDDO
304 ELSE
305 ENDIF
306 icpg = icp
307 IF (jhbe/=14.AND.jhbe/=17.AND.jhbe/=24) THEN
308 icpg = 1
309 ENDIF
310c
311 CALL s8zderict3(
312 1 offg, volg, ngl, x1,
313 2 x2, x3, x4, x5,
314 3 x6, x7, x8, y1,
315 4 y2, y3, y4, y5,
316 5 y6, y7, y8, z1,
317 6 z2, z3, z4, z5,
318 7 z6, z7, z8, pxc1,
319 8 pxc2, pxc3, pxc4, pyc1,
320 9 pyc2, pyc3, pyc4, pzc1,
321 a pzc2, pzc3, pzc4, px1h1,
322 b px1h2, px1h3, px1h4, px2h1,
323 c px2h2, px2h3, px2h4, px3h1,
324 d px3h2, px3h3, px3h4, px4h1,
325 e px4h2, px4h3, px4h4, hx,
326 f hy, hz, ajc1, ajc2,
327 g ajc3, ajc4, ajc5, ajc6,
328 h ajc7, ajc8, ajc9, smax,
329 i nel)
330 IF (icpg==2) THEN
331 CALL s8zsigp3(1 ,nel ,gbuf%SIG,e0 ,gbuf%PLA,
332 2 fac ,gbuf%G_PLA,nel )
333 DO i=1,nel
334 nu1(i)=nu(i)+(half-nu(i))*fac(i)
335 ENDDO
336 ELSEIF (jhbe==17.AND.icpg==1.AND.ismstr==10) THEN
337 DO i=1,nel
338 fac(i)=two*nu(i)
339 nu1(i)=nu(i)
340 ENDDO
341 ELSE
342 DO i=1,nel
343 nu1(i)=nu(i)
344 ENDDO
345 ENDIF
346 IF (mtn>=28) THEN
347 iadbuf = ipm(7,mxt(1))
348 ELSE
349 iadbuf = 1
350 ENDIF
351 nppt = npts*nptt*nptr*nlay
352 CALL mmats(1 ,nel ,pm ,mxt ,hh ,
353 . mtn ,ikorth ,ipm ,igeo ,gama ,
354 . bufmat(iadbuf) ,dm ,dgm ,gm ,
355 . jhbe ,gbuf%SIG ,bid ,nppt ,nel )
356 CALL mmat_h1(
357 1 hh, hh1, fac, icpg,
358 2 ipredu, nel, mtn, ismstr,
359 3 jhbe)
360
361C-----------Matrix B----------
362C-----------Begin integrating points-----
363 il = 1
364 DO ir=1,nptr
365 DO is=1,npts
366 DO it=1,nptt
367 IF (jhbe == 14.OR.jhbe == 12.OR.jhbe == 17) THEN
368 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
369 ELSE
370 lbuf => elbuf_str%BUFLY(il)%LBUF(1,1,1)
371 ENDIF
372C-----------
373 ip = ir + ( (is-1) + (it-1)*npts )*nptr
374 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*w_gauss(it,nptt)
375C
376 CALL s8zderii3(
377 1 offg, off, voln, ngl,
378 2 a_gauss(ir,nptr),a_gauss(is,npts),a_gauss(it,nptt),wi,
379 3 pxc1, pxc2, pxc3, pxc4,
380 4 pyc1, pyc2, pyc3, pyc4,
381 5 pzc1, pzc2, pzc3, pzc4,
382 6 px1h1, px1h2, px1h3, px1h4,
383 7 px2h1, px2h2, px2h3, px2h4,
384 8 px3h1, px3h2, px3h3, px3h4,
385 9 px4h1, px4h2, px4h3, px4h4,
386 a hx, hy, hz, px1,
387 b px2, px3, px4, px5,
388 c px6, px7, px8, py1,
389 d py2, py3, py4, py5,
390 e py6, py7, py8, pz1,
391 f pz2, pz3, pz4, pz5,
392 g pz6, pz7, pz8, pxy1,
393 h pxy2, pxy3, pxy4, pxy5,
394 i pxy6, pxy7, pxy8, pyx1,
395 j pyx2, pyx3, pyx4, pyx5,
396 k pyx6, pyx7, pyx8, pxz1,
397 l pxz2, pxz3, pxz4, pxz5,
398 m pxz6, pxz7, pxz8, pzx1,
399 n pzx2, pzx3, pzx4, pzx5,
400 o pzx6, pzx7, pzx8, pyz1,
401 p pyz2, pyz3, pyz4, pyz5,
402 q pyz6, pyz7, pyz8, pzy1,
403 r pzy2, pzy3, pzy4, pzy5,
404 s pzy6, pzy7, pzy8, bxy1,
405 t bxy2, bxy3, bxy4, bxy5,
406 u bxy6, bxy7, bxy8, byx1,
407 v byx2, byx3, byx4, byx5,
408 w byx6, byx7, byx8, bxz1,
409 x bxz2, bxz3, bxz4, bxz5,
410 y bxz6, bxz7, bxz8, bzx1,
411 z bzx2, bzx3, bzx4, bzx5,
412 1 bzx6, bzx7, bzx8, byz1,
413 2 byz2, byz3, byz4, byz5,
414 3 byz6, byz7, byz8, bzy1,
415 4 bzy2, bzy3, bzy4, bzy5,
416 5 bzy6, bzy7, bzy8, ajc1,
417 6 ajc2, ajc3, ajc4, ajc5,
418 7 ajc6, ajc7, ajc8, ajc9,
419 8 aj1, aj2, aj3, aj4,
420 9 aj5, aj6, smax, deltax,
421 a nu1, icpg, ideg, nel,
422 b mtn, ismstr, jhbe)
423 CALL mmstifs(
424 1 pm, mxt, hh1, voln,
425 2 icsig, dd, gg, dg,
426 3 g33, dm, gm, dgm,
427 4 ikorth, lbuf%SIG,ir, is,
428 5 it, nel, jhbe, mtn)
429Ccheck mpt are the same between ke and forc3
430 CALL s8zlke3(
431 1 px1, px2, px3, px4,
432 2 py1, py2, py3, py4,
433 3 pz1, pz2, pz3, pz4,
434 4 px5, px6, px7, px8,
435 5 py5, py6, py7, py8,
436 6 pz5, pz6, pz7, pz8,
437 7 pxy1, pxy2, pxy3, pxy4,
438 8 pxy5, pxy6, pxy7, pxy8,
439 9 pyx1, pyx2, pyx3, pyx4,
440 a pyx5, pyx6, pyx7, pyx8,
441 b pxz1, pxz2, pxz3, pxz4,
442 c pxz5, pxz6, pxz7, pxz8,
443 d pzx1, pzx2, pzx3, pzx4,
444 e pzx5, pzx6, pzx7, pzx8,
445 f pyz1, pyz2, pyz3, pyz4,
446 g pyz5, pyz6, pyz7, pyz8,
447 h pzy1, pzy2, pzy3, pzy4,
448 i pzy5, pzy6, pzy7, pzy8,
449 j bxy1, bxy2, bxy3, bxy4,
450 k bxy5, bxy6, bxy7, bxy8,
451 l byx1, byx2, byx3, byx4,
452 m byx5, byx6, byx7, byx8,
453 n bxz1, bxz2, bxz3, bxz4,
454 o bxz5, bxz6, bxz7, bxz8,
455 p bzx1, bzx2, bzx3, bzx4,
456 q bzx5, bzx6, bzx7, bzx8,
457 r byz1, byz2, byz3, byz4,
458 s byz5, byz6, byz7, byz8,
459 t bzy1, bzy2, bzy3, bzy4,
460 u bzy5, bzy6, bzy7, bzy8,
461 v k11, k12, k13, k14,
462 w k15, k16, k17, k18,
463 x k22, k23, k24, k25,
464 y k26, k27, k28, k33,
465 z k34, k35, k36, k37,
466 1 k38, k44, k45, k46,
467 2 k47, k48, k55, k56,
468 3 k57, k58, k66, k67,
469 4 k68, k77, k78, k88,
470 5 dd, gg, icpg, offg,
471 6 ikorth, dg, g33, nel)
472 IF (ikgeo<0.AND.jhbe==14) THEN
473 CALL s8zkgei3(
474 1 lbuf%SIG,voln, px1, px2,
475 2 px3, px4, py1, py2,
476 3 py3, py4, pz1, pz2,
477 4 pz3, pz4, px5, px6,
478 5 px7, px8, py5, py6,
479 6 py7, py8, pz5, pz6,
480 7 pz7, pz8, k11, k12,
481 8 k13, k14, k15, k16,
482 9 k17, k18, k22, k23,
483 a k24, k25, k26, k27,
484 b k28, k33, k34, k35,
485 c k36, k37, k38, k44,
486 d k45, k46, k47, k48,
487 e k55, k56, k57, k58,
488 f k66, k67, k68, k77,
489 g k78, k88, nel)
490 ENDIF
491 ENDDO
492 ENDDO
493 ENDDO
494C-------------
495 IF (ipredu > 0) THEN
496 CALL s8zlkep3(
497 1 pxc1, pxc2, pxc3, pxc4,
498 2 pyc1, pyc2, pyc3, pyc4,
499 3 pzc1, pzc2, pzc3, pzc4,
500 4 k11, k12, k13, k14,
501 5 k15, k16, k17, k18,
502 6 k22, k23, k24, k25,
503 7 k26, k27, k28, k33,
504 8 k34, k35, k36, k37,
505 9 k38, k44, k45, k46,
506 a k47, k48, k55, k56,
507 b k57, k58, k66, k67,
508 c k68, k77, k78, k88,
509 d hh, volg, fac, icpg,
510 e nel)
511 END IF
512C------due to the generic form, only the mean geometrical stiffness is considered
513 IF (ikgeo>0) THEN
514 CALL s8zkgeo3(
515 1 gbuf%SIG,volg, pxc1, pxc2,
516 2 pxc3, pxc4, pyc1, pyc2,
517 3 pyc3, pyc4, pzc1, pzc2,
518 4 pzc3, pzc4, k11, k12,
519 5 k13, k14, k15, k16,
520 6 k17, k18, k22, k23,
521 7 k24, k25, k26, k27,
522 8 k28, k33, k34, k35,
523 9 k36, k37, k38, k44,
524 a k45, k46, k47, k48,
525 b k55, k56, k57, k58,
526 c k66, k67, k68, k77,
527 d k78, k88, nel)
528 ENDIF
529C----------------------------
530C CONVECTE --> GLOBAL.
531C----------------------------
532 CALL s8zcumg3(
533 1 r11, r21, r31, r12,
534 2 r22, r32, r13, r23,
535 3 r33, k11, k12, k13,
536 4 k14, k15, k16, k17,
537 5 k18, k22, k23, k24,
538 6 k25, k26, k27, k28,
539 7 k33, k34, k35, k36,
540 8 k37, k38, k44, k45,
541 9 k46, k47, k48, k55,
542 a k56, k57, k58, k66,
543 b k67, k68, k77, k78,
544 c k88, nel)
545
546 IF (neig>0) CALL s8eoff(
547 1 1, nel, ixs(1,nf1), etag, offg)
548 CALL assem_s8(
549 1 ixs(1,nf1),nel ,iddl ,ndof ,k_diag,
550 2 k_lt ,iadk ,jdik ,k11 ,k12 ,
551 3 k13 ,k14 ,k15 ,k16 ,k17 ,
552 4 k18 ,k22 ,k23 ,k24 ,k25 ,
553 5 k26 ,k27 ,k28 ,k33 ,k34 ,
554 6 k35 ,k36 ,k37 ,k38 ,k44 ,
555 7 k45 ,k46 ,k47 ,k48 ,k55 ,
556 8 k56 ,k57 ,k58 ,k66 ,k67 ,
557 9 k68 ,k77 ,k78 ,k88 ,offg )
558C
559 RETURN
560 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
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
mmat_h1
subroutine mmat_h1(hh, hh1, fac, icp, ipredu, nel, mtn, ismstr, jhbe)
Definition mmat_h1.F:34
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
s8eoff
subroutine s8eoff(jft, jlt, ixs, etag, off)
Definition s8eoff.F:35
s8zcumg3
subroutine s8zcumg3(r11, r21, r31, r12, r22, r32, r13, r23, r33, 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, nel)
Definition s8zcumg3.F:44
s8zderict3
subroutine s8zderict3(off, det, ngl, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px1h1, px1h2, px1h3, px1h4, px2h1, px2h2, px2h3, px2h4, px3h1, px3h2, px3h3, px3h4, px4h1, px4h2, px4h3, px4h4, hx, hy, hz, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, smax, nel)
Definition s8zderict3.F:55
s8zderii3
subroutine s8zderii3(offg, off, vol, ngl, ksi, eta, zeta, wi, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, px1h1, px1h2, px1h3, px1h4, px2h1, px2h2, px2h3, px2h4, px3h1, px3h2, px3h3, px3h4, px4h1, px4h2, px4h3, px4h4, hx, hy, hz, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, jac4, jac5, jac6, jac7, jac8, jac9, smax, deltax, nu, icp, ideg, nel, mtn, ismstr, jhbe)
Definition s8zderii3.F:75
s8zke3
subroutine s8zke3(pm, geo, ixs, x, elbuf_str, nel, icp, icsig, etag, iddl, ndof, k_diag, k_lt, iadk, jdik, mpt, ipm, igeo, ikgeo, bufmat, nft, mtn, ismstr, jhbe, jcvt, igtyp, isorth, irep)
Definition s8zke3.F:55
s8zkgei3
subroutine s8zkgei3(sig, vol, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px5, px6, px7, px8, py5, py6, py7, py8, pz5, pz6, pz7, pz8, 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, nel)
Definition s8zkgei3.F:48
s8zkgeo3
subroutine s8zkgeo3(sig, vol, 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, nel)
Definition s8zkgeo3.F:45
s8zlke3
subroutine s8zlke3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px5, px6, px7, px8, py5, py6, py7, py8, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, 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, dd, gg, icp, off, iksup, dg, g33, nel)
Definition s8zlke3.F:77
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
s8zsigp3
subroutine s8zsigp3(lft, llt, sig, e0, defp, fac, g_pla, nel)
Definition s8zsigp3.F:37
srcoork
subroutine srcoork(x, ixs, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, offg, off, sav, r11, r12, r13, r21, r22, r23, r31, r32, r33, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, ngl, mxt, ngeo, 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, khbe, gama0, gama, nel, ismstr, jhbe, jcvt, irep, igtyp, isorth)
Definition srcoork.F:62