OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8_is17jac_i.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!|| s8_is17jac_i ../engine/source/elements/solid/solide8e/s8_is17jac_i.F
25!||--- called by ------------------------------------------------------
26!|| s8fupd11t12 ../engine/source/elements/solid/solide8e/s8fupd11t12.F
27!||--- calls -----------------------------------------------------
28!|| s8zjac_i3 ../engine/source/elements/solid/solide8e/s8zjac_i3.F
29!|| s8zjac_ic ../engine/source/elements/solid/solide8e/s8zjac_ic.F
30!||--- uses -----------------------------------------------------
31!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
32!||====================================================================
33 SUBROUTINE s8_is17jac_i(ELBUF_TAB,
34 . X1 ,X2 ,X3 ,X4 ,X5 ,X6 ,X7 ,X8 ,
35 . Y1 ,Y2 ,Y3 ,Y4 ,Y5 ,Y6 ,Y7 ,Y8 ,
36 . Z1 ,Z2 ,Z3 ,Z4 ,Z5 ,Z6 ,Z7 ,Z8 ,
37 . INDEX,NCH ,NEL )
38C-----------------------------------------------
39C M o d u l e s
40C-----------------------------------------------
41 USE elbufdef_mod
42C-----------------------------------------------
43C I m p l i c i t T y p e s
44C-----------------------------------------------
45#include "implicit_f.inc"
46C-----------------------------------------------
47C G l o b a l P a r a m e t e r s
48C-----------------------------------------------
49#include "mvsiz_p.inc"
50C-----------------------------------------------
51C D u m m y A r g u m e n t s
52C-----------------------------------------------
53 INTEGER INDEX(*),NCH ,NEL
54 DOUBLE PRECISION
55 . X1(*), X2(*), X3(*), X4(*), X5(*), X6(*),
56 . X7(*), X8(*), Y1(*), Y2(*), Y3(*), Y4(*), Y5(*), Y6(*), Y7(*),
57 . y8(*), z1(*), z2(*), z3(*), z4(*), z5(*), z6(*), z7(*), z8(*)
58 TYPE(elbuf_struct_), TARGET :: ELBUF_TAB
59C-----------------------------------------------
60C L o c a l V a r i a b l e s
61C-----------------------------------------------
62 INTEGER I ,J ,N
63 INTEGER IP,IR, IS, IT,JHBE,IREP,MPT,NLAY,NPTR,NPTS,NPTT,NUVAR,
64 . L_PLA,L_SIGB,NSPHDIR, NCELF, NCELL,LLPIJ,L_JAC,NNPT
65C REAL
66 TYPE(l_bufel_) ,POINTER :: LBUF
67 TYPE(g_bufel_) ,POINTER :: GBUF
68 my_real
69 . AJC1(MVSIZ) , AJC2(MVSIZ) , AJC3(MVSIZ) ,
70 . ajc4(mvsiz) , ajc5(mvsiz) , ajc6(mvsiz) ,
71 . ajc7(mvsiz) , ajc8(mvsiz) , ajc9(mvsiz) ,
72 . hx(4,mvsiz) , hy(4,mvsiz), hz(4,mvsiz),
73 . pxc1(mvsiz),pxc2(mvsiz),pxc3(mvsiz),pxc4(mvsiz),
74 . pyc1(mvsiz),pyc2(mvsiz),pyc3(mvsiz),pyc4(mvsiz),
75 . pzc1(mvsiz),pzc2(mvsiz),pzc3(mvsiz),pzc4(mvsiz)
76 my_real
77 . wt,zr,zs,zt,zz,wi
78C-----------------------------------------------
79 my_real
80 . w_gauss(9,9),a_gauss(9,9)
81 DATA w_gauss /
82c---
83 1 2.d0 ,0.d0 ,0.d0 ,
84 1 0.d0 ,0.d0 ,0.d0 ,
85 1 0.d0 ,0.d0 ,0.d0 ,
86 2 1.d0 ,1.d0 ,0.d0 ,
87 2 0.d0 ,0.d0 ,0.d0 ,
88 2 0.d0 ,0.d0 ,0.d0 ,
89 3 0.555555555555556d0,0.888888888888889d0,0.555555555555556d0,
90 3 0.d0 ,0.d0 ,0.d0 ,
91 3 0.d0 ,0.d0 ,0.d0 ,
92 4 0.347854845137454d0,0.652145154862546d0,0.652145154862546d0,
93 4 0.347854845137454d0,0.d0 ,0.d0 ,
94 4 0.d0 ,0.d0 ,0.d0 ,
95 5 0.236926885056189d0,0.478628670499366d0,0.568888888888889d0,
96 5 0.478628670499366d0,0.236926885056189d0,0.d0 ,
97 5 0.d0 ,0.d0 ,0.d0 ,
98 6 0.171324492379170d0,0.360761573048139d0,0.467913934572691d0,
99 6 0.467913934572691d0,0.360761573048139d0,0.171324492379170d0,
100 6 0.d0 ,0.d0 ,0.d0 ,
101 7 0.129484966168870d0,0.279705391489277d0,0.381830050505119d0,
102 7 0.417959183673469d0,0.381830050505119d0,0.279705391489277d0,
103 7 0.129484966168870d0,0.d0 ,0.d0 ,
104 8 0.101228536290376d0,0.222381034453374d0,0.313706645877887d0,
105 8 0.362683783378362d0,0.362683783378362d0,0.313706645877887d0,
106 8 0.222381034453374d0,0.101228536290376d0,0.d0 ,
107 9 0.081274388361574d0,0.180648160694857d0,0.260610696402935d0,
108 9 0.312347077040003d0,0.330239355001260d0,0.312347077040003d0,
109 9 0.260610696402935d0,0.180648160694857d0,0.081274388361574d0/
110c------------------------------------------------------------
111 DATA a_gauss /
112 1 0.d0 ,0.d0 ,0.d0 ,
113 1 0.d0 ,0.d0 ,0.d0 ,
114 1 0.d0 ,0.d0 ,0.d0 ,
115 2 -.577350269189625d0,0.577350269189625d0,0.d0 ,
116 2 0.d0 ,0.d0 ,0.d0 ,
117 2 0.d0 ,0.d0 ,0.d0 ,
118 3 -.774596669241483d0,0.d0 ,0.774596669241483d0,
119 3 0.d0 ,0.d0 ,0.d0 ,
120 3 0.d0 ,0.d0 ,0.d0 ,
121 4 -.861136311594053d0,-.339981043584856d0,0.339981043584856d0,
122 4 0.861136311594053d0,0.d0 ,0.d0 ,
123 4 0.d0 ,0.d0 ,0.d0 ,
124 5 -.906179845938664d0,-.538469310105683d0,0.d0 ,
125 5 0.538469310105683d0,0.906179845938664d0,0.d0 ,
126 5 0.d0 ,0.d0 ,0.d0 ,
127 6 -.932469514203152d0,-.661209386466265d0,-.238619186083197d0,
128 6 0.238619186083197d0,0.661209386466265d0,0.932469514203152d0,
129 6 0.d0 ,0.d0 ,0.d0 ,
130 7 -.949107912342759d0,-.741531185599394d0,-.405845151377397d0,
131 7 0.d0 ,0.405845151377397d0,0.741531185599394d0,
132 7 0.949107912342759d0,0.d0 ,0.d0 ,
133 8 -.960289856497536d0,-.796666477413627d0,-.525532409916329d0,
134 8 -.183434642495650d0,0.183434642495650d0,0.525532409916329d0,
135 8 0.796666477413627d0,0.960289856497536d0,0.d0 ,
136 9 -.968160239507626d0,-.836031107326636d0,-.613371432700590d0,
137 9 -.324253423403809d0,0.d0 ,0.324253423403809d0,
138 9 0.613371432700590d0,0.836031107326636d0,0.968160239507626d0/
139C
140C-----------------------------------------------
141 gbuf => elbuf_tab%GBUF
142C cas GBUF%JAC_I for all case
143 CALL s8zjac_ic(
144 . x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,x7 ,x8 ,
145 . y1 ,y2 ,y3 ,y4 ,y5 ,y6 ,y7 ,y8 ,
146 . z1 ,z2 ,z3 ,z4 ,z5 ,z6 ,z7 ,z8 ,
147 . ajc1 ,ajc2 ,ajc3 ,
148 . ajc4 ,ajc5 ,ajc6 ,
149 . ajc7 ,ajc8 ,ajc9 ,
150 . hx, hy, hz,
151 . pxc1, pxc2, pxc3, pxc4,
152 . pyc1, pyc2, pyc3, pyc4,
153 . pzc1, pzc2, pzc3, pzc4,
154 . gbuf%JAC_I,index,nch )
155 nptr = elbuf_tab%NPTR
156 npts = elbuf_tab%NPTS
157 nptt = elbuf_tab%NPTT
158 llpij = elbuf_tab%BUFLY(1)%L_PIJ
159C
160 DO ir=1,nptr
161 DO is=1,npts
162 DO it=1,nptt
163C-----------
164 lbuf => elbuf_tab%BUFLY(1)%LBUF(ir,is,it)
165c
166 zr = a_gauss(ir,nptr)
167 zs = a_gauss(is,npts)
168 zt = a_gauss(it,nptt)
169 wt = w_gauss(it,nptt)
170 ip = ir + ( (is-1) + (it-1)*npts )*nptr
171 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*wt
172C
173C cas LBUF%L_PIJ=24 global system w/o assumed strain for Isolid=17,18 only
174 CALL s8zjac_i3(
175 . zr,zs,zt,wi,
176 . hx, hy, hz,
177 . ajc1,ajc2,ajc3,
178 . ajc4,ajc5,ajc6,
179 . ajc7,ajc8,ajc9,lbuf%JAC_I,llpij,lbuf%PIJ,lbuf%SIG,lbuf%SIGL,
180 . index,nch,nel)
181c
182 ENDDO
183 ENDDO
184 ENDDO
185C
186 RETURN
187C
188 END
s8_is17jac_i
subroutine s8_is17jac_i(elbuf_tab, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, index, nch, nel)
Definition s8_is17jac_i.F:38
s8zjac_ic
subroutine s8zjac_ic(xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, hx, hy, hz, jac_i)
Definition s8zderi3.F:319
s8zjac_i3
subroutine s8zjac_i3(ksi, eta, zeta, wi, hx, hy, hz, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, jac_i, l_pij, pij, nel)
Definition s8zderi3.F:669