OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8bicp11.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!|| s8bicp11 ../engine/source/elements/solid/solide8e/s8bicp11.F
25!||--- called by ------------------------------------------------------
26!|| s8ederi_2 ../engine/source/elements/solid/solide8e/s8ederi_2.F
27!||====================================================================
28 SUBROUTINE s8bicp11(
29 1 PX1, PX2, PX3, PX4,
30 2 PX5, PX6, PX7, PX8,
31 3 PY1, PY2, PY3, PY4,
32 4 PY5, PY6, PY7, PY8,
33 5 PZ1, PZ2, PZ3, PZ4,
34 6 PZ5, PZ6, PZ7, PZ8,
35 7 PXC1, PXC2, PXC3, PXC4,
36 8 PYC1, PYC2, PYC3, PYC4,
37 9 PZC1, PZC2, PZC3, PZC4,
38 A BXY1, BXY2, BXY3, BXY4,
39 B BXY5, BXY6, BXY7, BXY8,
40 C BYX1, BYX2, BYX3, BYX4,
41 D BYX5, BYX6, BYX7, BYX8,
42 E BXZ1, BXZ2, BXZ3, BXZ4,
43 F BXZ5, BXZ6, BXZ7, BXZ8,
44 G BZX1, BZX2, BZX3, BZX4,
45 H BZX5, BZX6, BZX7, BZX8,
46 I BYZ1, BYZ2, BYZ3, BYZ4,
47 J BYZ5, BYZ6, BYZ7, BYZ8,
48 K BZY1, BZY2, BZY3, BZY4,
49 L BZY5, BZY6, BZY7, BZY8,
50 M BXX1, BXX2, BXX3, BXX4,
51 N BXX5, BXX6, BXX7, BXX8,
52 O BYY1, BYY2, BYY3, BYY4,
53 P BYY5, BYY6, BYY7, BYY8,
54 Q BZZ1, BZZ2, BZZ3, BZZ4,
55 R BZZ5, BZZ6, BZZ7, BZZ8,
56 S NU, NEL)
57C-----------------------------------------------
58C I m p l i c i t T y p e s
59C-----------------------------------------------
60#include "implicit_f.inc"
61C-----------------------------------------------
62C C o m m o n B l o c k s
63C-----------------------------------------------
64C-----------------------------------------------
65C D u m m y A r g u m e n t s
66C-----------------------------------------------
67 INTEGER, INTENT(IN) :: NEL
68C REAL
69 my_real
70 . PX1(*), PX2(*), PX3(*), PX4(*),
71 . PX5(*), PX6(*), PX7(*), PX8(*),
72 . PY1(*), PY2(*), PY3(*), PY4(*),
73 . PY5(*), PY6(*), PY7(*), PY8(*),
74 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
75 . PZ5(*), PZ6(*), PZ7(*), PZ8(*),
76 . PXC1(*), PXC2(*), PXC3(*), PXC4(*),
77 . PYC1(*), PYC2(*), PYC3(*), PYC4(*),
78 . PZC1(*), PZC2(*), PZC3(*), PZC4(*),
79 . BXY1(*),BXY2(*),BXY3(*),BXY4(*),
80 . BXY5(*),BXY6(*),BXY7(*),BXY8(*),
81 . BYX1(*),BYX2(*),BYX3(*),BYX4(*),
82 . BYX5(*),BYX6(*),BYX7(*),BYX8(*),
83 . BXZ1(*),BXZ2(*),BXZ3(*),BXZ4(*),
84 . BXZ5(*),BXZ6(*),BXZ7(*),BXZ8(*),
85 . BZX1(*),BZX2(*),BZX3(*),BZX4(*),
86 . BZX5(*),BZX6(*),BZX7(*),BZX8(*),
87 . BYZ1(*),BYZ2(*),BYZ3(*),BYZ4(*),
88 . BYZ5(*),BYZ6(*),BYZ7(*),BYZ8(*),
89 . BZY1(*),BZY2(*),BZY3(*),BZY4(*),
90 . BZY5(*),BZY6(*),BZY7(*),BZY8(*),
91 . BXX1(*),BXX2(*),BXX3(*),BXX4(*),
92 . BXX5(*),BXX6(*),BXX7(*),BXX8(*),
93 . BYY1(*),BYY2(*),BYY3(*),BYY4(*),
94 . BYY5(*),BYY6(*),BYY7(*),BYY8(*),
95 . bzz1(*),bzz2(*),bzz3(*),bzz4(*),
96 . bzz5(*),bzz6(*),bzz7(*),bzz8(*),nu(*)
97C-----------------------------------------------
98C L o c a l V a r i a b l e s
99C-----------------------------------------------
100 INTEGER I, J
101C REAL
102C 12
103 my_real
104 . BXHI,BYHI,BZHI,TEMP
105C-----------for the imcompressible laws
106C-------|BXX BYX BZX |
107C-------|BXY BYY BZY | same than H8A
108C-------|BXZ BYZ BZZ |
109 TEMP=half
110 DO i=1,nel
111c TEMP=NU(I)
112 bxhi=temp*(px1(i)-pxc1(i))
113 byhi=temp*(py1(i)-pyc1(i))
114 bzhi=temp*(pz1(i)-pzc1(i))
115c PX1(I)=PXC1(I)+TWO*BXHI
116c PY1(I)=PYC1(I)+TWO*BYHI
117c PZ1(I)=PZC1(I)+TWO*BZHI
118 bxy1(i)=-bxhi
119 bxz1(i)=-bxhi
120 byx1(i)=-byhi
121 byz1(i)=-byhi
122 bzx1(i)=-bzhi
123 bzy1(i)=-bzhi
124 ENDDO
125 DO i=1,nel
126c TEMP=NU(I)
127 bxhi=temp*(px2(i)-pxc2(i))
128 byhi=temp*(py2(i)-pyc2(i))
129 bzhi=temp*(pz2(i)-pzc2(i))
130c PX2(I)=PXC2(I)+TWO*BXHI
131c PY2(I)=PYC2(I)+TWO*BYHI
132c PZ2(I)=PZC2(I)+TWO*BZHI
133 bxy2(i)=-bxhi
134 bxz2(i)=-bxhi
135 byx2(i)=-byhi
136 byz2(i)=-byhi
137 bzx2(i)=-bzhi
138 bzy2(i)=-bzhi
139 ENDDO
140 DO i=1,nel
141c TEMP=NU(I)
142 bxhi=temp*(px3(i)-pxc3(i))
143 byhi=temp*(py3(i)-pyc3(i))
144 bzhi=temp*(pz3(i)-pzc3(i))
145c PX3(I)=PXC3(I)+TWO*BXHI
146c PY3(I)=PYC3(I)+TWO*BYHI
147c PZ3(I)=PZC3(I)+TWO*BZHI
148 bxy3(i)=-bxhi
149 bxz3(i)=-bxhi
150 byx3(i)=-byhi
151 byz3(i)=-byhi
152 bzx3(i)=-bzhi
153 bzy3(i)=-bzhi
154 ENDDO
155 DO i=1,nel
156c TEMP=NU(I)
157 bxhi=temp*(px4(i)-pxc4(i))
158 byhi=temp*(py4(i)-pyc4(i))
159 bzhi=temp*(pz4(i)-pzc4(i))
160c PX4(I)=PXC4(I)+TWO*BXHI
161c PY4(I)=PYC4(I)+TWO*BYHI
162c PZ4(I)=PZC4(I)+TWO*BZHI
163 bxy4(i)=-bxhi
164 bxz4(i)=-bxhi
165 byx4(i)=-byhi
166 byz4(i)=-byhi
167 bzx4(i)=-bzhi
168 bzy4(i)=-bzhi
169 ENDDO
170 DO i=1,nel
171c TEMP=NU(I)
172 bxhi=temp*(px5(i)+pxc3(i))
173 byhi=temp*(py5(i)+pyc3(i))
174 bzhi=temp*(pz5(i)+pzc3(i))
175c PX5(I)=-PXC3(I)+TWO*BXHI
176c PY5(I)=-PYC3(I)+TWO*BYHI
177c PZ5(I)=-PZC3(I)+TWO*BZHI
178 bxy5(i)=-bxhi
179 bxz5(i)=-bxhi
180 byx5(i)=-byhi
181 byz5(i)=-byhi
182 bzx5(i)=-bzhi
183 bzy5(i)=-bzhi
184 ENDDO
185 DO i=1,nel
186c TEMP=NU(I)
187 bxhi=temp*(px6(i)+pxc4(i))
188 byhi=temp*(py6(i)+pyc4(i))
189 bzhi=temp*(pz6(i)+pzc4(i))
190c PX6(I)=-PXC4(I)+TWO*BXHI
191c PY6(I)=-PYC4(I)+TWO*BYHI
192c PZ6(I)=-PZC4(I)+TWO*BZHI
193 bxy6(i)=-bxhi
194 bxz6(i)=-bxhi
195 byx6(i)=-byhi
196 byz6(i)=-byhi
197 bzx6(i)=-bzhi
198 bzy6(i)=-bzhi
199 ENDDO
200 DO i=1,nel
201c TEMP=NU(I)
202 bxhi=temp*(px7(i)+pxc1(i))
203 byhi=temp*(py7(i)+pyc1(i))
204 bzhi=temp*(pz7(i)+pzc1(i))
205c PX7(I)=-PXC1(I)+TWO*BXHI
206c PY7(I)=-PYC1(I)+TWO*BYHI
207c PZ7(I)=-PZC1(I)+TWO*BZHI
208 bxy7(i)=-bxhi
209 bxz7(i)=-bxhi
210 byx7(i)=-byhi
211 byz7(i)=-byhi
212 bzx7(i)=-bzhi
213 bzy7(i)=-bzhi
214 ENDDO
215 DO i=1,nel
216c TEMP=NU(I)
217 bxhi=temp*(px8(i)+pxc2(i))
218 byhi=temp*(py8(i)+pyc2(i))
219 bzhi=temp*(pz8(i)+pzc2(i))
220c PX8(I)=-PXC2(I)+TWO*BXHI
221c PY8(I)=-PYC2(I)+TWO*BYHI
222c PZ8(I)=-PZC2(I)+TWO*BZHI
223 bxy8(i)=-bxhi
224 bxz8(i)=-bxhi
225 byx8(i)=-byhi
226 byz8(i)=-byhi
227 bzx8(i)=-bzhi
228 bzy8(i)=-bzhi
229 ENDDO
230 DO i=1,nel
231 bxx1(i) = px1(i)
232 byy1(i) = py1(i)
233 bzz1(i) = pz1(i)
234 bxx2(i) = px2(i)
235 byy2(i) = py2(i)
236 bzz2(i) = pz2(i)
237 bxx3(i) = px3(i)
238 byy3(i) = py3(i)
239 bzz3(i) = pz3(i)
240 bxx4(i) = px4(i)
241 byy4(i) = py4(i)
242 bzz4(i) = pz4(i)
243 bxx5(i) = px5(i)
244 byy5(i) = py5(i)
245 bzz5(i) = pz5(i)
246 bxx6(i) = px6(i)
247 byy6(i) = py6(i)
248 bzz6(i) = pz6(i)
249 bxx7(i) = px7(i)
250 byy7(i) = py7(i)
251 bzz7(i) = pz7(i)
252 bxx8(i) = px8(i)
253 byy8(i) = py8(i)
254 bzz8(i) = pz8(i)
255 ENDDO
256C
257 RETURN
258 END
s8bicp11
subroutine s8bicp11(px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, 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, bxx1, bxx2, bxx3, bxx4, bxx5, bxx6, bxx7, bxx8, byy1, byy2, byy3, byy4, byy5, byy6, byy7, byy8, bzz1, bzz2, bzz3, bzz4, bzz5, bzz6, bzz7, bzz8, nu, nel)
Definition s8bicp11.F:57