OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i3cst3.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!|| i3cst3 ../engine/source/interfaces/inter3d/i3cst3.F
25!||--- called by ------------------------------------------------------
26!|| i5ke3 ../engine/source/interfaces/inter3d/i5ke3.F
27!|| i6main ../engine/source/interfaces/inter3d/i6main.F
28!|| i9main3 ../engine/source/interfaces/int09/i9main3.F
29!|| intti12a ../engine/source/interfaces/interf/intti12.F
30!|| intvo3 ../engine/source/interfaces/inter3d/intvo3.F
31!|| ud_int5 ../engine/source/implicit/upd_glob_k.F
32!||====================================================================
33 SUBROUTINE i3cst3(
34 1 X1, X2, X3, X4,
35 2 XI, Y1, Y2, Y3,
36 3 Y4, YI, Z1, Z2,
37 4 Z3, Z4, ZI, XFACE,
38 5 N1, N2, N3, SSC,
39 6 TTC, X0, Y0, Z0,
40 7 XX1, XX2, XX3, XX4,
41 8 YY1, YY2, YY3, YY4,
42 9 ZZ1, ZZ2, ZZ3, ZZ4,
43 A XI1, XI2, XI3, XI4,
44 B YI1, YI2, YI3, YI4,
45 C ZI1, ZI2, ZI3, ZI4,
46 D XN1, XN2, XN3, XN4,
47 E YN1, YN2, YN3, YN4,
48 F ZN1, ZN2, ZN3, ZN4,
49 G AREA, LFT, LLT)
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54C-----------------------------------------------
55C G l o b a l P a r a m e t e r s
56C-----------------------------------------------
57#include "mvsiz_p.inc"
58C-----------------------------------------------
59C C o m m o n B l o c k s
60C-----------------------------------------------
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 INTEGER, INTENT(INOUT) :: LFT
65 INTEGER, INTENT(INOUT) :: LLT
66 my_real, DIMENSION(MVSIZ), INTENT(IN) :: X1,X2,X3,X4,XI
67 my_real, DIMENSION(MVSIZ), INTENT(IN) :: Y1,Y2,Y3,Y4,YI
68 my_real, DIMENSION(MVSIZ), INTENT(IN) :: Z1,Z2,Z3,Z4,ZI
69 my_real, DIMENSION(MVSIZ), INTENT(IN) :: XFACE
70 my_real, DIMENSION(MVSIZ), INTENT(OUT) :: SSC,TTC,AREA
71 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: N1,N2,N3
72 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: X0,Y0,Z0
73 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: XX1,XX2,XX3,XX4
74 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: YY1,YY2,YY3,YY4
75 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: ZZ1,ZZ2,ZZ3,ZZ4
76 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: XI1,XI2,XI3,XI4
77 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: YI1,YI2,YI3,YI4
78 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: ZI1,ZI2,ZI3,ZI4
79 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: XN1,XN2,XN3,XN4
80 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: yn1,yn2,yn3,yn4
81 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: zn1,zn2,zn3,zn4
82C-----------------------------------------------
83C L o c a l V a r i a b l e s
84C-----------------------------------------------
85 INTEGER I
86 my_real
87 . a12(mvsiz), a23(mvsiz), a34(mvsiz), a41(mvsiz),
88 . b12(mvsiz), b23(mvsiz), b34(mvsiz), b41(mvsiz),
89 . ab1(mvsiz), ab2(mvsiz), an
90C-----------------------------------------------
91 DO 100 i=lft,llt
92 x0(i) = fourth*(x1(i)+x2(i)+x3(i)+x4(i))
93 y0(i) = fourth*(y1(i)+y2(i)+y3(i)+y4(i))
94 z0(i) = fourth*(z1(i)+z2(i)+z3(i)+z4(i))
95C
96 xx1(i) = x1(i)-x0(i)
97 xx2(i) = x2(i)-x0(i)
98 xx3(i) = x3(i)-x0(i)
99 xx4(i) = x4(i)-x0(i)
100 yy1(i) = y1(i)-y0(i)
101 yy2(i) = y2(i)-y0(i)
102 yy3(i) = y3(i)-y0(i)
103 yy4(i) = y4(i)-y0(i)
104 zz1(i) = z1(i)-z0(i)
105 zz2(i) = z2(i)-z0(i)
106 zz3(i) = z3(i)-z0(i)
107 zz4(i) = z4(i)-z0(i)
108C
109 xi1(i) = x1(i)-xi(i)
110 xi2(i) = x2(i)-xi(i)
111 xi3(i) = x3(i)-xi(i)
112 xi4(i) = x4(i)-xi(i)
113 yi1(i) = y1(i)-yi(i)
114 yi2(i) = y2(i)-yi(i)
115 yi3(i) = y3(i)-yi(i)
116 yi4(i) = y4(i)-yi(i)
117 zi1(i) = z1(i)-zi(i)
118 zi2(i) = z2(i)-zi(i)
119 zi3(i) = z3(i)-zi(i)
120 zi4(i) = z4(i)-zi(i)
121 100 CONTINUE
122C
123 DO 120 i=lft,llt
124 xn1(i) = yy1(i)*zz2(i) - yy2(i)*zz1(i)
125 yn1(i) = zz1(i)*xx2(i) - zz2(i)*xx1(i)
126 zn1(i) = xx1(i)*yy2(i) - xx2(i)*yy1(i)
127 n1(i)=xn1(i)
128 n2(i)=yn1(i)
129 n3(i)=zn1(i)
130 120 CONTINUE
131C
132 DO 140 i=lft,llt
133 xn2(i) = yy2(i)*zz3(i) - yy3(i)*zz2(i)
134 yn2(i) = zz2(i)*xx3(i) - zz3(i)*xx2(i)
135 zn2(i) = xx2(i)*yy3(i) - xx3(i)*yy2(i)
136 n1(i)=n1(i)+xn2(i)
137 n2(i)=n2(i)+yn2(i)
138 n3(i)=n3(i)+zn2(i)
139 140 CONTINUE
140C
141 DO 160 i=lft,llt
142 xn3(i) = yy3(i)*zz4(i) - yy4(i)*zz3(i)
143 yn3(i) = zz3(i)*xx4(i) - zz4(i)*xx3(i)
144 zn3(i) = xx3(i)*yy4(i) - xx4(i)*yy3(i)
145 n1(i)=n1(i)+xn3(i)
146 n2(i)=n2(i)+yn3(i)
147 n3(i)=n3(i)+zn3(i)
148 160 CONTINUE
149C
150 DO 180 i=lft,llt
151 xn4(i) = yy4(i)*zz1(i) - yy1(i)*zz4(i)
152 yn4(i) = zz4(i)*xx1(i) - zz1(i)*xx4(i)
153 zn4(i) = xx4(i)*yy1(i) - xx1(i)*yy4(i)
154 n1(i)=n1(i)+xn4(i)
155 n2(i)=n2(i)+yn4(i)
156 n3(i)=n3(i)+zn4(i)
157 180 CONTINUE
158C
159 DO 200 i=lft,llt
160 an= max(em20,sqrt(n1(i)*n1(i)+n2(i)*n2(i)+n3(i)*n3(i)))
161 n1(i)=n1(i)/an
162 n2(i)=n2(i)/an
163 n3(i)=n3(i)/an
164 area(i)=half*an
165 200 CONTINUE
166C
167 DO 210 i=lft,llt
168 x0(i)=(n1(i)*xn1(i)+n2(i)*yn1(i)+n3(i)*zn1(i))
169 y0(i)=(n1(i)*xn2(i)+n2(i)*yn2(i)+n3(i)*zn2(i))
170 z0(i)=(n1(i)*xn3(i)+n2(i)*yn3(i)+n3(i)*zn3(i))
171 xx1(i)=(n1(i)*xn4(i)+n2(i)*yn4(i)+n3(i)*zn4(i))
172 210 CONTINUE
173C
174 DO 220 i=lft,llt
175 xn1(i) = yi1(i)*zi2(i) - yi2(i)*zi1(i)
176 yn1(i) = zi1(i)*xi2(i) - zi2(i)*xi1(i)
177 zn1(i) = xi1(i)*yi2(i) - xi2(i)*yi1(i)
178 yy1(i)=(n1(i)*xn1(i)+n2(i)*yn1(i)+n3(i)*zn1(i))
179 220 CONTINUE
180C
181 DO 240 i=lft,llt
182 xn2(i) = yi2(i)*zi3(i) - yi3(i)*zi2(i)
183 yn2(i) = zi2(i)*xi3(i) - zi3(i)*xi2(i)
184 zn2(i) = xi2(i)*yi3(i) - xi3(i)*yi2(i)
185 zz1(i)=(n1(i)*xn2(i)+n2(i)*yn2(i)+n3(i)*zn2(i))
186 240 CONTINUE
187C
188 DO 260 i=lft,llt
189 xn3(i) = yi3(i)*zi4(i) - yi4(i)*zi3(i)
190 yn3(i) = zi3(i)*xi4(i) - zi4(i)*xi3(i)
191 zn3(i) = xi3(i)*yi4(i) - xi4(i)*yi3(i)
192 xx2(i)=(n1(i)*xn3(i)+n2(i)*yn3(i)+n3(i)*zn3(i))
193 260 CONTINUE
194C
195 DO 280 i=lft,llt
196 xn4(i) = yi4(i)*zi1(i) - yi1(i)*zi4(i)
197 yn4(i) = zi4(i)*xi1(i) - zi1(i)*xi4(i)
198 zn4(i) = xi4(i)*yi1(i) - xi1(i)*yi4(i)
199 yy2(i)=(n1(i)*xn4(i)+n2(i)*yn4(i)+n3(i)*zn4(i))
200 280 CONTINUE
201C
202 DO 300 i=lft,llt
203 zz2(i)=y0(i)*yy2(i)
204 xx3(i)=zz1(i)*xx1(i)
205 300 CONTINUE
206C
207 DO 320 i=lft,llt
208 IF(xface(i)==zero)GOTO 320
209 IF(zz2(i)+xx3(i)/=zero)THEN
210 ssc(i)=(zz2(i)-xx3(i))/(zz2(i)+xx3(i))
211 ELSE
212 ssc(i)=zero
213 ENDIF
214 IF(z0(i)/=zero)THEN
215 zz2(i)=yy1(i)*z0(i)
216 xx3(i)=xx2(i)*x0(i)
217 IF(zz2(i)+xx3(i)/=zero)THEN
218 ttc(i)=(zz2(i)-xx3(i))/(zz2(i)+xx3(i))
219 ELSE
220 ttc(i)=zero
221 ENDIF
222 ELSE
223 ttc(i)=(yy1(i)-x0(i))/x0(i)
224 ENDIF
225 320 CONTINUE
226C
227 RETURN
228 END
i3cst3
subroutine i3cst3(x1, x2, x3, x4, xi, y1, y2, y3, y4, yi, z1, z2, z3, z4, zi, xface, n1, n2, n3, ssc, ttc, x0, y0, z0, xx1, xx2, xx3, xx4, yy1, yy2, yy3, yy4, zz1, zz2, zz3, zz4, xi1, xi2, xi3, xi4, yi1, yi2, yi3, yi4, zi1, zi2, zi3, zi4, xn1, xn2, xn3, xn4, yn1, yn2, yn3, yn4, zn1, zn2, zn3, zn4, area, lft, llt)
Definition i3cst3.F:50
max
#define max(a, b)
Definition macros.h:21