OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
q4kel2.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!|| q4kel2 ../engine/source/elements/solid_2d/quad4/q4kel2.F
25!||--- called by ------------------------------------------------------
26!|| q4ke2 ../engine/source/elements/solid_2d/quad4/q4ke2.F
27!||--- calls -----------------------------------------------------
28!|| q4kelij2 ../engine/source/elements/solid_2d/quad4/q4kelij2.F
29!|| q4kelijs2 ../engine/source/elements/solid_2d/quad4/q4kelijs2.F
30!||====================================================================
31 SUBROUTINE q4kel2(
32 1 PY1, PY2, PY3, PY4,
33 2 PZ1, PZ2, PZ3, PZ4,
34 3 PYC1, PYC2, PZC1, PZC2,
35 4 AY, R22, R23, K11,
36 5 K12, K13, K14, K22,
37 6 K23, K24, K33, K34,
38 7 K44, K11U, K12U, K13U,
39 8 K14U, K22U, K23U, K24U,
40 9 K33U, K34U, K44U, K11L,
41 A K12L, K13L, K14L, K22L,
42 B K23L, K24L, K33L, K34L,
43 C K44L, DD, GG, DG,
44 D G33, IKSUP, ICP, OFF,
45 E NEL, JCVT)
46C-----------------------------------------------
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, INTENT(IN) :: JCVT
61 INTEGER, INTENT(IN) :: NEL
62 INTEGER ICP, IKSUP
63 my_real
64 . PY1(*), PY2(*), PY3(*), PY4(*),
65 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
66 . PYC1(*), PYC2(*), PZC1(*), PZC2(*), AY(*), R22(*), R23(*),
67 . K11(2,2,*), K12(2,2,*), K13(2,2,*), K14(2,2,*), K22(2,2,*),
68 . K23(2,2,*), K24(2,2,*), K33(2,2,*), K34(2,2,*), K44(2,2,*),
69 . K11U(2,2,*),K12U(2,2,*),K13U(2,2,*),K14U(2,2,*),K22U(2,2,*),
70 . K23U(2,2,*),K24U(2,2,*),K33U(2,2,*),K34U(2,2,*),K44U(2,2,*),
71 . K11L(2,2,*),K12L(2,2,*),K13L(2,2,*),K14L(2,2,*),K22L(2,2,*),
72 . K23L(2,2,*),K24L(2,2,*),K33L(2,2,*),K34L(2,2,*),K44L(2,2,*),
73 . DD(3,3,*), GG(*), DG(3,3,*), G33(3,3,*), OFF(*)
74C-----------------------------------------------
75C L o c a l V a r i a b l e s
76C-----------------------------------------------
77 INTEGER I, J, IS
78 my_real
79 . PYC3(MVSIZ), PYC4(MVSIZ), PZC3(MVSIZ), PZC4(MVSIZ)
80C-----------------------------------------------
81C S o u r c e L i n e s
82C-----------------------------------------------
83 DO I=1,nel
84 pyc3(i) = -pyc1(i)
85 pyc4(i) = -pyc2(i)
86 pzc3(i) = -pzc1(i)
87 pzc4(i) = -pzc2(i)
88 ENDDO
89C
90 IF (iksup==0) THEN
91 is = 1
92 CALL q4kelijs2(
93 1 py1, pz1, py1, pz1,
94 2 pyc1, pzc1, pyc1, pzc1,
95 3 ay, r22, r23, dd,
96 4 gg, k11, k11u, k11l,
97 5 is, nel, jcvt)
98 CALL q4kelijs2(
99 1 py2, pz2, py2, pz2,
100 2 pyc2, pzc2, pyc2, pzc2,
101 3 ay, r22, r23, dd,
102 4 gg, k22, k22u, k22l,
103 5 is, nel, jcvt)
104 CALL q4kelijs2(
105 1 py3, pz3, py3, pz3,
106 2 pyc3, pzc3, pyc3, pzc3,
107 3 ay, r22, r23, dd,
108 4 gg, k33, k33u, k33l,
109 5 is, nel, jcvt)
110 CALL q4kelijs2(
111 1 py4, pz4, py4, pz4,
112 2 pyc4, pzc4, pyc4, pzc4,
113 3 ay, r22, r23, dd,
114 4 gg, k44, k44u, k44l,
115 5 is, nel, jcvt)
116 is = 0
117 CALL q4kelijs2(
118 1 py1, pz1, py2, pz2,
119 2 pyc1, pzc1, pyc2, pzc2,
120 3 ay, r22, r23, dd,
121 4 gg, k12, k12u, k12l,
122 5 is, nel, jcvt)
123 CALL q4kelijs2(
124 1 py1, pz1, py3, pz3,
125 2 pyc1, pzc1, pyc3, pzc3,
126 3 ay, r22, r23, dd,
127 4 gg, k13, k13u, k13l,
128 5 is, nel, jcvt)
129 CALL q4kelijs2(
130 1 py1, pz1, py4, pz4,
131 2 pyc1, pzc1, pyc4, pzc4,
132 3 ay, r22, r23, dd,
133 4 gg, k14, k14u, k14l,
134 5 is, nel, jcvt)
135 CALL q4kelijs2(
136 1 py2, pz2, py3, pz3,
137 2 pyc2, pzc2, pyc3, pzc3,
138 3 ay, r22, r23, dd,
139 4 gg, k23, k23u, k23l,
140 5 is, nel, jcvt)
141 CALL q4kelijs2(
142 1 py2, pz2, py4, pz4,
143 2 pyc2, pzc2, pyc4, pzc4,
144 3 ay, r22, r23, dd,
145 4 gg, k24, k24u, k24l,
146 5 is, nel, jcvt)
147 CALL q4kelijs2(
148 1 py3, pz3, py4, pz4,
149 2 pyc3, pzc3, pyc4, pzc4,
150 3 ay, r22, r23, dd,
151 4 gg, k34, k34u, k34l,
152 5 is, nel, jcvt)
153 ELSE
154 is = 1
155 CALL q4kelij2(
156 1 py1, pz1, py1, pz1,
157 2 pyc1, pzc1, pyc1, pzc1,
158 3 ay, r22, r23, dd,
159 4 dg, g33, k11, is,
160 5 nel, jcvt)
161 CALL q4kelij2(
162 1 py2, pz2, py2, pz2,
163 2 pyc2, pzc2, pyc2, pzc2,
164 3 ay, r22, r23, dd,
165 4 dg, g33, k22, is,
166 5 nel, jcvt)
167 CALL q4kelij2(
168 1 py3, pz3, py3, pz3,
169 2 pyc3, pzc3, pyc3, pzc3,
170 3 ay, r22, r23, dd,
171 4 dg, g33, k33, is,
172 5 nel, jcvt)
173 CALL q4kelij2(
174 1 py4, pz4, py4, pz4,
175 2 pyc4, pzc4, pyc4, pzc4,
176 3 ay, r22, r23, dd,
177 4 dg, g33, k44, is,
178 5 nel, jcvt)
179 is = 0
180 CALL q4kelij2(
181 1 py1, pz1, py2, pz2,
182 2 pyc1, pzc1, pyc2, pzc2,
183 3 ay, r22, r23, dd,
184 4 dg, g33, k12, is,
185 5 nel, jcvt)
186 CALL q4kelij2(
187 1 py1, pz1, py3, pz3,
188 2 pyc1, pzc1, pyc3, pzc3,
189 3 ay, r22, r23, dd,
190 4 dg, g33, k13, is,
191 5 nel, jcvt)
192 CALL q4kelij2(
193 1 py1, pz1, py4, pz4,
194 2 pyc1, pzc1, pyc4, pzc4,
195 3 ay, r22, r23, dd,
196 4 dg, g33, k14, is,
197 5 nel, jcvt)
198 CALL q4kelij2(
199 1 py2, pz2, py3, pz3,
200 2 pyc2, pzc2, pyc3, pzc3,
201 3 ay, r22, r23, dd,
202 4 dg, g33, k23, is,
203 5 nel, jcvt)
204 CALL q4kelij2(
205 1 py2, pz2, py4, pz4,
206 2 pyc2, pzc2, pyc4, pzc4,
207 3 ay, r22, r23, dd,
208 4 dg, g33, k24, is,
209 5 nel, jcvt)
210 CALL q4kelij2(
211 1 py3, pz3, py4, pz4,
212 2 pyc3, pzc3, pyc4, pzc4,
213 3 ay, r22, r23, dd,
214 4 dg, g33, k34, is,
215 5 nel, jcvt)
216 ENDIF
217C
218 RETURN
219 END
q4kel2
subroutine q4kel2(py1, py2, py3, py4, pz1, pz2, pz3, pz4, pyc1, pyc2, pzc1, pzc2, ay, r22, r23, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, k11u, k12u, k13u, k14u, k22u, k23u, k24u, k33u, k34u, k44u, k11l, k12l, k13l, k14l, k22l, k23l, k24l, k33l, k34l, k44l, dd, gg, dg, g33, iksup, icp, off, nel, jcvt)
Definition q4kel2.F:46
q4kelij2
subroutine q4kelij2(pyi, pzi, pyj, pzj, pyci, pzci, pycj, pzcj, ay, r22, r23, dd, dg, g33, kij, is, nel, jcvt)
Definition q4kelij2.F:34
q4kelijs2
subroutine q4kelijs2(pyi, pzi, pyj, pzj, pyci, pzci, pycj, pzcj, ay, r22, r23, dd, gg, kij, kiju, kijl, is, nel, jcvt)
Definition q4kelijs2.F:34