OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s10jacob1.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!|| s10jacob1 ../engine/source/elements/solid/solide10/s10jacob1.F
25!||--- called by ------------------------------------------------------
26!|| s10deri3 ../engine/source/elements/solid/solide10/s10deri3.F
27!||====================================================================
28 SUBROUTINE s10jacob1(ALPH,BETA,W,
29 . X1B,X2B,X3B,X4A,X5B,X6B,X7B,X8B,X9B,X10B,X8A,X9A,X10A,
30 . Y1B,Y2B,Y3B,Y4A,Y5B,Y6B,Y7B,Y8B,Y9B,Y10B,Y8A,Y9A,Y10A,
31 . Z1B,Z2B,Z3B,Z4A,Z5B,Z6B,Z7B,Z8B,Z9B,Z10B,Z8A,Z9A,Z10A,
32 . PX1,PX2,PX3,PX4,PX5,PX6,PX7,PX8,PX9,PX10,
33 . PY1,PY2,PY3,PY4,PY5,PY6,PY7,PY8,PY9,PY10,
34 . PZ1,PZ2,PZ3,PZ4,PZ5,PZ6,PZ7,PZ8,PZ9,PZ10,
35 . NX1,NX2,NX3,NX4,NX5,NX6,NX7,NX8,NX9,NX10,
36 . VOL,NNEGA,INDEX,VOLDP)
37C-----------------------------------------------
38C I m p l i c i t T y p e s
39C-----------------------------------------------
40#include "implicit_f.inc"
41C-----------------------------------------------
42C G l o b a l P a r a m e t e r s
43C-----------------------------------------------
44#include "mvsiz_p.inc"
45C-----------------------------------------------
46C D u m m y A r g u m e n t s
47C-----------------------------------------------
48 INTEGER NNEGA,INDEX(*)
49C REAL
50 DOUBLE PRECISION
51 . X1B(MVSIZ),X2B(MVSIZ),X3B(MVSIZ),X4A(MVSIZ),X5B(MVSIZ),
52 . X6B(MVSIZ),X7B(MVSIZ),X8B(MVSIZ),X9B(MVSIZ),X10B(MVSIZ),
53 . X8A(MVSIZ),X9A(MVSIZ),X10A(MVSIZ),
54 . Y1B(MVSIZ),Y2B(MVSIZ),Y3B(MVSIZ),Y4A(MVSIZ),Y5B(MVSIZ),
55 . Y6B(MVSIZ),Y7B(MVSIZ),Y8B(MVSIZ),Y9B(MVSIZ),Y10B(MVSIZ),
56 . y8a(mvsiz),y9a(mvsiz),y10a(mvsiz),
57 . z1b(mvsiz),z2b(mvsiz),z3b(mvsiz),z4a(mvsiz),z5b(mvsiz),
58 . z6b(mvsiz),z7b(mvsiz),z8b(mvsiz),z9b(mvsiz),z10b(mvsiz),
59 . z8a(mvsiz),z9a(mvsiz),z10a(mvsiz),voldp(mvsiz)
60
61 my_real
62 . px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),px5(mvsiz),
63 . px6(mvsiz),px7(mvsiz),px8(mvsiz),px9(mvsiz),px10(mvsiz),
64 . py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),py5(mvsiz),
65 . py6(mvsiz),py7(mvsiz),py8(mvsiz),py9(mvsiz),py10(mvsiz),
66 . pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),pz5(mvsiz),
67 . pz6(mvsiz),pz7(mvsiz),pz8(mvsiz),pz9(mvsiz),pz10(mvsiz),
68 . nx1(mvsiz),nx2(mvsiz),nx3(mvsiz),nx4(mvsiz),nx5(mvsiz),
69 . nx6(mvsiz),nx7(mvsiz),nx8(mvsiz),nx9(mvsiz),nx10(mvsiz),
70 . vol(mvsiz),alph,beta,w
71C-----------------------------------------------
72C L o c a l V a r i a b l e s
73C-----------------------------------------------
74 INTEGER I,J
75 DOUBLE PRECISION
76 . D,A4MB4,A4,B4,BB,AB,B2,
77 . DXDR,DXDS,DXDT,DYDR,DYDS,DYDT,DZDR,DZDS,DZDT
78
79 DOUBLE PRECISION
80 . AA,A4M1(MVSIZ),B4M1(MVSIZ)
81
82 DOUBLE PRECISION
83 . drdx(mvsiz), dsdx(mvsiz), dtdx(mvsiz),
84 . drdy(mvsiz), dsdy(mvsiz), dtdy(mvsiz),
85 . drdz(mvsiz), dsdz(mvsiz), dtdz(mvsiz),
86 . det(mvsiz)
87C-----------------------------------------------
88c AA = (TWO*ALPH - ONE)*ALPH
89c BB = (TWO*BETA - ONE)*BETA
90c B2 = FOUR*BETA*BETA
91c AB = FOUR*ALPH*BETA
92C-----------------------------------------------
93#include "vectorize.inc"
94 DO j=1,nnega
95 i = index(j)
96 aa = x5b(i) + x6b(i) + x7b(i)
97 . - x4a(i) - x8b(i) - x9b(i) - x10b(i)
98 dxdr = x1b(i) +x8a(i) - x6b(i) + aa
99 dxds = x2b(i) +x9a(i) - x7b(i) + aa
100 dxdt = x3b(i) +x10a(i)- x5b(i) + aa
101C
102 aa = y5b(i) + y6b(i) + y7b(i)
103 . - y4a(i) - y8b(i) - y9b(i) - y10b(i)
104 dydr = y1b(i) +y8a(i) - y6b(i) + aa
105 dyds = y2b(i) +y9a(i) - y7b(i) + aa
106 dydt = y3b(i) +y10a(i)- y5b(i) + aa
107C
108 aa = z5b(i) + z6b(i) + z7b(i)
109 . - z4a(i) - z8b(i) - z9b(i) - z10b(i)
110 dzdr = z1b(i) +z8a(i) - z6b(i) + aa
111 dzds = z2b(i) +z9a(i) - z7b(i) + aa
112 dzdt = z3b(i) +z10a(i)- z5b(i) + aa
113C
114 drdx(i)=dyds*dzdt-dzds*dydt
115 dsdx(i)=dydt*dzdr-dzdt*dydr
116 dtdx(i)=dydr*dzds-dzdr*dyds
117C
118 drdy(i)=dzds*dxdt-dxds*dzdt
119 dsdy(i)=dzdt*dxdr-dxdt*dzdr
120 dtdy(i)=dzdr*dxds-dxdr*dzds
121C
122 drdz(i)=dxds*dydt-dyds*dxdt
123 dsdz(i)=dxdt*dydr-dydt*dxdr
124 dtdz(i)=dxdr*dyds-dydr*dxds
125C
126 det(i) = dxdr * drdx(i)
127 . + dydr * drdy(i)
128 . + dzdr * drdz(i)
129 voldp(i) = w * det(i) * one_over_6
130 vol(i) = voldp(i)
131C
132 ENDDO
133C
134#include "vectorize.inc"
135 DO j=1,nnega
136 i = index(j)
137 d = one/max(em30,det(i))
138C
139 a4 = four * alph
140 b4 = four * beta
141 a4m1(i) = a4 - one
142 b4m1(i) = b4 - one
143 b4m1(i) = d * b4m1(i)
144 a4m1(i) = d * a4m1(i)
145 b4 = d * b4
146 a4 = d * a4
147 a4mb4 = a4 - b4
148C
149 px1(i) = b4m1(i) * drdx(i)
150 py1(i) = b4m1(i) * drdy(i)
151 pz1(i) = b4m1(i) * drdz(i)
152C
153 px2(i) = b4m1(i) * dsdx(i)
154 py2(i) = b4m1(i) * dsdy(i)
155 pz2(i) = b4m1(i) * dsdz(i)
156C
157 px3(i) = b4m1(i) * dtdx(i)
158 py3(i) = b4m1(i) * dtdy(i)
159 pz3(i) = b4m1(i) * dtdz(i)
160C
161 px4(i) =-a4m1(i) * (drdx(i)+dsdx(i)+dtdx(i))
162 py4(i) =-a4m1(i) * (drdy(i)+dsdy(i)+dtdy(i))
163 pz4(i) =-a4m1(i) * (drdz(i)+dsdz(i)+dtdz(i))
164C
165 px5(i) = b4 * (drdx(i)+dsdx(i))
166 py5(i) = b4 * (drdy(i)+dsdy(i))
167 pz5(i) = b4 * (drdz(i)+dsdz(i))
168C
169 px6(i) = b4 * (dsdx(i)+dtdx(i))
170 py6(i) = b4 * (dsdy(i)+dtdy(i))
171 pz6(i) = b4 * (dsdz(i)+dtdz(i))
172C
173 px7(i) = b4 * (dtdx(i)+drdx(i))
174 py7(i) = b4 * (dtdy(i)+drdy(i))
175 pz7(i) = b4 * (dtdz(i)+drdz(i))
176C
177 px8(i) = a4mb4 * drdx(i) - px6(i)
178 py8(i) = a4mb4 * drdy(i) - py6(i)
179 pz8(i) = a4mb4 * drdz(i) - pz6(i)
180C
181 px9(i) = a4mb4 * dsdx(i) - px7(i)
182 py9(i) = a4mb4 * dsdy(i) - py7(i)
183 pz9(i) = a4mb4 * dsdz(i) - pz7(i)
184C
185 px10(i)= a4mb4 * dtdx(i) - px5(i)
186 py10(i)= a4mb4 * dtdy(i) - py5(i)
187 pz10(i)= a4mb4 * dtdz(i) - pz5(i)
188C
189 ENDDO
190C
191 RETURN
192 END
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21
s10jacob1
subroutine s10jacob1(alph, beta, w, x1b, x2b, x3b, x4a, x5b, x6b, x7b, x8b, x9b, x10b, x8a, x9a, x10a, y1b, y2b, y3b, y4a, y5b, y6b, y7b, y8b, y9b, y10b, y8a, y9a, y10a, z1b, z2b, z3b, z4a, z5b, z6b, z7b, z8b, z9b, z10b, z8a, z9a, z10a, px1, px2, px3, px4, px5, px6, px7, px8, px9, px10, py1, py2, py3, py4, py5, py6, py7, py8, py9, py10, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pz9, pz10, nx1, nx2, nx3, nx4, nx5, nx6, nx7, nx8, nx9, nx10, vol, nnega, index, voldp)
Definition s10jacob1.F:37